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Chapter 1 

Introduction 



1.1 Motivations and main results 
1.1.1 Motivations 

Graphene is a newly discovered material, and its discovery [70] in 2004 by the group 
of A. Geim and K. Novoselov elicited an enormous interest in the physical commu- 
nitjQ. Consisting in a single monoatomic layer of graphite, graphene can be seen as 
the first realization of a two- dimensional crystal; its stability is a remarkable and 
not yet completely understood issue, since, according to Mermin- Wagner theorem, 
two-dimensional systems are not expected to exist in nature. Before its experimen- 
tal discovery, graphene was known to theorists as an interesting academical problem; 
the first research paper on graphene dates back to the work of Wallace in 1947, |81| . 
where graphene was studied as the "building brick" of graphite. Some of the very 
peculiar features of graphene were brought to the attention of the physical com- 
munity many years later, see [76l EH EH HSl 00] for instance; already at the time 
of these works it was understood that a similar physical compound, if existing in 
nature, would exibit unique physical properties. 

From the pioneering experimental work of Novoselov, Geim et al, |70| 
graphene has shown immediatly a rich variety of interesting electronic and struc- 
tural properties. For instance, the hallmark of graphene is the shape of its electronic 
dispersion relation, which vanishes linearly at the Fermi surface, given by only two 
inequivalent points at neutrality (i.e. when there is only one electron per site in 
average, that is at half-filling); see Fig. 11.11 for a sketch of the energy dispersion 
relation in graphene. As expected from the early works, |81] . the conical shape of 
the low-energy excitations in graphene is a consequence of the hexagonal geometry 
of its two-dimensional lattice, see Fig. 11.21 Such property implies that the low 
energy excitations of graphene are effectively described by massless Dirac fermions, 
with velocity v much smaller than the speed of light (experimentally, v ~ c/300); 
this fact stimulated a fruitful cross-fertilization between different areas of Physics, 
namely high energy physics and condensed matter. For instance, in the work of Se- 
menoff [76] graphene has been proposed as a condensed matter realization of 2 -|- 1 
dimensional Quantum Electrodynamics (QED); and since then, a lot of theoretical 

^Geim and Novoselov have been awarded the 2010 Nobel prize for Physics, "for groundbreaking 
experiments regarding the two-dimensional material graphene". 
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Figure 1.1. Sketch of the energy dispersion relation of graphcne; the plot reports the energy 
£{k) as function of the momentum k. The upper/lower manifolds correspond respectively 
to the conduction and valence bands; the bands meet at six points, among which only two 
are independent (the others lie outside the first Brillouin zone) , and are called Fermi points 
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Figure 1.2. Sketch of the honeycomb lattice of graphene; it can be seen as the superposition 
of two shifted triangular sublattices, namely (black dots) and (red dots). On each 
site sits a carbon atom, which is sp^-bonded to its nearest neighbours; the length of the 
carbon-carbon bonds is about 0.142 nm. 
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research in graphene has been motivated by fundamental questions in quantum field 
theory, e.g. the understanding of quantum anomalies |76| SB], or more recently of 
electron fractionalization in two dimensional systems, [501 \ST\ EH [52] . 

Other very interesting physical properties of graphene are for instance the pres- 
ence of a minimal value of the electric conductivity |69| , of the order of the conduc- 
tance quantum a ~ e'^/h (the exact prefactor is still debated), reached in presence 
of a vanishing density of charge carries; or the insensitivity to localization effects 
which are usually induced by disorder, a phenomenon which was already guessed 
by Fradkin in [25\ [26] . Graphene also shows an unusual quantum Hall effect: the 
spacing between the Hall plateaus is Ae^/h [MIIES], i.e. bigger than for the usual 
quantum Hall effect, and the steps occur at half integer multiples of this value. 
Finally, graphene is considered by many as the ideal canditate for interesting tech- 
nological applications [33], due to presence of extremely high electric and thermal 
condictivities (higher than copper and silver, respectively), and of the overall high 
electronic quality (the electrons have very long mean free paths); moreover, it is 
substantially stronger than steel and at the same time very stretchable. 

Aim of this Thesis is to contribute to the understanding of the effect of the 
electron-electron interactions in graphene; this is a highly debated topic, also in view 
of the recent experimental realization of suspended graphene samples, |14| [23l [57] , a 
playground in which electron-electron interactions - so far obscured by the presence 
of substrates - are expected to play a major role. This is confirmed by recent 
measurements of fractional quantum Hall energy gaps, [231 El l35] . 

We can consider two types of interactions: short range and long range interac- 
tions, corresponding respectively to screened or unscreened Coulomb interactions; 
in this Thesis we shall consider the latter, however let us brefly mention what it is 
known for the former. It is widely believed that the presence of weak short range 
interactions shoud not affect too much the low energy properties of electrons on 
the honeycomb lattice; this was argued on the basis of one-loop computations, see 
for instance [ID] and references therein. Remarkably, this belief has been recently 
rigorously proved by Giuliani and Mastropietro in [361 EZ]; they considered the 
Hubbard model on the honeycomb lattice and proved that the free energy and the 
Schwinger functions are analytic functions of the coupling constant uniformly in the 
temperature and in the system size, provided the coupling is small enough. From 
the mathematical physics point of view, this is one of the two cases in which the 
ground state of a two dimensional interacting fermionic system can be rigorously 
constructed; the other one is |24| . where two dimensional fermions with highly asym- 
metric Fermi surface were considered. 

Regarding unscreened Coulomb interactions, a general picture of their effect on 
the low energy physics of graphene is still lacking; however, there is a wide consensus 
that they should play a nontrivial role. For instance, it has been first proposed by 
Gonzalez, Guinea and Vozmediano, [401 HZl HI] that the Fermi velocity is strongly 
renormalized by the interaction; in particular, it is believed that the presence of 
unscreened Coulomb interaction between the electrons produces an unbounded log- 
arithmic growth of the effective Fermi velocity close to the Fermi points. This 
behavior has been claimed by many authors, see |66t WE[ EB] for instance; all these 
arguments heavily rely on lowest order perturbation theory, and are affected by 
crude approximations: for instance, the presence of the lattice is always neglected 
and a momentum cutoff is imposed by hand. Attempt to include higher order correc- 
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tions have been done in [M], and the same logarithmic behavior has been found at 
two-loops in perturbation theory. The effect of such a strong renormalization of the 
Fermi velocity on physical observables is still puzzling and not yet fully understood; 
for instance, in [M] it was claimed that because of this phenomenon graphene in 
the vacuum is an insulator, a claim which was subsequently proven to be wrong, 
[751118]. 

A quite unsatisfactory feature of all the above results is that they are regulariza- 
tion dependent; in fact, because of the absence of the lattice an ultraviolet regular- 
ization scheme has to be introduced. This led to an explicit dependence of physical 
quantities on the cutoffs, and to controversies in the literature: see [671ll8l[53], where 
three different expressions for the optical conductivity of graphene were found using 
momentum |67| US] or dimensional (53] regularization. And no one of these results 
succeeds in explaining why experimentally [68] the optical conductivity of graphene 
appears to be universal, that is independent of electron-electron interactions and 
lattice parameters. 

Another interesting issue is whether interaction-dependent anomalous exponents 
are present in the many body correlations of graphene. The presence of anomalous 
scaling was first argued in [12] on the basis of lowest order perturbation theory, 
and considering a continuum model in presence of dimensional regularization; in 
fact, in |42j a log-correction to the two point correlation function was found, and 
this was interpreted as the lowest order contribution to an anomalous power law 
decay. The presence of anomalous scaling in the correlations would suggest that the 
system behaves as a Luttinger liquid [IS] , a quite common feature of one dimensional 
systems, which has never been established in two dimensions. A proof of this would 
be of great physical interest, for instance in view of the role that two dimensional 
Luttinger liquids are expected to play in high temperature superconductivity, as 
argued by Anderson [2]. Recently, |56t 175] . on the basis of "large N expansions", 
where N is the number of fermionic "species", anomalous exponents in the scaling 
of various physical quantities were found; but still, the contact with actual graphene 
is not clear, since in graphene = 4. 

Large N expansions are widely used tools in perturbative analyses of quantum 
field theory models; in particular, in the eighties they allowed to argue the emer- 
gence of spontaneous chiral symmetry breaking in 2 -|- 1 dimensional QED, [3l 0]. 
Recently, the same mechanism has been revisited and applied to graphene models 
[5^ [59l H3l [55l [56] , in order to predict the occurrence of a metal-insulator transi- 
tion parametrized by the strength of the many body coupling; this issue is of great 
interest for possible technological applications. The presence of a mass gap has also 
been claimed on the basis of Monte Carlo calculations, see [201 121|. But, again, from 
a theoretical point of view it is not clear how to control the error involved when 
choosing A^ = 4 in the large A^ expansion. Finally, another debated point is the ef- 
fect of the strong renormalization of the Fermi velocity in the gap generation; it has 
been recently proposed in [73] that the growth of the velocity due to the presence of 
unscreened Coulomb interactions depresses the phenomenon of gap opening, while 
in [77] it is claimed that Coulomb interactions support the gap once it is open. 

Therefore, at present time no unambiguous prediction on the effect of electron- 
electron interactions in graphene can be made. The main way to overcome these 
difficulties is to "start from scratch", that is from the definiton of the model that 
one wish to study. In fact, the typical model which has been considered so far to 
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describe Coulomb interactions in graphene has the following Hamiltonian: 

n = n^ + n^, (frnn 

Ti^ = free Hamiltonian of 2 + 1 dim. massless Dirac fermions with velocity v, 
= three dimensional instantaneous Coulomb interaction; 

therefore, in the honeycomb structure of the crystalline lattice is completely 

lost; moreover, the coupling with the electromagnetic field is not complete, since the 
vector potential does not enter in %. Usually, this omission is justified by arguing 
that the contribution from the vector potential should be negligible in view of the 
fact that the ratio v/c\s small; however, this argument cannot be consistent at all 
energy scales, since in this kind of models the interacting Fermi velocity as function 
of the momentum of the quasi-particles blows up at the Fermi surface. Moreover, the 
omission of the vector potential breaks gauge invariance, which has to be respected 
in order to correctly define the coupling with the electromagnetic field. 

1.1.2 Main results 

In this Thesis we shall introduce a new model for graphene in presence of electro- 
magnetic interactions, which takes into account the honeycomb lattice, and where 
the coupling with the electromagnetic field is gauge symmetry preserving; in the 
language of high energy physics, the model that we shall consider is a "lattice 
gauge theory model". Then, we will study the model using rigorous Renormaliza- 
tion Group (RG) methods, which have been developed in the last 25 years by the 
Roma school of Constructive Quantum Field Theory, see [6l |3ll |6T] for extensive 
reviews; this approach is based on the functional Renormalization Group developed 
in the 1980's starting from [281 1291 \7T\ . The methods have been originally devised 
to prove the ultraviolet stability of various quantum field theories, like ip^ theory in 
four dimensions, see |28l 1291 [27] : subsequently, the same ideas have been extended 
to develop a general and rigorous RG approach to interacting fermionic systems, [S]. 
Since then, these techniques have been proved very effective to study many inter- 
acting quantum physical models; in particular, they allowed to rigorously prove the 
Luttinger liquid behavior of various one dimensional systems, see [71 [151 [HI |63] for i^" 
stance, and more recently to investigate the ground state and low energy properties 
of two dimensional systems [IHl [13 [Ml [371 EE] • 

Our results, contained in [SS], can be briefly and informally summarized as 
follows; see Section [1.41 for a more extensive discussion. 

• Anomalous scaling of the two point Schwinger function. We show that the 
scaling properties of the "dressed propagator" is dramatically changed by the 
interaction; in particular, interaction dependent anomalous exponents appear, 
as usual in Luttinger liquids. 

• Emergent Lorentz symmetry. The effective Fermi velocity, that is the one 
appearing in the leading contribution to the two point Schwinger function at 
momenta close to the Fermi surface, tends to the speed of light. 

• Enhancement of the "excitonic" and "charge density wave" response functions. 
We compute various response functions, and we show that their space-time 
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decays are governed by anomalous exponents. In particular, we find that the 
response functions associated to special lattice distortions, the Kekule ones (see 
Fig. II. or to a periodic charge inequivalence on the sites of the honeycomb 
lattice, a charge density wave, are amplified by the interaction. 

• Enhancement of Kekule distortions. We show that small Kekule distortions 
of the honeycomb lattice are dramatically enhanced by the electron-electron 
interactions. 

• A mechanism for spontaneous lattice distortion. Considering the lattice dis- 
tortion as a dynamical variable, we show that the energy of the system is 
extremized in correspondence of a Kekule distortion; the amplitude of the dis- 
tortion is determined by a self-consistence non-BCS equation, from which we 
argue that strong electromagnetic interactions favor the emergence of a spon- 
taneous Kekule distortion and the opening of a gap in the fermionic energy 
spectrum. 

Our results are true at all orders in renormalized perturbation theory; in par- 
ticular, we provide explicit ('W!") bounds on the coefficient of the series of all the 
quantities that we compute, which are uniform in the temperature and in the sys- 
tem size. These bounds are not enough to prove absolute convergence of the series; 
however, here we shall prove that perturbation theory is consistent to all orders, in 
the sense of [5] . 

A key role in the derivation of the above results is played by gauge invariance and 
lattice Ward identities (WI); in particular, following a strategy similar to the one 
introduced in [T2] to prove Luttinger liquid behavior in one dimensional systems, the 
Wis allow us to prove that the Beta function of the effective charge is asymptotically 
vanishing, which means that the renormalized charge is close to the bare one. 



1.2 A honeycomb lattice gauge theory 

Let us start by defining the model that we introduced and studied; our goal is 
to describe electrons hopping on the hexagonal lattice interacting through a quan- 
tized three dimensional electromagnetic field. For this purpose, the Hamiltonian we 
consider has the following structure, see [39] : 

where: (i) TH!^^ is the gauge invariant hopping term; (ii) "H^ contains the contri- 
bution due to the instantaneous Coulomb interaction between the electrons on the 
hexagonal lattice; (iii) H.^ is the energy of the free photon field. Choosing units 
such that h = c = 1, these three terms are given by: 

^gA (T=t4- 
j=l,2,3 

•= y E {nx-'^)f{x-y,0){n^^-l) , 

xGAaUAb 
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where: 



1. A.A = A is a periodic triangular lattice, defined as A = B/LB, where L G N and 
B is the triangular lattice with basis vectors di = ^(3, \/3), 0.2 = ^(3, — \/3); 
Af^ is the area of A. We denote by A* the dual of A; its basis vectors are 
hi = f (1, V3), 62 = f (1, -V3). 

2. The vectors := (1,0), 82 := ^(—1, \/3), 5^ := ^(—1, — \/3) connect each site 
X to its three nearest neighbours. 

3. The operators a^^, 6^ are real space fermionic creation/annihilation op- 

erators, acting respectively on A^, A^ := A^ + 5i and satisfying periodic 
boundary conditions in x G A a', they satisfy the anticommutation relations 

Setting V = (u, W3) with iT G M^, the operators c^^ are momentum space 
bosonic creation/annihilation operators, defined on Vl '■= Vl^Ips = ^f^, n S 
Z} with Vl '■= {p = T^bi + : ^ G Z^}, and satisfying the commutation 
relations 

[Cp,r>Cp,,.,] = -L^A^p,p"^r,r' , [Cp,r,(y,r'] = • ([1214) 

4. j4(^ j) := Jq ds 5i ■ A^^_^^g Qy and is the three dimensional quantized elec- 
tromagnetic vector potential in the Coulomb gauge and in presence of an 
ultraviolet and an infrared cutoff: 



-y AaL ^ \\ 2\p 



X[h*fi]{\p\) 



^k,r ' ^k,r' ~ ^r,r' j S.k,r ' k. — 



X[h*,o]{'t) '■= xii) ~ x{^~^ t) with h* G and x{t) is a smooth compact 
support function equal to 1 for t < and equal to for t > a^M, with 
M > 1 and ag constant to be chosen below. The presence of a finite h* plays 
the role on an infrared cutoff, and it will be removed in the computation of 
physical quantities; instead, the ultraviolet cutoff will be kept fixed. Definition 
(|1.2I5|) implies that ^(^.^s) is periodic in x according to the periodicity of the 
triangular lattice A^, and in X3 along the line [0, L]. 

5. (p{y) is a regularized version of the Coulomb potential, 

■^A^ ~ \P\ 



Definition (|1.216p implies that X3) is periodic in x according to the peri- 
odicity of the triangular lattice A^, and in X3 along the line [0, L]. 
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6. is the density operator, and it is given by J2a '^t a^x a J^a^^ a^x 

depending on whether x £ A a ot x £ Ab- 

Notice that the Hamiltonian (|1.2I1|) is invariant under the gauge transformation 

this is a property which has to be fulfilled in order to correctly describe the inter- 
action of matter with the electromagnetic field. As far as we know, this is the first 
time that a lattice gauge theory model is considered in the analysis of condensed 
matter systems; usually, lattice gauge theory models are introduced in quantum 
field theory as gauge symmetry preserving regularizations of continuum models, to 
be recovered in the limit in which the bond length is sent to zero. Here the bond 
length is kept fixed, and we shall be interested in the infrared, i.e. large distance, 
properties of the model. Moreover, notice also that the Hamiltonian is particle-hole 
symmetric, that is it is invariant under the exchange 

this invariance implies in particular that, if we define the average density of the 
system to be 

where := y^^f^ (at aZ -\- b't p bZ p ) is the total particle number operator, 
then one has p = 1 for all |A| and /3. This is the so-called half-filling condition. 



1.3 The free theory 

Before discussing the properties of the model in presence of interaction, that is for 
e 7^ 0, we briefly review what it is known in the case e = 0. In absence of interaction 
the fermionic and bosonic degrees of freedom are completely decoupled, and we shall 
be interested in the fermionic sector only. The Hamiltonian of the system is given 
by: 

xsA (T=t4- 

i=l,2,3 

the ground state of this model can be explicitly determined, in the sense that the 
free energy and the n-point Schwinger functions can be explicitly computed. In 
fact, being the Hamiltonian quadratic in the fermionic creation/annihilation oper- 
ators, this last fact simply follows from the knowledgle of the two point Schwinger 
functions and from the Wick rule. Let us call ^^o-p' P ~ Fourier trans- 

form of the imaginary time evolutions of , 5^ p , namely e^*^°ai q-'^a^o ^nd 

•^1^ x-\-Oi,a -^1^ 

g'HA2-'05±^_, g-'^AX'o. jj- — ^j^^^ jZ^ g 'P/3,L := x'Z^L, where 2?^ is the set of fermionic 
Matsubara frequencies and Vl is the first Brillouin zone, that is 

Vj3 := |fco = — (no -h -) : ?io G Z| , 

Vl := {a?= + ^62 : n G < 7ii,n2 < L- l} , ^11) 
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where bi = ^(1, \/3), 62 := 3^(1, -\/3) form a basis of the dual lattice A*. We 
define the non-interacting two point Schwinger function in momentum space, or free 
propagator, as 

where denotes the grand-canonical average with respect to the Hamiltonian 

Ti/^. An explicit computation shows that, see [SB] or Appendix IB.l. 11 

where n{k) = ^f^^ e'^^^^^""^!); its modulus \^{k)\ is the energy dispersion relation, 
and it is plotted in Fig. Il.li Interestingly, the function Q{k) is vanishing if and only 
if A; = p^, where fy^ are the two Fermi points, 

j_ /2n 27r \ , , 

the fact that the energy dispersion relation vanishes only at two points is very 
unusual for two dimensional systems, where the Fermi surface typically consists 
of a closed curve. Close to the Fermi points n{k' + p^) = |(iM ± + O(I^T); 
therefore, setting p'p = (0,pp) and k' = k — p'^, we can rewrite 

Q ^^.| ^r.^\ - LI ~^^o -v{-\k[ + ujk'2) + r^{k')\ 



(II2115) 

where Z = 1 is the bare wave function renormalization, u = |t is the bare Fermi 
velocity, and < C|A;'p for some C > 0. Therefore, from (|1.2115p we see that 

the free propagator is asymptotically the same of massless Dirac fermions in 2 -|- 1 
dimensions. 



I. 4 Results and discussion 

Below we shall discuss more extensively our results, which have been presented in 
|39] , obtained with methods similar to those of [S5] ; as already mentioned in Section 

II. l.2| they concern the effect of the electromagnetic interaction on the two point 
Schwinger function, various response functions, the effect of lattice distortions and 
a possible mechanism for gap generation. A sketch of the proof together with the 
main ideas underlying the methods is given in Section 11.51 In what follows we will 
assume that the infrared limit h* — )• —00 has been taken. 

1.4.1 The two point Schwinger function 

We start by discussing the result on the two point Schwinger function; as we are 
going to see, the interaction dramatically modifies its scaling properties. In fact, 
the interacting two point Schwinger function is given by, for L — >■ -|-oo and close 
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to the singularities p'^.: 
5(k' + p^) = 

-1 / iko v{k'){-ik[ + ujk'^) + f^{k' 

Z{k') [v{k'){ik[+ojk'2)+fZik') iko 

with \r{k')\ < C|A?'|2, R(k') = 0{e^) and 

Z(k')=^|k'r\ ^;(k')~l-(l-?;)|k'|'^, (11312) 

where v := |t, r/(>2) = O(e^), f/(>^) = O(e^); -R(k'), r/(>^) and fy^^^) are expressed 
as renormahzed series with finite coefficients admitting A^! bounds, in the sense of 
Theorem 13.3.11 The functions Z(k') and u(k') are cahed respectively the wave func- 
tion renormalization and the effective Fermi velocity; as we see from (|1.4l3p . both 
these objects have a non-trivial dependence on the quasi-momentum k': they scale 
with interaction dependent anomalous exponents. 

This results tells us that the scaling of the two point Schwinger function is dra- 
matically modified by the presence of the electromagnetic interaction; in particular, 
in absence of interaction 5*0 (k') ~ |k'|~^, see (|1.2[15p . while in the interacting case 
^(k') ~ |k'|-i+'' with T] = 0{e^) > 0, therefore 

k'^o S{k' + p%) 

This is due to the fact that the quasi-particle weight Z(k')~^ vanishes as a power 
law at the singularity; this suggests that the interacting system is a Luttinger liquid, 
in the sense of |45| . If convergence of the series is proved, this would be the first 
rigorous proof of Luttinger liquid behavior in more than one spatial dimension. 
Regarding one dimensional systems, these have been widely studied in the last 15 
years, starting from the work of Benfatto, Gallavotti, Procacci and Scoppola in |7], 
where the first rigorous proof of Luttinger liquid behavior for a non solvable one 
dimensional system was given; in particular, quantities of physical interest like the 
Schwinger functions, the anomalous exponents and the free energy were explicitly 
computed as convergent series in the coupling constant. A crucial ingredient in the 
proof of [7j is the vanishing of the Beta function of the effective coupling; in [7] 
this was shown using informations coming from the exact solution of the Luttinger 
model. Recently, an independent strategy to prove this remarkable cancellation has 
been proposed and developed by Benfatto and Mastropietro, [15] , which completely 
avoids the use of the exact solution; this new method is based on the rigorous 
implementation of Ward identities in the RG, and similar ideas will be adopted 
here. 

Another remarkable feature of the two point Schwinger function is the non- 
trivial behavior of the effective Fermi velocity f(k'); in fact, in the limit |k'| — >■ 
the effective Fermi velocity tends to 1, which is the speed of light in our units, for 
any value of the bare Fermi velocity v. 
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Finally, we stress that our result does not neglect the lattice, and does not 
assume unphysical regularization procedures, like the dimensional one. The first 
attempt to investigate the effect of electromagnetic interactions in graphene is due 
to Gonzalez, Guinea and Vozmediano in [30]; the model they considered described 
2+1 dimensional Dirac fermions interacting with a 3 + 1 dimensional photon field 
in the Feynman gauge; both the fields lived in the continuum, and dimensional 
regularization was used to avoid ultraviolet divergences. One-loop computations 
suggested that logarithmic corrections to the scaling of the two point function were 
present, and these corrections were interpreted as the first order contributions to 
the expansions of |k'|~'', |k'|''. For the same model, anomalous scaling has been 
recently established at all orders by Giuliani, Mastropietro and Porta in [38], using 
momentum regularization (in order to mimic the presence of the lattice) instead of 
the dimensional one. 

Remark 1 In general, we shall represent the grand- canonical average (•)^ us- 
ing the functional integral representation in the Feynman gauge, see Section \2.2\: 
if the averaged observable is gauge invariant the result of the computation is equal 
to the corresponding grand- canonical average with the Hamiltonian T-Lj^, see Appen- 
dicesOO All the quantities that we shall compute in this Thesis, except the two 
point Schwinger function, will be gauge invariant; however, the two point Schwinger 
function appears as "dressed propagator" in the perturbative series of physical ob- 
servables, and therefore it is interesting it its own right. 

1.4.2 The response functions 

The same methods used to evaluate the two point Schwinger function can be used to 
compute other correlations; in particular, the analysis of various response functions 
allows us to investigate the effect of the electromagnetic interactions on possible 
quantum instabilities which may take place at strong coupling. Let us define: 

^ x+{0,<5,),cr x+(0,<5,).o- ' 

cr=t4- 

EaX rrOZ „ — b'^ , p , b~ , -. ^ , 
X,(T x,cr x+(0.(5,),cr x+(0,(5,).o- ' 

cr=U 

where a^^g., 6y g., ^(x,i) are respectively the imaginary time evolutions under the 
Hamiltonian T-L\ of a^^, 6^^, A(^j); these correlations are called the excitonic, 
charge density wave or density- density susceptibilities, depending on whether a = 
E±, a = CDW or a = D. The correlations C^-^^^ measure the tendency of the 
system to form particle-hole pairs between nearest neighbours on the honeycomb 
lattice, and in particular C^^+^ measures the response of the system to Kekule lat- 
tice distortions, see next section; C^'^^^^ measures the tendency of the system to 
produce charge asymmetries between the two triangular sublattices; finally, C^^^ 
measure the correlations of on-site electronic densities. 



C!?(x-y) 

(CDW) 



(D) 
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It follows that, for (3,L ^ +00 and |x| » 1, the Cfi correlations are given bjH 

Q"i(x) = + cos(p^ • ^)7§S- + r-?,i(x) , (II216) 

|x| '^++ |x| 

where: (i) the functions Gf (x) are given by, if 6 = 8tt'^/27: 

0^2 I I ^1 2 ^2 I ^2 ^2 

Gf ^^(x) = , G^^'^(x) - ""^0 + 



b ' ^ ^ ' 6|x 



2 



Gf- (X) = + , G^ (X) = ; C217) 



(ii) the anomalous exponents are given by renormalized series in the renor- 
malized charge, and are equal to: 

C = + C'^>'^ , V(a,^) : {a,oj) + +-) , (CW, ++) , 

where ^S'^"^^^ are expressed as renormalized series starting from fourth order, ad- 
mitting A^! bounds in the sense of Theorem I3.3.H 

(iii) rf]^(x) contains corrections which are either bounded by Ixl"^^^^"™*^^— for 
^ < < 1 or by e |x| ^ i^'^i^/^ 

The above results show that the electromagnetic interaction has the effect of de- 
pressing the space-time decay of the correlations C^^+\ c(C!DW).^ ^j^jg g^gggg^g tj^at 
the interaction may favor excitonic or charge density wave instabilities at strong 
coupling. The depression of the decay of the oscillating part of the a = sus- 
ceptibilities is associated with a particular lattice distortion, the Kekule one; as we 
are going to see with the next result, the amplitude of Kekule distortions is greatly 
enhancedhy the interactions. The relevance of charge density wave and Kekule insta- 
bilities in graphene have been first discussed in [SU [59] and [501 ED El] , respectively. 
As far as we know, this is the first time that these correlations are computed for 
graphene; regarding one dimensional systems, analogous computations have been 
performed in [H [8]. 

1.4.3 Lattice distortions 

In this Section we shall discuss the effect of the interaction on some special lattice 
distortions, the Kekule ones. First, we will report the computation of the two point 

^The C"j(x) correlations are obtained from C°i(x) using that C"^ij^.i(x) = C°',j{xo,Tx). 
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Schwinger function in presence of Kekule distortion; then, we will discuss a possible 
mechanism for spontaneous lattice distortion in presence of strong enough electron- 
electron interactions. Remarkably, the amplitude of the Kekule distortion behaves 
as a bare mass for the fermion propagator; therefore, the problem of spontaneous 
lattice distortion is equivalent to the one of mass generation in quantum field theory. 

The two point function in presence of Kekule distortion 

The presence of a small distortion can be taken into account by assuming that the 
hopping parameter t is replaced by a lattice site and bond dependent one; namely 
by replacing t in (ll.2lip with t + AQ(j)g^i, where is the distortion of the bond 
connecting x to x + 5i, and Aq is a small parameter. This approximation has been 
widely used to study the Peierls instability in one dimensional systems, see [SDl 155] . 
for instance. In this framework, the Kekule distortion is obtained by replacing t in 
(jl.2i2|) with t-|-|Ao cos{pp{5j — Sj^ —x)); this choice corresponds to the dimerization 
pattern represented in Fig. 11.31 




Figure 1.3. The Kekule dimerization pattern for jq = 1. Double and single bonds are 

3 ' 3^ 



respectively "long" and "short", and correspond to the hopping parameters t+ IAq. 



To express the result for the two point Schwinger function is convenient to adopt 
a "relativistic" notation; let 

^l,T ._ { 

'A'k'.a •- \^ k'+p+,(7,l k'+p+,(T,2 k'+p^,<T,2 k'+p-,cr,l 

j: .— ^+ , , 

H^]i',U ■- \ k'+P~,(T,2 k'+p-,(7,l k'+P+,(T,l k'+p+,CT,2 

where ^^o-p' P = li^, are the Fourier transform of the imaginary time evolutions 
of , 6^ ^ , and define 

It follows that, for /3, L — t- +00, and |k'| small, the limiting propagator is given by: 
5A(k') = -^ -—4 , , , (1 + i?(kO) , (EllO) 
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where: (i) 70,71,72 are Euclidean gamma matrices (defined in (|3.2l3ip . (|2.4[18|) ). 
and 

(ii) setting 

A 

(>2) 

where r]j^ is given by a series in the renormahzed charge starting from the fourth 
order and admitting A^! bounds, the functions A(k'), v(k'), Z(k') are given by: 

Z(k') ~max{|k'|,A}"'' , v{k') ~ 1 - (1 - t;) max{|k'|, A}^ , 

A(k') ~ Aomax{|k'|, A}-''^^- ; ^13) 

(iii) R{k') = Ri{k') + i?2(k') with Ri{k') = 0{e^) and i?2(k') = 0(|k'|). 

Therefore, the presence of a Kekule distortion with a non-vanishing amphtude 
Aq produces an effective mass A(k') for the fermion field, which is strongly renor- 
malized by the interaction; in fact, 

lim A = . (0114) 

Ao^O Ao 

In this sense we say that the presence of electromagnetic interactions enhances the 
Kekule distortion of the honeycomb lattice; as far as we know, this is the first time 
that the effect of electromagnetic interaction on a preexisting Kekule distortion is 
discussed. 



-i/e'PF(5.o-^i) A 

i/e'PF('^Jo~'^i) / 



A mechanism for spontaneous lattice distortion 

Finally, we conclude the summary of the results of this Thesis by discussing a 
mechanism for spontaneous lattice distortions. We consider the lattice distortion 
as a dynamical variable; the full Hamiltonian of the model is: 

^a(W)=^a(W)+/Ca({</)}):=?^a(W) + ^ ^ 4>%, (mis) 

1=1,2,3 

where: 'Ha({(^}) is given by (jl.2iip after the replacement t ^ t + (p^^i] the second 
term is the elastic energy of the distortion and k, g are respectively the stiffness 
constant and the phonon coupling. Let us define the specific free energy in the 
Born-Oppenheimer approximation as: 



g-/3/CA({0})-/3|A|F;5,i. ({</,}) ^ ([13116) 
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where is the specific free energy of the model with Hamiltonian ^{/^{{(j)}) 

in presence of a fixed distortion. It follows that 9,^^ .E'^j^ ({</)}) = in correspon- 
dence of a Kekule distortion 

4:?"^ = -^0 + ^Ao cos{pp{6j - 4 - x)) (11118) 

with (j)Q and Aq self-consistently determined by the following equations: 

3 K |A| Pf'I l 

where: the symbol (• • • )^^^ denotes the limit j3 — t- +00 of the grand-canonical average 
(■ ■ ■ )/3 L computed with the Hamiltonian 'Ha({0*''-^H); the operator P'^^'^ is the 

( F \ 

Fourier transform in x of p^^g) j defined in (|1.415|) . Moreover, for g small enough 
the solution (I1.4ll8p . ()1.4ll9p corresponds to a local minimum of the energy. 

The first equation gives 00 = 0{g^); this amounts to a small renormalization 
of the hopping parameter. The second equation is the most interesting one; it can 
be seen as a non-BCS gap equation for the fermion field. To understand it from a 
qualitative viewpoint, we replace it with: 

Ao 8g^ f dk' 1 A(k') „ , 

" ^ I ' ([1120) 



3 K J D Z(k') fc^ + 'L.(k')2|k'|2 ' 

A<|k'|<l 

where D = (27r)|i3i| and Bi = 87r^/{3^/3) is the volume of the first Brillouin zone; 
Eq. (|1.4i2U|) is very similar to the gap equation first derived by Mastropietro in [02] , 
in the context of certain Luttinger superconductors. It follows that: 

• for small e, that is when ijk — V — + • • • is negligible with respect to 1, 
the equation ()1.4l20p admits the solution 

^0 ^ (1 - ^) , for \g\ >go = 0{Vkv) , ([1121) 

where g^ is positive and essentially independent of e; moreover, since the 
exponent in ()1.4l2ip is greater than 1, dgA{g) is continuous at l^l = go: the 
transition is smoothened by the interaction. 

• For large e, that is when rjK — rj ^ 1, then g^ — )• 0. In particular, if ry;^ — 77 — 1 
exceeds 0, then 

Ao ^ (1 + ^^n^i^l^) ^ , for \g\ > 0. ([1122) 

Putting these results together with those regarding the excitonic susceptibility 
and the two point Schwinger function in presence of Kekule distortion, we 
argue that strong enough electron-electron interactions may favor the occurrence 
of a possible "excitonic quantum instability" at strong coupling, in the form of a 
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spontaneous Kekule distortion of the honeycomb lattice. This issue is interesting not 
only because of a possible interaction driven metal-insulator transition in graphene 
(the amplitude of the distortion produces a gap in the energy spectrum of the 
fermions); but also because, as proposed first by Hou, Chamon and Mudry in [50], 
and by Jackiw and Pi in , the Kekule distortion may be a prerequisite for electron 
fractionalization, that is for the existence of collective excitations with net charge 
equal to a half the electric charge. This would be the two dimensional analogous 
of the Peierls instability studied in one dimensional systems by Su, Schrieffer and 
Heeger in |80] and by Jackiw and Rebbi in |52j . 

1.5 Sketch of the proof 

Here we shall give a sketch of the proof of the above results, trying to explain in 
an informal way the ideas behind our RG methods. Our framework consists in a 
rigorous formulation of the Wilsonian Renormalization Group, |83l , based on 
the functional RG developed in the 1980's in [251 [Ml E] , see also [27]. In the last 
25 years these methods have been successfully applied to different fields, such as 
quantum field theory, condensed matter and classical statistical mechanics. See 
[B] |3S1 IH] for extensive reviews. 

The underlying physical idea is very simple: we look for a sequence of effective 
theories that describe the system at larger and larger scales, whose parameter are 
determined by the integration of the degrees of freedom corresponding to smaller 
scales. Each theory is determined starting from the ones on smaller scales by ap- 
plying a well defined procedure, the "RG map"; the goal is to find a fixed point of 
this map, and to determine the scaling behavior of the corresponding many-body 
correlations. 

First of all, we rewrite our many-body problem in terms of a quantum field 
theory model, in the sense that we rewrite the partition function and the generating 
functional of the correlations as suitable functional integrals, involving fermionic, i.e. 
Grassmann, and bosonic, i.e. real, Gaussian fields. The equivalence is established 
noting that the perturbative series are equal order by order; see Appendix IbI where 
the case of the partition function is discussed. Our goal is to evaluate such functional 
integrals, at least in terms of perturbative series with coefficients bounded uniformly 
in /3, L, which may or may not converge (the problem of proving convergence will 
be not adressed in this Thesis). Let us consider for simplicity the partition function 
'^13, L', the generating functional of the correlations can be studied in a similar way, 
just more technically involved, and we will not discuss it here. It follows that: 



where V{^,A) is the interaction and P{d^), P{dA) are Gaussian fermionic and 
bosonic measures respectively, with propagators ^(k), w{p); the interaction and 
the gaussian measures are defined in Section 12.21 The identity has to be understood 
in the sense that the perturbative series in e of the l.h.s. and of the r.h.s. are equal 
order by order for any /3, L, h* fixed; moreover, by gauge invariance, see Appendix 
IC] we are free to express the photon propagator in the Feynman gauge, which is 
more convenient for our purposes. The momenta k = {ko, k), p = {po,p) appearing 
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at the argument of the propagators are constrained by the fact that: (i) k has to 
belong to the first Brillouin zone; (ii) \p\ cannot be too big, because of the bosonic 
ultraviolet cutoff in the definition of the model, which will be kept fixed. However, 
the "Matsubara frequency" variables /cq, Po are unbounded; therefore, to have a well 
defined problem we fix an ultraviolet cutoff on both of them, by imposing that the 
arguments of the propagators must satisfy /cq < , pQ < , with M > 1 a 
scaling parameter and K a positive integer, which at the end will be sent to +oo. 

One may naively think of evaluating (jl.Slip by expanding the exponential in 
(|1.5iip and taking averages; as it is well known, the result can be expressed in 
terms of a sum over Feynman graphs. However, an easy computation shows that, 
already at lowest order, these graphs are not bounded uniformly in /?, L, h*: this 
is so because the fermionic and bosonic propagator are singular at k = and 
p = 0. The hope is that these infinities cancel and produce a finite result. Roughly 
speaking, to see these cancellations we shall write each graph as a sum of many 
pieces, where in every piece each propagator carries a momentum which is "close" 
to some prefixed value, labelled by some momentum scale label. The goal is to show 
that classes of finite pieces cancel, and produce a result which is summable over the 
scale labels. 

To see these cancellations directly from the perturbative series would be an 
almost desperate task; the startegy that we shall adopt to reorganize in a convenient 
way the series and to keep track of all the contributions has been introduced by 
Gallavotti in 1984, see [27]. First of all, setting k' = k — p'^, we rewrite the 
propagators as 

m = J2 5i-°HkO +5(>°)(k) , ^(p) = u;(^o)(p) +zS(>o)(p) , (11312) 

where g!^^\'k'), w^-^\p) are supported respectively on momenta close to p'^. and 
0; correspondingly, we rewrite the fermionic and bosonic fields as 

UJ=± 

where the fields labelled by (< 0) and (> 0) are independent Gaussian variables, 
with propagators given by (|1.5[2p . After this, we integrate the fields labelled by 
(> 0) and we get a theory which depends only on the (< 0) fields; this is done by 
writing 

VW(., .) = log J P>o(d^(>°))P>o((i^^>°^) [l + V{- + ^(>°), • + A(>°)) + . . 

where the integration measures labelled by < 0, > have propagators given by 
(|1.5[2p . and the dots stand for the higher order terms coming from the expansion 
of the exponential in (jl.Sllj) . The expectations in the second line of (I1.514P can 
be graphically represented in terms of Gallavotti - Nicold (GN) trees, [281 EHl [27] . 
which in turn can be evaluated as sums over connected Feynman graphs, with an 
arbitrary number of external lines denoting the fields labelled by (< 0); the values of 
these graphs, and therefore of the trees, are bounded uniformly in K, see Appendix 
[D1 Therefore, the limit K — t- +oo can be safely taken. 
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So far we discussed the "easy" part of the work; now comes the "hard" one, since 
singularities k' = 0, p = belong to the domains of the fields labelled by (< 0). 
We rewrite: 

h=—oo h=—oo 

where ^W, ^("i), A^-^\ ... are independent Gaussian fields supported on 
momenta |k'| ~ M^, |p| ~ {h is the scale label), whose propagators will be de- 
termined inductively. We integrate the fields iteratively expanding the exponential 
as in (ll.5[4p . starting from scale going down to — cxd; each step of the integra- 
tion can be organized graphically in terms of GN trees. It follows that, after the 
integration of the first \h\ scales, the partition function (|1.5I4|) can be rewritten as: 

Ef,,L = e-'^l^l^'^ / P<,(^ivI/(^^))P<,(d^(^'^))e^"^Hv^*^-''^'^'-''') , (II316) 
where: 

• Ffi is the specific free energy on scale h. 

• The integration measures have propagators given by g!^'^\k'), w'^-^^p), both 
supported on momenta < M^~^^. For these values of the momenta, the photon 
propagator is essentially equal to w^-'^\p), while ^^-^^(k') is morally equal 
to (I1.2I15P after having relaced Z and v with Zh and Vh, the wave function 
renormalization and the effective Fermi velocity on scale h. 

• The interaction on scale h has the form 

where J^^'* is a "relativistic" fermionic current, e^^/i, M'^Uf^^h are respectively 
the effective charges and the effective mass terms for the photon field, and 
contains irrelevant terms in the RG sense. 

The parameters Zh, Vh, e^,h, M^u^^h are called running coupling constants; it 
follows that, collecting all these numbers in a vector oth, the iterative integration 
implies an "evolution equation" for q/j, namely: 

a/i = a/i+i + /3h+i(a/i+i) • • • )"o) ; ([I3l8) 

the vector in (jl.5l8|) is called the Beta function of the theory, and it can 

be computed as series in the running coupling constant on scales > h + \. It 
turns out that the Beta function, the free energy and the correlation functions 
can be expressed as perturbative series in the renormalized couplings {e^^fc, fp,fc}) 
with coefficients depending on {Z^^/ Z}^_i, Vk}; in particular, the coefficients of the 
series are finite if the ratio Zk/Zk_i is close to 1 and if is bounded away from 
0. Therefore, the problem of computing the partition function (I1.5I6P is reduced 
to the study of a finite dimensional dynamical system, evolving under p.5l8|) : an 
explicit computation shows that if e^^h — e and v^^h — then by truncating (jl.SlSp 
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to any finite order p > 0, clioosing e small enough, we get that — M~'^^ , Vh — 
V-oo — {l—v)M^'^, where rj, fj are O(e^) and positive, and V-oo — 1 = 0(eo,-oo— ei^_oo)- 
Clearly, one would show that the remainder of the truncation is small; however to 
do this one has to prove the convergence of the series, which is notoriously a difficult 
task for bosonic field theories, and it is outside the purposes of this Thesis. However, 
in the language of [5], here we prove that perturbation theory is consistent to the 
order p for any p > 0. 

But then how to check the assumptions on the renormalized couplings? No- 
tice that now we cannot rely on the evolution equation (jl.SiSp . because a single 
non- vanishing order in the components of /?/,, contributing to e^^^h, f^^/^ would gener- 
ically produce an unbounded flow as h — ^ — oo. The informations that we need 
are provided by Ward identities; the Ward identities are non perturbative identi- 
ties between the Schwinger functions of a field theory, which are implied by gauge 
invariance. To get informations on e^^/,, z/^^/i we will implement Ward identities at 
each RG step (that is at each single scale integration), following a strategy recently 
proposed and developed by Benfatto and Mastropietro, [15], in the context of one 
dimensional Luttinger liquidsH The strategy can be summarized as follows. 

• Assume inductively that our assumptions on the renormalized couplings are 
true on scales > h. 

• At the step — /i of the RG we introduce a new model, the reference model, 
which differs from the original model because of the presence of a cutoff on 
the bosons that suppresses all momenta smaller that M^; the idea of putting 
a cutoff only on the bosonic sector of the theory is borrowed from Adler 
and Bardeen, [1]. This new model can be investigated using the multiscale 
integration described above, and it turns out that on scales > h it has the 
same running coupling constants of the full model. After the scale h all the 
bosonic fields have been integrated out, and we are left with a purely fermonic 
theory; but this theory is much easier to investigate with respect to the initial 
one, and in particular it turns out to be superrenormalizable. This implies 
that the running coupling constants essentially "cease to flow" on scale h. 

• The bosonic cutoff does not break gauge invariance, and Ward identities for the 
Schwinger functions of the reference model can be derived; by suitably choos- 
ing the external momenta of the Schwinger functions involved in the Ward 
identities we will be able to check our inductive assumptions on e^i^h, ^n,h- 

Finally, another interesting consequence of Ward identities is that eo,-oo = ^i-oo 
(while ei^h = (^2,h for all h, as it follows from the 27r/3 rotational symmetry of the 
model); therefore, since in general f-_oo = l + 0(eo,-oo — ei^^oo), this last fact implies 

^This strategy allowed to prove the so called vanishing of the Beta function for one dimensional 
fermionic systems. At the formal level, this was known since the 1970's, see [221 119| : however, the 
earlier works neglected the presence of cutoffs, which necessarily break gauge invariance. Therefore, 
the problem of establishing whether gauge invariance and formal Ward identities were recovered 
in the limit of cutoff removal was not considered by the Authors of the original proof. In [T^] 
this problem has been considered, and the Authors proved that actually the Ward identities found 
after the removal of the cutoff are different from the formal ones: this is the phenomenon of chiral 
anomaly, well-known in the context of similar models used in Relativistic Quantum Field Theory, 
e.g. the Schwinger model [86| . 
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that u_oo = 1) which is the speed of Hght in our units. This concludes the sketch 
of the strategy that allows to compute the free energy, the Schwinger functions and 
the correlation functions of the model. 

Now, let us conclude the Section by briefly discussing what happens in presence 
of Kekule distortions; the effect of Kekule distortions can be taken into account 
by replacing the hopping parameter with j = t + |Ao cos{p^{5j — Sj^ — x)) with 
jo = 1, 2, 3. The main difference with respect to the case Aq = is that Aq behaves 
as a hare mass for the fermion propagator, and it gives rise to a new running coupling 
constant A/j, which grows as AqM~^'^^, with rjx > 0. As for the case Aq = we 
perform a multiscale integration to evaluate the partition function, but here we 
have to distinguish two regimes: in the first regime the mass Af^ is smaller than 
M^, which is the size of the momentum of the single scale field ^k/^g-pw' ^^'^ 
analysis is qualitatively the same discussed before. The second regime is defined by 
the scale h(A) such that A/^^^) = M^^^\ where the fermion propagator becomes 
massive; the main consequence of this fact is that the running coupling constants 
on scales < h(A) remain close to their values at the threshold /i(A). 

1.6 Summary 

In Chapter [2] we introduce a functional integral representation for our model and, 
as a pedagogical example, we start to evaluate the free energy by discussing the 
integration of the ultraviolet degrees of freedom; the analysis will be similar, but 
not equal, to the one discussed in [36j for the Hubbard model on the honeycomb 
lattice with short range interactions. The outcome of the ultraviolet integration 
is an effective infrared theory, whose action has a precise structure determined by 
remarkable lattice symmetries. 

In Chapter [3] we continue the evaluation of the free energy by discussing the 
integration of the infrared degrees of freedom; the procedure will be similar to the 
one perfomed in [38], in the context of an effective continuum model for graphene. 
Here we shall define the localization and renormalization operations, which allow 
to safely integrate the infrared scales. In this way we get a sequence of infrared 
effective theories, involving fields depending on momenta closer and closer to the 
singularities and a number of running coupling constants; the various contributions 
to the effective actions will be graphically represented in terms of Gallavotti-Nicolo 
trees. This graphical representation will be crucial to derive infrared-stable bounds 
on the effective potentials of the effective actions. After this, we discuss the flow of 
the running coupling constants, and in particular we introduce the reference model 
on which we derive the Ward identities that we need in order to control the flow of 
the effective couplings. 

In Chapter H] we adapt the strategy developed in the previous two Chapters to 
the generating functionals of the Schwinger functions and of the respose functions; 
here shall we prove our results on the two point Schwinger function (ll.4[ip and 
(|1.4l3p . and on the response functions (|1.4l6p - (|1.4[8p . 

In Chapter Owe change the definition of the model in order to take into account 
Kekule distortions. We first study the effect of the electromagnetic interaction on 
a preexisting Kekule distortion; then, we investigate a mechanism for a possible 
spontaneous Kekule instability. Here we will prove the result on the two point 
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Schwinger function in presence of Kekule distortion (|1.4il0|) - (jl.4il3|) . and derive 
and discuss the gap equation (I1.4I18P - (ll.4[2ip . 

In the Appendices we collect all the technical ingredients and explicit computa- 
tions that are needed in our proofs; in particular, in Appendix [B] we show how the 
Hamiltonian model is mapped in a Quantum Field Theory, in AppendixIClwe prove 
the perturbative equivalence of Coulomb and Feynman gauges, and in Appendix IFl 
we collect all the explicit computations of this Thesis. 



Chapter 2 

The ultraviolet integration 



2.1 Introduction 

In this Chapter and in the next one we introduce our rigorous RG framework, by 
performing as a pedagogical example the computation of the specific free energy. 
The same methods with some modifications will be used in Chapter H] to prove our 
results on the two point Schwinger function and on the response functions, and in 
Chapter [5] to take into account the presence of lattice distortions. Here we start by 
discussing the integration of the ultraviolet degrees of freedom; this is the "easiest" 
part of the work, since the bosons are regularized by an ultraviolet cutoff, which will 
be kept fixed, and the hexagonal lattice provides a natural ultraviolet cutoff for the 
fermions. As an outcome of the integration we will get an effective infrared theory, 
where the various contributions to the effective action will be expressed as series 
in the bare charge e with finite coefficients, the A^-th order bounded proportionally 
to (A/2)!. In particular, as a consequence of some remarkable lattice symmetry 
properties of our model, the effective action will have a precise structure, which will 
be preserved by the subsequent infrared integration. 

The Chapter is organized in the following way: in Section [2.21 we introduce the 
functional integral representation of the model; in Section [2. 31 we discuss the symme- 
tries of our model, while in Section 12.41 we describe the outcome of the integration 
of the ultraviolet regime; the details of the ultraviolet integration are contained in 
Appendix iDl 



2.2 Functional integral representation 

Here we shall introduce the quantum field theory associated to the model defined 
in Section [1.2t in fact, as shown in Appendix [Bl for any L, h* fixed the pertur- 
bative series in e of the partition function of our model, which can be computed 
using Trotter formula and Wick theorem, turns out ot be equal order by order to 
the perturbative series of a suitable functional integral, involving Grassmann (i.e. 
fermionic) and real {i.e. bosonic) Gaussian fields. The same is true for the average 
of physical observables which can be expressed as series in the fermionic and bosonic 
creation/annihilation operators. 
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2.2.1 Integration measures and interaction 

In this Section we shall introduce the elements that define the quantum field theory 
in which we will map our Hamiltonian model; namely, the free fermionic and bosonic 
measures, and the fermion-boson interaction. 



Fermionic gaussian measure. Let K G'N, and 



■■= ^/3,L n {ko : x(M-^|fco|) > 0} . dZZll) 



We consider the Grassmann algebra generated by the Grassmann variables ^^o-; 
with a =ti) /9 = 1) 2 and k E ^, and a Grassmann integration 



p=l,2 



we define the Fourier transform of the fermionic field as: 



I3,L 

where has been defined in (ll.2llip . In the following, we shall use the notation 
/ ^ = (/3|A|)-i Eke©,3 ^^^^^ ^ ■= and Bi = 87rV(3^/3) is the volume 

of the first Brillouin zone. Let us define the free fermionic propagator g{h) as 

=: XK{ko)[B{k)r\ XK{ko):=x{M~''\ko\) . 
Then, we introduce the Grassmann gaussian integration P{d^) as 

•exp{ - (/3|A|)"i ''^^ ^.^,.,i(k)^^^k..,.} ; 

by the discussion of Appendix lAl it turns out that 

/ ^(^^)^k,.,p^k',.',p' = /3|A|5,,<x'4,k' [^(k)]^^^, . me) 

Bosonic gaussian measure. Let 

„ ^ r 2-Km ^^ 

Vp^L ■■= Vl U |po = m G Z| , 

'PIl ■■= n {p = (po,l^ : X[h',o]i\p\)XKiPo) > 0} , a?) 

- U {po = m G Z+} n {p = (po,p) : Xlh*,o]{\p\)XK{po) > 0}, 
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where the set Vl has been defined in the lines before (|1.2i4p . and consider the 
complex Gaussian variables ^^,p with = 0, 1, 2, p G Vp^i and A*^^^ = this 
last condition implies that the Fourier transform of the ^^,p field, namely 



A 



IS 



real. In the following, we shall use the notation / = (/3^a) ^ Yl,peVp l' 
^ G [0, 1] and no = 1, = 0, we define the free photon propagator in the ^-gauge as 

dp3 XK{po)X[h*fi]{\p\) 



(27r) p2+p2 



A« (p,p3) , dog) 



where G {0, 1, 2}, and the gaussian integration P^{dA) as 



310) 



P^dA) 



n h 



pep-t H^/3^A)'detz^?(P/- ^=0,1,2 



13, L 



•exp{-(2MA)-i Yl Ap[*^(p)] '^-.-p}; 
notice that if ^ G [0, 1] then p G V^^^ =^ det #(p) > 0, and 



J P«(dA)i^,pi^,_p/ = MA5p,p,tl;J ,,(p) 



(Haii) 



12) 



Interaction. Finally, we define the interaction V{'^, A) as 



jG[l,3] 

-ie ^ nxAo,x 



313) 



where ^{x,i) := Jo ds5i ■ ^x+s(0(5 )" useful to rewrite (j2.2[13|) in momentum 

space; it follows that 



V{^,A) = 



(fc+p)(<5,-5i) 



-ie / E ['I'k+p,.,l^k,.,l^,P + ^^^+p,.,2^k,.,2^,pe"'^^'^ 



([22114) 
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where k + p is constrained to live in V^^l and 

dxQ 



Fj,p := I dxo 5^ e^P''[exp{ieA(^,j)} - l] ([22115) 

s^/^...|^^(p-Ep.)[n..(«)%-4] 

A'>1 \ I \ I J I— I 



with 77j(p) := — TTTT^: and (5(p) := (^^^1^1 ^)'^(p) with 5(p) := S{po)6{p) where 

5ik):=\A\ Y: %nAW.' m)--=mo,o- (122116) 

2.2.2 Averages 

Consider a generic function 0{'^, A) of the fields A; we shall consider only func- 
tions 0{'^,A) which are integrable in the sense of [S], i-e. that can be expressed as 
series in the fields ^, A, with bounded coefficients. For any positive (3, L, we define 
the average in the .^-gauge of C'(^', A) as 

{0{^,M,L ■■= ^\}^^J^^^JpmPHdAV^''^^MA,^)^ 

Sg;-^! := I P{d^)P^dA)e^^''''^\ ([22117) 

where S^^^ := limi^_>._|_oo lini/i*_j.„oo ^{3*/,'^' is the partition function of the model. 
As explained in Appendix |B1 the quantity ()2. 21170 can be formally expressed as 
a perturbative series in the electric charge e; in particular, if O is the operator 
obtained by suitably replacing the fermionic and bosonic fields in O with fermionic 
and bosonic creation/annihilation operators, it follows that 

^A^P^ = (0(^,A))J,,. ([22118) 

Moreover, as shown in Appendix [Ul 0{'^,A) is gauge invariant, that is if it is left 
unchanged by the transformation '^^,a,p ~^ ^x.o-.p^'^'"'') ^m.x ^m,x + f^^Ox, then: 

d^{O{^,A))l^ = 0. ([22119) 

This means that in the functional integral (|2.2117p we are free to choose the gauge 
that we prefer; for convenience, in this Thesis we shall work in the Feynman gauge, 
corresponding to the choice ^ = 0. We shall set 

P{dA) := P^=\dA) , {... ),r :=(••• )^ . ([22120) 



2.3 Symmetries 

Before discussing the multiscale integration, it is important to note that both the 
Gaussian integrations and the interaction V{'^,A) are invariant under the action 
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of suitable symmetry transformations; this invariance will be preserved by the sub- 
sequent iterative integration procedure, and will guarantee the vanishing of some 
running coupling constants. In the following Lemma we collect all the symmetry 
properties that will be needed in the following. 

Lemma 2.3.1 For any choice of K , h* , j3, L, the fermionic Gaussian integration 
P{d^), the bosonic Gaussian integration in the Feynman gauge P{dA) and the 
interaction V^^, A) are invariant under the following transformations: 

(1) spin exchange: ill^^^ o ^k-a.p/ 

(2) global U{1): ^'k,a,p ~^ ^''^"''^k a.p' G 1^ independent ofk; 

(3) spin S0{2): { f ^'^'^ | ^Re\ X^'^A , with Rg = ( """^^^ ''"^^ | and 9 e T 

—\n,iJ \n,iJ i^-sm^ cos^; 

independent of k, where T is the one- dimensional torus; 

(4) discrete spatial rotations: ^ geiTfc(52-5i)(p-i)^± ^ ^^^^^^ _^ 
(^o,(po,Tp)) ^~^^(po,Tp)); where T is the 27r/3 rotation matrix. 

(5) complex conjugation: 'I'ko-p ~^ '^-ko-p' ^m.p ~^ ~^m.-P' ^ ~^ ^ ^ where c is 
generic constant appearing in P{d^), P{dA) and/or in V{'^,A); 

(6. a) horizontal reflections: ^'^^^ ^^ o ^\k,,^^k^^k2),a,2' ^Lp ^ -e'^'''^i,(po -pi,P2)' 
(6.b) vertical reflections: il^^^^ %oM,-k2),a,p' ^2,p ^ -i2,(po,Pi,-P2)' ^^,p ^ 

^I^,(po,Pl,-P2) ^/'^ / 2. 

(7) particle-hole: ^^^^^ ^ ^^lko,-k),a,p' ^ 'A-po,?)' ^ ^,(-po,p)- 

(8) inversion: ^^^^^ ^ i("l)^^(-fco,fc),<x,p' ^ A-po,p}' ^ -^,(-po,p)- 

Proof. For notational simplicity, here we shall consider the limit if — )■ +00; the case 
of a finite K can be worked out in exactly the same way. 

The invariance of P{d^), P{dA), V{'^,A) under (1), (2), (3) is obvious, and so is 
the invariance of P{dA) under (4)-(8). 

Symmetry (4). Let us prove the invariance of 
k 

= ^k,.,1^0^k,.,l - E ^'k,.,l*^^*(^)^k,.,2 - E ^'k,.,2i^^(^)^k,.,l 

k k k 
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which imphes the invariance of P{d^). The first and the fourth term in the second 
hne of ()2.3[ll) are obviously invariant, while the sum of the second and the third 
term is changed into 

-ty^t n*(k)e'^^^^~^''^i/- ^j., +C.C. (12312) 

k 

k k 

Using that n{T-^k) = e'^^^^-^^'^n{k), as it follows by the definition of n{k), we 
find that the last line of (I2.312P is equal to the sum of the second and third term 
in (|2.311|) . as desired. Let us now prove the invariance of the temporal part of the 
interaction, that is: 

- ^k+p,.,i^k,.,i^o,p - '^k+p,.,2^k,.,2^,pe-^^^'^ • as) 

k,p k,p 

The first term in (j2.313|) is obviously invariant under (4) , while the second is changed 
into 

- E ^.^+p,.,2^^.,2^,pe-^^"^'^e^^^-(^^-^^) ; m^) 
k,p 

using that T5j = (^j+i, the invariance of (|2.313p follows. Finally, let us prove the 
invariance of 

''Z^ ^k+p,cr,l^k,(T,2-^J.P^ ^ ^k+p,(T,2^k,(T,l-^J -P^ 
k,p j k,p j 

(12315) 

Notice that under (4) ^{x,j) — ^ A{{xo,Tx),j+i)^ which means that Fj^p — )• -Fj^i j^p^; 
therefore, (|2.3[5p is changed into 

''Z^Z^ ^k+p,cr,l^k,cr,2-^J + l,P*3 ^ ^ 

k,p j 

E E ^k+p,.,2^k,.,i^+i,-pe^^'+^^^'^-^^-'^^e^^^^+p^(^^-^^^) , me) 

k,p j 

which shows the invariance of ()2.3[6p and concludes the proof of the invariance of 
Vi'^jA) under (4). 

Symmetry (5). The invariance of P{d^) follows simply by noting that Q{—k) = 
Q*(k). Concerning V{'^,A), ()2.3l3p is obviously invariant; the invariance of (I2.3l5p . 
and therefore of V{'^,A), follows by noting that under (5) ^(x,j) — ^ — ^(x,j) and 
hence Fj^p — >■ -Fj,-p. 

Symmetry (6. a). To check the invariance of P{d^) notice that under (6. a) the sum 
of the first and the fourth term in the second line of (|2.3iip is obviously invariant, 
while the sum of the second and the third is changed into 

E ^Jo,-;..fe),.,2^*(^)^r^o,-fci,;..),.,i + ^-^^ 07) 

k 

= -*E^k,.,2^*((-^^i'^2))^'k,.,i -tE^k,.,i^((-^i'^2))^'k,.,2 ; 

k k 
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noting that 0,{{—ki,k2)) = ^*{k), one sees that (|2.317p is equal to the sum of the 
second and third term in the second hne of ()2.3[ip . as desired. The invariance of the 
temporal part of the interaction is obvious. To conclude the proof of the invariance 
of A) notice that F^, (p^^.p^^^^) ^ e'^'-^o' F*, with 1' = 1, 2' = 3, 3' = 2. 
In fact, setting x := (xo,-xi,X2), after (6.a) ( ~ ^i,x+(0A)' ^2,x+(o,5i)) 

and -/o^t^'5'^j'^x-{o,s5*.,-5i)! therefore, setting dx = J dxoJ^neA^ with 

/ dxe'P'^exp {ieA(xj)} e*^''i / dx e^^'' exp { - / ds 6j'A^_,^ 

J A. J B Jo 'J 

= e^^^(^-^^)/^dxe^P-exp{ -ze £ S,,A^^^^^^^__^)s,)} 

= e'P^^' ) dx e*P^ exp { - ieA^^j,) } , (I018) 

where in the first equality we used that = Aa+Sj/ and that p-(0, 5j) = — p-(0, 6j/). 
Formula (|OTHl) . together with the fact that / dxe^P'' = e*P(^i-^j) / tixe^P'' impHes 
that Fj^p e*P(^i-^j)F*, _~; therefore, dOTS]) changes into 

E E ^k+p,.,i^'k,.,2e-*('^+^^^'^'-'^^e^^^(^.'-^^^)F^^ (12319) 

k,p j 

which proves the invariance of V{'^,A) under (6. a). 

Symmetry (6.h). The invariance of P{d'^) follows by noting that r2((A;i, — /C2)) = 
r2(/c). The invariance of the temporal part of V{'^,A) is obvious. To show the 
invariance of (j2.3[5|) , and therefore to conclude the proof of the invariance oiV{'^ ^ A), 
notice that under (6.b) Fj^^ — >■ Fji^p with p := (poiPi)— P2); then, (|2.3[5p changes 
into 

''Z^ ^k+p,cr,l^k,cr,2-^j',P*^ ''Z^ ^k+p,cr,2^k,cr,l-f^j' ,-p6 ' 
k,p j k,p j 

(123110) 

and (|2.3ll0p concludes the proof of the invariance under (6.b). 

Symmetry (7). To show the invariance of P{d'^)^ notice that the sum of the first 
and the fourth term in (j2.3[ip is left unchanged, while the sum of the second and 
the third term is changed into 

^ (ko,-k),a,l ^ ' {ko,-k),a,2 ^ {ko-k),a,2 ^ ' {kQ,-k),a,l 
k k 

= -*E^k,.,2^*(-^)^k,.,i -*E^k,.,i^(-^)^k,.,2 ; Qii) 

k k 

using that Q{—k) = Q*{k), we see that (|2.311ip is transformed into the sum of 
the second and third term of (|2.3[ip . as desired. Regarding the interaction, (|2.3l3p 
changes into, setting k := {ko, —k): 

-ieYi>t , ~ 1 ^ -P - ie y n^r ~ -v^"'^^^ , (123112) 

^ k,(T,l k+p,o-,l ' P ^ k,(T,2 k+p,cr,2 ^' P ' ' ^ 

k,p k,p 
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which after a change of variables gives (j2.3l3|) . To show the invariance of the spatial 
part of V{'^,A) (f23T5]) first notice that under (7) F^-p F*p; therefore, (I23l5]) is 
changed into 

k,p J 
k,p i 

which is equal to (|2.3[5p . 

Symmetry (8). Under (8) all the terms appearing in (|2.3lip are separately invariant; 
this proves the invariance of the same is true for V{^,A). This concludes 

the proof of Lemma 12.3.11 ■ 



2.4 The ultraviolet integration 

We start by studying the partition function ' =: e~ ' ' ''■^ . Note that in 

our model the fermionic momenta k have an intrinsic ultraviolet cut-off induced by 
the lattice, while the bosonic momenta p have an ultraviolet cut-off because of the 
definition of the photon field; on the contrary, the ko,pQ variables are not bounded 
uniformly in K. A preliminary step to our infrared analysis is the integration of the 
ultraviolet degrees of freedom corresponding to large values of ko,pQ. We proceed 
in the following way. We decompose the free propagators ^(k), w)(p) into sums of 
two propagators supported in the regions of ko,po "large" and "small". The regions 
of kQ,pQ large and small are defined in terms of the smooth support function x{t) 
introduced in Section II. 2t the constant ao entering in its definition is chosen so that 

the supports of + 1^ ~ -PfP) ^^"^ + 1^ ~ ^'fP) ^^'^ disjoint (here | • | 

is the Euclidean norm over \ A*). In order for this condition to be satisfied, it is 
enough to choose 2aoM < 47r/(3-v/3); in the following, for reasons that will become 
clearer later (see discussion after Lemma I2.4.ip , we shall assume the slightly more 
restrictive condition uqM < tt/G. We rewrite ^(k) and w{p) as 

5(k) =5("-^-)(k)+5(^-)(k) , «;(p) =7^("-)(p)+^I;(--)(p) (gail) 

where, setting pp = (0,p^) with uj = ziz: 

g^--^-\k) = m - g^'-^-\k) , <7(--)(k) = X(|k - p^|)5(k) , (12312) 

UJ=± 

^(u.v.)( N ^/ N -^(i.r.)( N ^(i.r.)/ N .X[/1*,0](|P|) ^^'^^^^ (ijf) 
'{p) = w{p) — '{p) , '{p) = I— — j-^-j — — , 

IpI ^ 

with w{p) := w^{p) is the photon propagator in the Feynman gauge and / is the 
3x3 identity matrix; notice that arctan (^^f^) = f + 0(|p|)- We now introduce 

four sets of independent Gaussian fields {^k"o"p^}, {^k ^p^^ ^^"^ {^a^p '*}) {^m,p''} 
defined by, if # = n.v., i.r.: 

j P(d^(#))^L*j;^ktSp' = /3|A|<5.,.'5k,k'5^5(k) , 

|p(dA(#))4#]iLt,\, = /3|A|Vp'45(p) • (12313) 
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Similar to P(ci\I'), P{dA), the gaussian integrations P(rf^'(#)), also admit 

an explicit representation analogous to (I2.2l5p . (I2.2[lip . with ^(k), w{p) replaced 
by 5^*''(k), t()(*)(p) respectively, and the sums over k, p restricted to the values in 
the support the corresponding propagator. The definition of Grassmann integration 
implies the following identity ("addition principle"): 

P(d^)P(d^)e^(*'^) = (IM4) 



so that we can rewrite the partition function as 
^[h*,K] ^ f p(rf^(-))p(d^(-0)exp{^lC.(m^""^ + ^^^""^ + 0;^ 

J ^ Tl. 



where the truncated expectation £^y^ is defined, given any series X A^'^''"'^) 
in ^("•''•)^ depending on A^'-''-\ ^'(^•'■■), as 

C.(^(-);n) := ^log / P(d*("-))P(d^(«-))e^^(*^-^-^'^^"-"-^)|^^^ , (IM6) 

and Vk is fixed by the condition Vk{0,0) = 0. It can be shown, see Appendix [Pl 
that Vk can be written as 

Vi^(^,A)= (IM7) 



„ rt m 



:.r.) 

n,m>0 "i=l i=l 

n+m>l 

n X 

,2n,in,p,ii 

({k,},{pa)5(E(k 

j=i i=i 

where / = /(^...^^^...^,p:= (pi, . . . , P2n), g •= (q- i, • • • , cin), := 
(//I, . . . , fj,m), and the symbol 5(k) has been defined in Eq. (j2.2116p : the possibility 
of representing Vk in the form (|2.4i7p . with the kernels independent of the spin 
indeces ai, follows from the symmetry listed in Lemma 12.3.11 As an outcome of the 
discussion reported in Appendix [Dl the constant Fq^k and the kernels WK,2n,m,p,^i 
are given by formal power series in e, 

^0,K - 2^ ^o^K e 



K ^ ' 

Af>0 



T^X,2n,^,p,^({k,}, {p.}) = ^f2n,™,p,^({k,}, {pj) dMO) 



with coefficients bounded uniformly in /?, L, /i* independent of e; in particular, 
the following result holds. 
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Lemma 2.4.1 There exist some positive constants Cn,m, C independent of (3, L, 
K, h* such that 



(TV) 



0,K 



< C 



N 



N 



\e 



max 



([^ 10) 



moreover, the following limits exist: 



lim F: 



(N) 



0,K ' 



aiii) 



(TV) 



K->-+oo 



Proof. See Appendix [Dl 



Regarding the choice of the constant ao in the definition of this is motivated 
as fohows. The quadratic part in the fermionic fields of V^^*'^\'^,A) is 

E / §^ ^(q)^£i^^l:;i;p^ w^x,2,o,p.p'(k, k + q) . mi2) 

From the definitions (|2.4lip . (j2.4l2|) . the support of consists of two disjoint 

regions around and p^; because of the assumption oqM < tt/G if k and k + p 
belong to the support of then |p| < 47r/3, which implies that the only non- 

zero contribution in the sum over p in (|2.4ll2p is the one corresponding to p = 0. 
Analogously, the contributions to Vk with n = l,m = 1 and n = 0,m = 2 are, 
respectively: 

J :^:^(^'^(^)^k+p+q,<.,p^Up' ^M,P^i^,2,i,p,p>(k,k+p+q,p) 

p,p' 

E / ^(q)^i:-i-.i^;>KA2,p,p'(P, -P - q) ; »i3) 

the only non-vanishing contribution in the sums over q are those associated to q = 0. 
Moreover, in the first line of (j2.4[13p . because of the assumption oqM < vr/G, if 
k & = {q : |^— #p| < aoM} and k-|-p is in the support of ^f^*-^-) then k+p G B^^. 



Remark 2 From now, with a little abuse of notation we shall only write the inde- 
pendent momenta at the argument of the kernels. 

Finally, it is important for the forthcoming discussion to note that the symmetries 
listed in Lemma 12.3.11 imply some non-trivial invariance properties of the kernels. 

Lemma 2.4.2 Let Wk,2,o(}^), W^A',2,i,/i(k, p) be the matrices with entries given by 
Wk,2,o,p,p'0^), 1^-R',2,i,p,p',p(k, p), respectively; in the limit f3 — >■ -|-oo, |A| — )• -|-oo the 
following properties are true: 

WkM^' + pf) = -^pk;r- + o(|kf) , dzaiM) 

WK,0,2,^^AP) = V^P + 0(|p|'), (123115) 

M^j^,2,i,p(k' + p'j^,p) = A^r;i + o(|k'| + |p|) , (123116) 

W^ir,0,3,p(Pl,P2) = 0(|P1| + |P2|) , (123117) 
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where: (i) z^, X^, are real and such that zi = Z2, Xi = X2, ^'i = 1^2; (ii) the 
matrices are given by: 

Proof. See Appendix lEl 



Chapter 3 

The infrared integration 



3.1 Introduction 

As an outcome of the ultraviolet integration discussed in the previous Chapter and 
in Appendix [Dl we obtained that the low energy physics of our model is described 
by an infrared effective theory. In this Chapter we analyze this theory by Renormal- 
ization Group; this is the "hardest" part of the work, since the remaining degrees of 
freedom that we are left with can be arbitrarily close to the singularities k = p*^ and 
p = 0. First of all, we will rewrite the fermion field as a sum over two independent 
fields with disjoint supports, each one supported around one of the two singulari- 
ties; these new fields will be called quasi-particle fields, [5]. Then, we rewrite the 
quasi-particle fields and the photon fields as sums over independent fields, depend- 
ing on momenta closer and closer to the infrared singularities; the propagators of 
these fields will be determined inductively. As for the ultraviolet integration, see 
Appendix ini we will integrate scale after scale these new fields, and at each step of 
the integration we will be left with a new effective theory, whose parameters (wave 
function renormalization, effective masses, effective couplings, and effective Fermi 
velocity), also called running coupling constants, are determined by the integration 
of the higher energy degrees of freedom. Roughly speaking, the single step is per- 
formed by splitting the effective potential in the sum of two contributions, the local 
part and the renormalized part; the local part of the effective potential is absorbed 
in the redefinition, i.e. in the renormalization, of the fermionic integration mea- 
sure and in the definition of the effective couplings: this step is the "heart" of the 
KG, and it corresponds to an infinite resummation of Feynman graphs. The new 
fermionic measure is similar to the old one, except for the fact that now the wave 
function renormalization and the effective Fermi velocity appearing in the fermion 
propagator are slightly changed with respect to the ones before the redefinition. 
Then, we rescale the fermionic fields in a suitable way, and we perform the single 
scale integration; the outcome of the integration has the same form of the effective 
potential from which we started, and we can iterate the process. 

The outcome of the integration, that is the new effective potential, can be ex- 
pressed graphically in terms of Gallavotti-Nicold trees, [SSlEn]; each tree has a value 
which can be computed as a sum over renormalized Feynman graphs, and it turns 
out that if the running coupling constants verify suitable bounds then the value of 
each tree is finite. To fully appreciate this fact one should compare with the naive 
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perturbative series, which is plagued by infrared divergences. 

Because of the iterative integration scheme, the running coupHng constants on 
a given scale are determined starting from those on higher scales, and ultimately 
from the bare ones, exploiting a non-trivial recursion relation; in fact, the running 
coupling constants on a given scale can be obtained starting from those on the 
previous scale plus a well-defined function of all the running coupling constants 
on higher scales: this function is called the Beta function. The Beta function is 
expressed as a series in the running coupling constants, with finite coefficients; by 
truncating the Beta function to any finite order in the couplings and choosing the 
bare charge e small enough uniformly in the momentum scale we will be able to 
check our assumptions on the running coupling constants and to close the single 
scale integration. Clearly, one would be able to prove that the remainder of the 
truncation is small; however to prove this one should be able to prove convergence 
of the series, which is notoriously a very hard task in bosonic theories, and it is 
outside the purposes of this Thesis. 

However, the fact that even truncating the Beta function the flow of the running 
coupling constants is bounded is remarkable, and follows from the implementation 
of Ward identities at each integration step; this should be compared with what 
happens for instance in "asymptotically slave" theories {e.g. the 93^ theory in 4 
dimensions, |86| I27|). where the "wrong" sign of the lowest order contribution to 
the Beta function makes the flow of the effective coupling unbounded. Here the 
phenomenon taking place is well known in one dimensional systems, and it is called 
vanishing of the Beta function [3 [HJ [12] . This is crucial to establish Luttinger liquid 
behavior; as far as we know, this is the first time that the same cancellation is found 
in a two dimensional system. 

The Chapter is organized in the following way: in Section 13.21 we describe the 
infrared multiscale analysis; in Section [33] we discuss the graphical representation 
of the effective potentials in terms of Gallavotti-Nicolo trees; in Section 13.41 we 
discuss the flow of the running coupling constants; in Section 13.51 we derive the 
Ward Identities that allow to control the flows of the effective charge and of the 
effective photon mass, and finally in Section 13.51 we show that the infrared fixed 
point of the RG procedure is described by a Lorentz invariant effective theory. 

Remark 3 From now on we shall assume that the limits h* — )■ —00, K — t- -|-oo 
have been taken; this can be done safely since, as it is implicit from the discussion 
of AppendixWi and from the one that follows, the perturbative series that we shall 
get in presence of finite h* , K would converge uniformly order by order to a limit 
as h* —7- —00, K —7- -|-oo. 

3.2 Multiscale analysis 

In order to compute the partition function (|2.4[5p we will use standard functional 
Renormalization Group methods, [271 [Sj [61]. The integration of (I2.415P will be 
performed in an iterative way, moving from high to small momentum scales; the 
procedure will be similar to the one discussed for an effective continuum model of 
graphene in [38]. At the n-th step of the iteration the functional integral ()2.4[5p 
is rewritten as an integral involving only the momenta at a distance proportional 
to M~" from the singularities k = p"^ for the fermions and p = for the bosons. 
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and both the propagators and the interaction will be replaced by "effective" ones, 
which are renormalized by the integration of the momenta on higher scales. We will 
start by setting the notations of the zero-th scale, corresponding to the outcome 
of the ultraviolet integration; in particular, we rewrite the fermionic field as the 
sum of two quasi-particle fields, supported on disjoint sets centered around the two 
singularities pp. Then, we will discuss the single scale integration, and finally we 
will show that the outcome of the integration can be rewritten using a "relativistic" 
notation, that will allow us to show an emergent Lorentz invariance. 



The zero-th scale. As a starting point, we rewrite the propagators as 

^("•■)(k' + P-) =: g(<0)(k') , u;(*-)(p) =: w^^^\p) , dS^l) 

where 





Zo koTO + vok' ■ 

^(<o)(p) ^(l + R'ip)), 

^IPI 

with xo(k') := x(|k'|), ||i?o,^(k')|| < (const.)|k'|, [ fi^p)! < (const.)|p|, vq = |i, 
Zq = 1 and the matrices have been defined in (|2.4il8|) . Correspondingly, we 
define 

moreover, the following notations will be useful in the following: 

Then, we rewrite (|2.4i5p as, in the limit K — )■ +co, h* — t- — oo, 

where: (i) v(°H^^-°^ ^^-°^) is equal to V(^'(^-''-), ^(*-^-)) once that ^'(^■'■•) and 
have been rewritten according to ()3.2I3|) that is: 

/n m 

n,m^u i^.f; i=l i=l 
n+m>l 

(m 2n \ 

^p,+^(-irkn , doe) 
j=i i=i / 

where u = (ui, . . . ,uj2n), QL = (fi, • • • ,cr.«), /_f = (fii, . . . ,iJ.m), and the kernels are 
defined as, see (|2.4l7p . 

Wl]^,,,^jm, {Pj}) := W2n,n.,pA{K + Pf'}, {Pj}) ; 



38 



3. The infrared integration 



(ii) the integration measures are: 



" k'G© 



n n dAi^^uA^^'^ 



•exp{-(2MA)-i ^ i.(|°)[#(p)]-M3}, (13218) 



where A/o,!!!, A/q^a are normahzation factors. 

Single scale integration and renormalization. Setting XhO^') '■= x{^~^\^'\)^ 
we start from the following identity: 

X(|k'|)= A(k'), A(k'):=x^(k')-Xh-i(k'); m^) 

h=—oo 

let '^k',<x,p,a; = El=-oo^'£l,p,u, and 4,p = EL-ooAI%, where {^-C^)}, {i^^)} are 
independent free fields with the same support of the functions fh introduced above. 
We evaluate the functional integral (j2.4[5|) by integrating the fields in an iterative 
way starting from A^^\ We want to inductively prove that after the integration 

we can rewrite: 



where: (i) P<h{d'^^-^^) and P<hidA^-'^^) have propagators 

Xh(kO 1 



10) 



^(1 + ^'(P)), 



E2111) 



with ||i?/,,^(k')|| < C|k'|; (ii) V^'*) has the form 

/n m 
W^t K TT^MP- 

11 k^-_-^,(Ti,p2i-l,'<J2i-l k^i,o-i,P2i,1^2i 11 '^"P» 

-t,m^u iL.;: i=l 1=1 



n+m>l 
.(ft) 



2n 



\j=i i=i / 

and F/i, Z/i(k'), ?}/i(k') and the kernels W2n]m,p,fi,u) '^i^^ defined recursively. For- 
mulas (l3.:jllUl) - are trivially true for K~=0, see dMHI, (15:^ and ^M^. 
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In order to inductively prove (|3.2[lU|) . we split V^'') as £V^^^ +7^V(^^ where 
TZ = 1 — C and C, the localization operator, is a linear operator on functions of the 
form ()3.2ll2p . defined by its action on the kernels M^i^^^y.^^^^ in the following way: 

'CW^S!p,-(k') ■■= <o!p,a,(0) + fc;5«<o!p,c,(0) , (132113) 

^W^%iP) ■■= <4(0) +PadM'^JO) , £<4(pi,p2) := <4(0,0) > 

and CW^^ := otherwise. By the symmetries listed in Lemma 12.3. H which are 
preserved by the multiscale analysis, it turns out that, see Appendix lEl 

<i>) = o, <3>,o)=o, wiS,,jo) = o, 

W^'SA^) = -6,u,,M%,,, , a^Tl^i'^Uo) = (132114) 
and, moreover, that 



12115) 

with Zf,^h, K,h, real, and zi^h = Z2,h, >^i,h = h,h, i^i,h = i^2,h- We can renormal- 
ize P<h{d'il)^-^'') by adding to the exponent of its gaussian weight the local part of 
the quadratic terms in the fermionic fields; we get that 

j P<;,(d^(^'^))P<,,((i^(^^))e^''^(v%:*,A) ^ (J32J16) 



where th takes into account the different normalization of the two functional inte- 
grals, ■pC') is given by 

=: V^^\^,A) - C<s,V^''\^,A) , (132117) 
and P<h{d'4>^-^^) has propagator equal to 

Zh^i{\^')koTl + Vh-i{\^')k' -T^ 

with 

Z/,„i(k') = Zh{\^) + ZhZ^^hXhO^) , (132119) 
Z^_i(k')i)/^_i(k') = Zh{\^)vh{\^) + ZhZi^hXh{\^) 

and, setting (1 + i?/,_i,^(k'))^^ =: 1 + ^;^_l,a;(k'): 
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After this, defining Zh~i := Zh~i{0), we rescale the fermionic field so that 

V^^\V^^,A) = V^''\,/zZ~i'^,A) ; (132121) 

therefore, setting 

Zh Z^ 

Vh_i := Vh-iiO) , eo,h := — Ao,h , ei^hVh-i = e.2,hVh-i ■= ^\,h , 



122) 
we have that: 



32123) 



where 



.(<'^) ._ ^ f^^i<h)+>TY^O^{<h)- 
:= .»-,i:/=|^4ii^;,-ij:f„4i.ft- 6324) 

After this rescahng, we can rewrite (j3.2[16|) as 

P<;,(d^(^'^)d^(^'^))e^'''(^'v^*) = e^^l^l*^ J P<^„i(d^I/(^'^-i)d^(^'^-i)) • 

where have propagators given by (|3.2[lip (with h replaced by h—1) 

and ^(^),A^^^ have propagators given by 

9^'\^') ^ A(kO (1 + Rh-iA^')) - , ^ . M 

^.-1 fcorO+t;,_i(k')^'-f^ ' ^/.-i(k')^'^ ^ 

^(h)(p) ^ Mpl^i + /?'(p)) . ([32J26) 
^iPl 

At this point, we can integrate the scale h and, defining 

(122127) 

our inductive assumption (|3.211Up is reproduced at the scale h — 1 with F^-i := 
F'h + ^/i + -^/i- Notice that (|3.2i27p can be seen as a recursion relation for the 
effective potential, since from (|3.2117p . (|3.2i2ip it follows that 

V^'^^A, v^^) = v^^\a, ^/Z^i^) = v'^^Xa, ^h^p) - C^V^^\A, ^hi^) . 

(IS2128) 

The integration in (I3.2127P is performed by expanding in series the exponential in 
the r.h.s. (which involves interactions of any order in ^ and A, as apparent from 
(I3.2I12P ). and integrating term by term with respect to the gaussian integration 
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P{d'i>^''^)P{dA^^^). This procedure gives rise to an expansion for the effective po- 
tentials V^'*^ (and to an analogous expansion for the correlations) in terms of the 
renormalized parameters {e^^k,i^fi,k-i, Zk-i,Vk-i}h<k<i, where ei = e, which can 
be conveniently represented as a sum over Feynman graphs according to rules that 
will be explained below. We will call {e^^fc, f^,fc}h<fc<o effective couplings or running 
coupling constants while {e^^k}h<k<o are the effective charges 

Note that such renormalized expansion is significantly different from the power 
series expansion in the bare coupling e; while the latter is plagued by logarithmic 
divergences, the former is order by order finite. 

By comparing fl3:2T5]l with (j:i2ll()ll . (13.21121) and (13.21231) . we see that the inte- 
gration of the fields living on momentum scales > M'^ produces an effective theory 
very similar to the original one, modulo the presence of a new propagator, involving 
a renormalized velocity Vh and a renormalized wave function Zh, and the presence 
of a modified interaction V^'^^. 



Relativistic notations. We conclude this Section by showing that the effective 
theory on a given scale h can be rewritten as a "relativistic" theory, where the 
fermions propagate with a velocity Vh-i- Let us define the 4- component Grassmann 

fields ip^i^\ V'k^a'' ^ follows: 

j^i<h),T ._ (^{<h)- ^{<hy *{<h)- \ (13 2129) 

■- [^k',c7,l,+ ^k',<7,2,+ ^k',<7,2,- ' U^"^^) 

Wk',a — ['ifk',a,2,-^ ^k',<7,l,- -'^k',<7,l,+ -^k',a,2, + J ' (13^3Uj 

and the Euclidean gamma matrices 7^, /i = 0, 1,2 as: 

it is easy to see that {7^,71/} = —25^^y. With these notations, the fermionic inte- 
gration measure can be rewritten as: 

p<,{d^^^'^) = M,}, n #k',.,rfV'k',.,exp{-^/ ^^Lfj^^^'^HkO^Vg;^} 

(132132) 

where Mh,ii) is a normalization factor, and the fermion propag ator ^(^'')(k') IS given 
by 

g(.<h){y^) .- ^^lM 1 ^ (1 + i?^k')) , (132133) 

^ ' Z,(k')ifco7o + ^/.(k')fc'-7 
with ||i?/i(k')|| < (const.) |k'|; analogously, we can rewrite the fermionic current 
defined in (I3.2I24|) as: 



.{<h) s:^ f dk' —(<h) (<h) -r{<h) ^ f dk' -(<h) ^ A<h) 



™^ 134) 

In a fully relativistic theory, that is with bare propagators given by (|3.2133|) , (|3.2126|) 
with /i = and Zq = vq = 1, -Ro(p) = -^'(p) = 0) a-iid bare interaction given by 
(|3.2i24|) with /i = and e^^o = e, f^^o = ^1 Lorentz symmetry implies that Vh = 1 
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{the speed of light is not renormalized) e^^h = eo,/i, i^^^h = '^o,h] see [35]. On the 
contrary, the lack of Lorentz symmetry in our model, due to the presence of the 
lattice, has two main effects: (i) the Fermi velocity Vh has a non trivial flow; (ii) the 
marginal terms in the effective potential are defined in terms of two charges, namely 
eo,h and ei^h = ^2,h-, which are different, in general. 



3.3 Tree expansion 

The iterative integration procedure described above leads to a representation of 
the effective potentials in terms of a sum over connected Feynman diagrams, as 
explained in the following. The key formula, which we start from, is (13.21270 . which 
can be rewritten as 



n>l \ / 



with the truncated expectation on scale h, defined as 

4^(X(AW,^W);n) := ^ log / P,(dM/W)P,(dAW)e^^(^"^'''^"^^) |^^^ (12312) 

If X is graphically represented as a vertex with external lines corresponding to A^'^^ 
and the truncated expectation (j3.3[2|) can be represented as the sum over the 
Feynman diagrams obtained by contracting in all possible connected ways the lines 
exiting from n vertices of type X. Every contraction corresponds to a propagator 
on scale h, as defined in (|3.2126p . Since V^'') is related to V'^'*^ by a rescaling and 
a subtraction, see (I3.2[17|) and ()3.2[21|) . Eq. (l3.3[ip can be iterated until scale 0, 
and V(^-^) can be written as a sum over connected Feynman diagrams with lines 
on all possible scales between h and 0. The iteration of ()3.3[ip induces a natural 
hierarchical organization of the scale labels of every Feynman diagram, which will 
be conveniently represented in terms of tree diagrams. In fact, let us rewrite V^'^^ 
in the r^h.s. of (13^11]) as V^^\./Zi~[^ , A) = CV^^^^/Zj^'i' , A) + TZV^^^^/Zi^'^' , A), 
where C := C — Cq/, see (|3.2ll7p . Let us graphically represent V^'*^ CV^'^^ and 
as in the first line of Fig. 13. H and let us represent Eq. ()3.3[ip as in the second 
line of Fig. 13. 1[ in the second line, the node on scale h represents the action of 
S'^. Iterating the graphical equation in Fig. 13.11 up to scale 0, we end up with a 





Figure 3.1. Graphical interpretation of Eq. (|3.3lip . The graphical equations for £V'^'' ^\ 
are obtained from the equation in the second line by putting an £, TZ label, re- 
spectively, over the vertices on scale h. 
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representation of V^'^^ in terms of a sum over Gallavotti-Nicold trees r [27\ \5\ I34| : 



N>lTeTh,N 

where Th,N is the of rooted trees with root r on scale hr = h and N endpoints, see 
Fig. E21 The tree value V'^^\t) can be evaluated in terms of a sum over connected 



V*-i) 



= S 



h-l 



1 



Figure 3.2. The effective potential V'''"^-' can be represented as a sum over Gallavotti - 
Nicold trees. The black dots will be called vertices of the tree. All the vertices which are not 
endpoints except the first {i.e. the one on scale /i) have an TZ label attached, which means 
that they correspond to the action of T^S^^, while the first represents The endpoints 
on scales < 1 correspond to the first two graph elements in Fig l3.31 which are associated 
to the two terms in CV'-'^\ see p.2l23p : instead, the endpoints on scale 1 correspond to 



Feynman diagrams, defined by the following rules. 

With each endpoint u of r on scale < 1 we associate a graph element of type 
e or v, corresponding to the two terms in the r.h.s. of (I3.2l23p . see Fig. 13. 3| if v 
is on scale 1 we also have to consider the graphs arising in the integration of the 
ultraviolet degrees of freedom, contributing to the kernels of TZV^^^ , see Appendix 



e) V) 




Figure 3.3. The graph elements corresponding to the two terms in the r.h.s. of (|3.2l23p . 



We introduce a field label f to distinguish the fields associated to the graph 
elements (any field label can be either of type A or of type the set of field labels 
associated with the endpoint v will be called ly. Analogously, if v is not an endpoint, 
we call the set of field labels associated with the endpoints following the vertex 
V on T. 

We start by looking at the graph elements corresponding to endpoints on scale 
1: we group them in clusters, each cluster G„ being the set of endpoints attached 
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to the same vertex v on scale 0, to be graphically represented by a box enclosing 
its elements. For any on scale (associated to a vertex v on scale that is 
not an endpoint), we contract in pairs some of the fields in Ua.GG^-^uji in such a 
way that after the contraction the elements of are connected; each contraction 
produces a propagator g!;^^ or w^^\ depending on whether the two fields are of type 
^' or of type A. We denote by Zy the set of contracted fields inside the box Gy and 
by = Iy\ ly the set of external fields of G^; if f is not the vertex immediatly 
following the root we attach a label TZ over the box Gy, which means that the TZ 
operator, defined after (|3.2112p . acts on the value of the graph contained in Gy. 
Next, we group together the scale-0 clusters into scale- (-1) clusters, each scale- (-1) 
cluster Gy being a set on scale-0 clusters attached to the same vertex v on scale — 1, 
to be graphically represented by a box enclosing its elements, see Fig. 13.41 




-2-1 1 



Figure 3.4. A possible Feynman diagram contributing to and its cluster structure. 



Again, for each v on scale —1 that is not an endpoint, if we denote hy vi, . . . , Vg^ 
the vertices immediately following v on r, we contract some of the fields of U^l-^Py^ 
in pairs, in such a way that after the contraction the boxes associated to the scale- 
clusters contained in Gy are connected; each contraction produces a propagator 
g(~^) or w^~^\ We denote by X„ the set of fields in U^^^Py^ contracted at this second 
step and by Py = Ull-^Py^ \ Iy the set of fields external to Gy] if v is not the vertex 
immediatly following the root we attach a label TZ over the box Gy. 

Now, we iterate the construction, producing a sequence of boxes into boxes, hi- 
erarchically arranged with the same partial ordering as the tree r. Each box Gy is 
associated to many different Feynman (sub-)diagrams, constructed by contracting 
in pairs some of the lines external to Gy^, with Vi, i = 1, . . . ,Sy, the vertices imme- 
diately following V on r; the contractions are made in such a way that the clusters 
Gyj^, . . . ,Gy^^ are connected through propagators of scale hy. We denote by P^ 
and by Py the set of fields of type A and ^, respectively, external to Gy. The set 
of connected Feynman diagrams compatible with this hierarchical cluster structure 
will be denoted by T[t). Given these definitions, we can write: 

= E / n § n T^vaKc) . ea4) 



Val(a) = [ n ^1tAp,][ n V^^Lfir/5fj(/)M/)]^(-o)Vai(g) 



vai(g)=(-ir/ n (y^) ^ (n^^n( n 

J.. „„i „_ ^^h,.~2^ Syl ^ ^ ^ _ _ 



V not c.p. 



Idv V* e.p 
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where: ( — 1)'^ is the sign of the permutation necessary to bring the contracted 
fermionic fields next to each other; in the product over / G P*, e(/) = =b depending 

on the specific field label /; S{vo) = s(^J2feP^ P/ ~Z]/gp* {~^Y^'^^^'f^ '■, the integral 
in the third line runs over the independent loop momenta; Sy is the number of 
vertices immediately following v on t; TZ = 1 — C is the operator defined in (|3.2113|) 
and preceding lines); ay = if f = fo, and otherwise a„ = 1; g^J"^ is equal to g^'^'^ or 
to w^'^^ depending on the fermionic or bosonic nature of the line £, and i G v means 

(k) 

that £ is contained in the box Gy but not in any other smaller box; finally, K^, is 
the matrix associated to the endpoints v* on scale k + 1: setting 

eo,k ■= eo,fc , ej^k ■= Vk-iej,k , ([3315) 

it is given by e^,A:r^ if v* is of type e, by —M^v^j^^k if v* is of type v, or, if v* is 
on scale 1 and is not of type e, v, it is equal to one of the kernels contributing to 
7^y(°), see Appendix Id! In (j3.314|) it is understood that the operators TZ act in the 
order induced by the tree ordering (i.e., starting from the endpoints and moving 
toward the root); moreover, the matrix structure of g^l^^ is neglected, for simplicity 
of notations. 

3.3.1 An example of Feynman graph 

To be concrete, let us apply the rules described above in the evaluation of a simple 
Feynman graph Q arising in the tree expansion of V^*-^). Let g be the dia gram m 
Fig. 13.51 associated to the tree r drawn in the left part of the figure; let us assume 
that the sets Py of the external lines associated to the vertices of r are all assigned. 




h-l h h+l h+2 



Figure 3.5. A possible Feynman diagram contributing to V^'^ and its cluster structure. 



We can write: 



Val(g) 



'4!2! Zh^2 '"'''"^^ ^/^i.'^^k.a.o; 



3316) 



•7^ 



(2^ 



r^i^5if+')(k + p + q)r^i^u;('^+i)(q)l5('^)(k + p)r^i [^gi" , 



where 7^[F(k + p)] = F(k + p) - F(0) - (k + p) • VF(0) = + p^){K + 

Pu)df^duF(k*). Notice that the same Feynman graph appears in the evaluation of 
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other trees, which are topologically equivalent to the one represented in the left 
part of Fig. I3.5l and that can be obtained from it by: (i) relabeling the fields in P^^, 
Pyg, (ii) relabeling the endpoints of the tree, (iii) exchanging the relative positions 
of the topologically different subtrees with root vq. If one sums over all these trees, 
the resulting value one obtains is the one in Eq. (|3.316|) times a combinatorial factor 
2^ • 3 • 4 (2^ is the number of ways for choosing the fields in P^^ and in P^^,; 3 is 
the number of ways in which one can associate the label v to one of the endpoints 
on scale h + 1; 4 is the number of distinct unlabelled trees that can be obtained by 
exchanging the positions of the subtrees with root vq). 

3.3.2 Dimensional bounds 

We are now ready to derive a general bound for the Feynman graphs produced by 
the multiscale integration; before stating our result, we define the order of a tree 
as follows. Let Tk-h.n be the set of trees with endpoints obtained expanding 
into trees the effective potentials TZV'^^^ attached to the endpoints of the elements of 
U^=i Th^n- Let us define the order of an endpoint v in the following way: the order 
of the endpoints of type e,v is 1; the order of an endpoint on scale + 1 is equal 
to the number of wavy lines exiting from the corresponding graph element. We 
define the order 0(r) of a tree r € TK;h,N as the sum of the orders of its endpoints. 
In general, trees belonging to TK;h,N may be of arbitrary order > N ^ since the 
endpoints on scale K + 1 correspond to vertices with an arbitrary number of wavy 
lines (because of the expansion of the exponential of the photon field performed 
before the ultraviolet integration, see Eq. (|2.2il5p ). 

Now, let W^2^!m!^p,a;,^ ^e the contribution due to the trees of order N to the kernel 
m puj fi which in the following we shall simply refer as "the contribution of 
order N"), that is 

win!^,p..,;.({ka,{p.}) = E E E E vai(g) 

0{t)=N 1^41=™' 
l<l=2n 

oo 

= ■■ EW^in?i,^,M(ika,{P.}), »7) 
N=l 

where the * on the sum indicates the constraints that: UjgpAjp/} = U"i;^{pj}; 

U/eP* {k/} = Uf^jk'j; U^eP* W/)} = p; U^^p* {a;(/)} = cv; U^^pA M/)} = !£■ 
The A^-th order contribution to the kernel of the effective potential admits the 
following bound. 

Theorem 3.3.1 (A^! bound) Let Sf,. = max/i<fc<i{|e^^fc|, be small enough. 

If < e^'^h and < Vk_i < C, for all h < k < and a suitable constant 

C > 0, then 

< (const.)^s;:^(f )! MM3— 2n) , 

where \\W^n,m,p,uiJ\ ■= ^^V{k'J,{p,}\Wi'^]m,p,c^A{'^'i}APj})\- 
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The factor 3 — 2n — m in (|3.318p is referred to as the scaling dimension of the 
kernel with 2n external fermionic fields and m external bosonic fields; according 
to the usual RG therminology, kernels with positive, vanishing or negative scaling 
dimensions are called relevant, marginal or irrelevant operators, respectively. Notice 
that, if we tried to expand the effective potential in terms of the bare coupling e, the 
A^-th order contributions in this "naive" perturbation series could not be bounded 
uniformly in the scale h as in (|3.3[8p . but rather by the r.h.s. of (|3.3[8p times \h\^ , 
an estimate which blows up order by order as /i — )■ — oo. 

Proof. To begin, let us assume that all the endpoints are on scales < 1; the general 
case can be worked out in a completely analogous way, and it will be discussed later. 
Therefore, for the moment we shall consider the contribution of trees r G Th,N- 
Plugging the bounds 

||ff(f )(k')|| < const • M-'' , J dk'||g(,^)(k')|| < const • M^^ , 
||i2;('')(p)|| < const • , J dk\\w^''\p)\\ < const • M^'^ , (I33l9) 

into (I3.3l4p . we find that, if r G Th,N and Q G r(r). 



|Val(C?)| < (const.) 



n 

V not o.p. 



C=2|p*| 



10) 



V not e.p. 

v>vo 



where: is the number of propagators £ £ v, i.e., of propagators £ contained in 
the box Gy but not in any smaller cluster; Sy is the number of vertices immediately 
following V on r; is the number of end-points of type v immediately following v on 
T (i.e., contained in Gy but not in any smaller cluster); v' is the vertex immediately 
preceding u on r and Zy = 2 if |P*| = \Py\ = 2 or \P^\ = \Pv\ = 2, Zy = 1 
if |P*| = 2\P^\ = 2 or IP^^^I = and \P^\ = 3, and Zy = otherwise. The 
last product in (|3.3[lUp is due to the action of TZ on the vertices v > vq that 
are not end-points. In fact, the operator TZ, when acting on a kernel W^\^\p,k') 
associated to a vertex v with \Py\ = 2, \P^\ = 
of first order in its Taylor expansion around p 



1, extracts from 1^2'}"^ the rest 
= k' = 0: if |^';^)(p,k')| < C, 
then \nW^'l^\p,k')\ = ^\{pdp + k'ff^)W^'^{\p*,k*)\ < {const. ) M-''-+^-' G, where 
M~'^" is a bound for the derivative with respect to momenta on scale hy and M^^' 
is a bound for the external momenta p, k'; i.e., TZ is dimensionally equivalent to 
]\^-{K-h^i) ^ The same is true if TZ acts on kernels Wq']^\p,(\). Similarly, if TZ acts 
on a terms with |P„| = 2, it extracts the rest of second order in the Taylor expansion 
around k' = 0, and it is dimensionally equivalent to k'^9^, ~ M~'^^^''~^^'\ As a 
result, we get (I3.3ll0p . 

Now, let n% (n^) be the number of vertices of type e (of type v) following v on 
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r. If we plug in (|3.3[10|) the identities 



{K - h){s^ - 1) = (^^ - hv'){nl + nl - 1) 

V not c.p. V not c.p. 



IP,, 



V not c.p. V not e.p. 

Y {hv-h)m'',= Y {hv-h,')K OH) 

V not e.p. V not e.p. 

we get the bound 

C^2|p*| 



|Vai(g)| < (const.)^e;y^M^(3-|P„J) -Q V)(3-|n|-..) . 

^ Si) . 



,1 J-J- ,1 

"■"O- V not e.p. 

(133112) 

In the latter equation, 3 — |P^| is the scaling dimension of the cluster G^,, and 
3 — |P^| — Zy is its renormalized scaling dimension. Notice that the renormalization 
operator TZ has been introduced precisely to guarantee that 3 — — z,, < for 
all f, by construction. This fact allows us to sum over the scale labels h < < 1, 
and to conclude that the perturbative expansion is well defined at any order of 
the renormalized expansion. More precisely, the fact that the renormalized scaling 
dimensions are all negative implies, via a standard argument (see, e.g., [51IM])) the 
following bound, valid for a suitable constant C: 



iP^f)^^. I 



< (const. )^e^^M''(3-m-2n) . ([^ 13) 

£e2|P*| 

. ^ ^ Y[ " " i^^(h^-h,,)(3-\PJ-z^) ^ 

TeTh,N Ger(T) V not e.p. 

l<l=2n 



from which, after counting the number of Feynman graphs contributing to the sum 
in ^MM, (IH^ follows. 

An immediate corollary of the proof leading to (|3.3[8p is that contributions from 
trees in Th,N with a vertex v on scale hv = k > h admit an improved bound with 
respect to ([3:318]) . of the form < (const. )^ej^(Ar/2)! M'^(3-|p.o ^j, ^^^^y 

< < 1; the factor M^^^'^'^ can be thought of as a dimensional gain with respect 
to the "basic" dimensional bound in (j3.3[8p . This improved bound is usually referred 
to as the short memory property (i.e., long trees are exponentially suppressed); it 
is due to the fact that the renormalized scaling dimensions = 2> — \P^\ — z.^ m. 
(|3.3ll2p are all negative, and can be obtained by taking a fraction of the factors 
j^{hv-h^,)Dv associated to the branches of the tree r on the path connecting the 
vertex on scale k to the one on scale h. 

We conclude by discussing the contribution of the trees r G TK;h,N with end- 
points on scales > 1. In this case, the bound ()3. 31120 is replaced by the analogous 
bound obtained by taking into account the contribution of vertices on scales > by 
matching the bound ()3.3ll2p with the ultraviolet bound ()D[16P : the final result is 
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that 

|Vai(g)| < (const. )^e^ — ^m'*(3-|p„oI) tt ^±J_^m^K-K,)d, _ ((33114) 

^0' V not e.p. 

v>vo 

where the renormahzed scahng dimensions Dy are always negative by construction, 
a„ = 1 if hy < or otherwise. ■ 



Symbol 


Description 


r 
r 

Vo 
hy 

Th,N 
Iv 

Gy 

J 
J^y 

Sy 

p# 
^v 

r(r) 

rriy 

v' 

nf 

Dy 
tty 


Gallavotti - Nicolo (GN) tree. 
Root label of the tree. 

First vertex of the tree, immediately following the root. 
Scale label of the tree vertex v. 

Set of GN trees with root on scale hr = h and with N endpoints. 
Set of field labels associated with the endpoint of the tree v. 
Cluster associated with the tree vertex v. 

Set of contracted fields inside the box corresponding to the cluster 

Gy. 

Set of external fields of Gy . 
i-th vertex immediately following v on the tree. 
Number of vertices immediately following the vertex v on the tree, 
oex oi neiQS oi xype y/= — w exteinai xo \jy. 
Set of connected Feynman diagrams compatible with the hierar- 
chical cluster structure of the tree r. 

Number of propagators contained in Gy but not in any smaller 
cluster. 

Number of end-points of type v immediately following v on the 
tree. 

Vertex immediately preceding v on the tree. 
Number of vertices of type = e^v following v on the tree. 
Improvement on the scaling dimension due to the renormalization. 
Renormahzed scaling dimension (equal to the scaling dimension if 
hy > 0). 

Integer equal to 0, 1 depending on whether hy <0,hy> 0. 



Table 3.1. List of symbols introduced in Section [3?3l 



3.4 The flow of the running couphng constants 

As a consequence of the interative integration scheme the effective parameters 
Zh, Vh} obey to a nontrivial recursion relation; in fact, the difference of 
renormahzed parameters on scales h and h + 1 can be expressed in terms of a well 
defined function of all the renormahzed parameters on scales > h + 1. This function 
is called the Beta function, and it can be written as a series in {e^^k, f^,fc}fc>h+i with 
coefficients depending on ffc}fc>/i+i, admitting A^! bounds in the sense 

of Theorem 13.3.11 A crucial point for the consistency of our approach is that the 
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running coupling constants e^^/^, Ufj^^h are small for all h < 0, that the ratios Zh/Zh^i 
are close to 1, and the effective Fermi velocity Vh does not approach zero; even if we 
do not prove the convergence of the series but only A^! bounds, we expect that our 
series gives meaningful information only as long as the running coupling constants 
satisfy these conditions. In this Section we describe how to control their flow; we 
shall proceed by induction: we will first assume that Eh = niaxk>h+i{\efi,k\, Wfi,k\} 
is small, that Z^/Zk-i < e'^'^fc and C^^ < < C for all /c > /i + 1 and a suitable 
constant C > 0, and we will study the flow of Zf^^ and under these assuptions. 
We will show that, asymptotically as /i — >■ — oo, 

Zh ~ M-^'^ , Vh ~ v.oo - {v-oo - v)M'i^ , daai) 



„2 n^2 



+ 0(e^), fi='^ + 0{e^), v.^ = l + 0(l 



where 1 is the speed of light in our units. These behaviors are true at all orders, in 
the sense that they are obtained by truncating the Beta functions of Z^ and to 
any finite order in the running coupling constants and taking the bare charge small 
enough uniformly in the momentum scale h. Finally, as discussed in Section [3.51 to 
show the boundedness of the flow of the running coupling constants we shall use 
nonperturbative informations coming from lattice Ward identities; in particular, we 
will get that eo,-oo = ^i-oo, which implies that V-oo = 1. 

3.4.1 The Beta function 

The multiscale integration, and in particular formulas (I3.2ll4p . ()3. 21150 . ()3.2ll9p . 
(|3.2[22p . imply the following flow equations: 



h 



1 + z^^h =:1 + I3h, (13312) 



Vh-i = -^{vh + z^^h)=:Vh + PK, 

v^^H = -M-^W^'^l^^^{0)=:Mv^^h+i+P''^,h+i, (13314) 
eo,h = ^ Ao,h =: eo,h+i + /3o,h+i > ([2315) 

ei,h = /^^=:ei,,+i + /3iVi, (11316) 

and 62, h = eijj. The Beta functions appearing in the r.h.s. of flow equations are 
related, see (|3.2ll4p . (|3.2I15|) . to the kernels W^in^rn^p /i' that they are expressed 
by series in the running coupling constants admitting the bound (|3.3I8|) . For the 
explicit expressions of the one-loop contributions to the Beta function, see below. 



3.4.2 The flow of Fermi velocity 

In this section we show that, under proper assumptions on the flow of the effective 
charges, the effective Fermi velocity tend to a limit value V-oo, which will be 
explicitly computed in Section 13. 6t remarkably, we will show that, no matter how 
small the bare velocity v is, the effective Fermi velocity f_oo is equal to 1, which is 
the speed of light in the units we have chosen. 



3.4 The flow of the running coupling constants 



51 



Let us make the following assumption on effective charges and on the effective 
photon masses: 



with < ^ < 1 and e^^_oo = e + O(e^); we shall refer to the first line of (|3.4I7|) 
by saying that the charge tens to a line of fixed points. The remarkable properties 
(|3.4i7|) will be proven order by order in perturbation theory by using Wis, see 
Section 13.51 We start by studying the flow of the Fermi velocity; at lowest order 



(see Appendix IF.ip . its beta function reads, setting ^/j '■= yv^ ~ 



Ph 



logM 



eo,h^/i ^ arctan in /„ 2 ^o./i^h S - arctan 



Note that if eo,fc := &i,k^ then the r.h.s. of (|3.4i8p is strictly positive for all > 
and it vanishes quadratically in at = 0. The higher order contributions to 
have similar properties. This can be proved as follows: we observe that the Beta 
function is a function of the renormalized couplings and of the Fermi velocities 
on scales > h, i.e.: 

/^h = /3h({(eo,fc,ei,fc,e2,fc), (i^o,fc,i^i,fc,'^2,fc)>^fc}fc>J • (13319) 

We can rewrite fij^ as fij^ o+Z^h >0' where takes into account all the contributions 
coming from trees with at least one end-point on scale > 0; by the short memory 
property (see discussion after Eq. (|3.3[13p ) , the bound on admit an improve- 

ment of a factor M^^ with respect to the basic one. Moreover, we rewrite q as 

K^Q^'^^'h 0' "^here /3^'o'^' is given by q after replacing all the propagators with their 
relativistic approximations, namely those obtained setting Rh-i^^{k!) = R'{p) = 
in (|3.2l26p : therefore, q admits an improvement of a factor M^^ in its dimensional 

1 J rpi nv,rel nV,reLO , nV,reLl , nV,reL2 , nV,reL3 , av 

bound. Then, we can rewrite Pf^'^ as Pj^'^Q + Pj^Q + p^Q + Pf^'^ + p^ Q, 
where = ©(e^M^^) and ' ' 

^«,re/,o ^ /3j^;o'''({(eo,fc,eo,fc,eo,fc),(i^o,fc,t'o,fc,^'o,fc),l}fc>J , ([33110) 
^^,re/,i ^ /3j^;o^'({(eo,fc,eo,fc,eo,fc),(^'o,fc,'^o,fc,i'o,fc),''^fc}fc>/j) - 

^^,reZ,2 ^ /3l'^Q''\{{eo^k,eo,k,eo,k),{^0,k,^l,k,l'2,k),Vk},,^fJ - 
~f^hfi'^ i{i'^0,k, eo,fc, eo,A:)) {l^O,k, ^0,k, l^O,k) , Vk} , 

Ph}}^'^ = /3h'o^'({(eo,fc5ei^fc,e2,fc),(fo,fc,'^i,fc,t'2,fc),^^fc};,>^) - 

-Ph'fi\{^'^0,k,eo,k,eQ^k), iyo,k, l^l,k, ^2,k)-,Vk]j^yy) ■ 
By the short memory property we get: 

/3h 0'''' = O(eo,/.(eo,h - ei,,)) • (133111) 
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Using (|3.4i4p and a telescopic argument similar to the one leading to Eq. (|3.4llip . 
we also find that i/o,/i ~ '^i,h can be written as a sum of contributions of order 
eo,h(eo,/i - ei,h) and of order el f^{l - Vh). 

Now, let us consider /S^'q^^'^; this is the Beta function of the continuum relativis- 
tic theory studied in [38], where both photons and fermions live in the continuum 
and are regularized by an ultraviolet cutoff function XO; but where the effective 
charges and the photon masses have been replaced by those of our lattice model. 
However, we can rewrite eo,fc = eo,fc + eo,fc - eo,fc, t'o.fc = i^o.fc + I'o.fc - i^o.fe, where 
6o,fc) ^o,fc are the running coupling constants of the theory of [35] with bare Fermi 
velocity v = 1 and bare parameters e, i> chosen so that eo,/i = eo,/i, vq^k = z^o./i! 
it is straightforward to see that eo,fc — co^fc, z^o.fc ~ ^o,k can be written as sums of 
contributions 0{e^ M^^), 0(e^(l — ffc)), 0{e{ei^k+i — eo,fc+i))- Moreover, by rela- 
tivistic invariance we know that, see [38], if we replace at the argument of p,'"''^'^^''^ 



h,0 

the couplings eo,h) '^o,h with eo,/i, t'o.fc we get 0; therefore, at the end we find that, 
using the short memory property: 

PZo'^' = 0{el,M''^) + 0{eHl - v,)) + ©(e^ ,(eo,,+i - ei,,+i)) . ([23112) 

Therefore, we can write: 



Vh-i ^ _^ log M 



47r2 



8 4 

-e2(l - vh){l + A'^) + -e(l + S;,)(eo,h - ei,,,) 
5 3 



([33113) 



where: the numerical coefficients are obtained from the explicit lowest order com- 
putation ()3.4[8p : A'^ is a sum of contributions that are finite at all orders in the 
effective couplings, which are either of order two or more in the effective charges, 
or vanishing at Vk = 1; similarly, i?^ is a sum of contributions that are finite at 
all orders in the effective couplings, which are of order two or more in the effective 
charges. From (|3.4113p it is apparent that Vh tends as /i — )■ — oo to a limit value 

^-oo = 1 + ^(eo,-oo - ei,_oo)(l + C-oo) (IM114) 
be 

with C^oo a sum of contributions that are finite at all orders in the effective cou- 
plings, which are of order two or more in the effective charges. The fixed point 
(|3.4ll4p is found simply by requiring that in the limit /i — ?■ — oo the argument of the 
square brakets in (13.41130 vanishes. 

Using Eq. ()3.4[7p . we find that the expression in square brackets in the r.h.s. of 
(|3.4ll3p can be rewritten as {8e^/b){v-oo - Vh + -R'^)(l + A^), where: (i) A^ is a 
sum of contributions that are finite at all orders in the effective couplings, which are 
either of order two or more in the effective charges, or vanishing at = V-oo', (ii) 
R'f^ is a sum of contributions that are finite at all orders in the effective couplings, 
which are of order two or more in the effective charges and are bounded at all orders 
by M^^, for some < 9 < 1. Therefore, ()3.4ll3p can be rewritten as 

-Vh-i = {v-oo - Vh)(l - '"""^ logM^(l + ad) , ([23115) 

from which, using the fact that i?^ = 0{e^M^^), we get that there exist two positive 
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constants Ci , C2 such that 0: 



Vh 



3116) 



with 



fj = — log 



M 



l-t;_oologM^(l + ^': 



17) 



at lowest order, Eq. (l3.4ll7p gives 7?^^^ = 2e^/(57r^). 



3.4.3 The flow of the wave function renormahzation 

In contrast to what happens in the case of the other running coupling constants, 
the flow of Zh diverges as /i — )■ — 00 with an interaction-dependent power law. In 
particular, as it will be clear from the analysis of Section 14.21 the divergence of 
the wave function renormahzation, that is the vanishing of the quasi-particle weight 
ZJ^^, implies the anomalous scaling of the Schwinger functions. 

As for the Fermi velocity, assume that (|3.4I7|) holds; from (|3.4i2p it follows that, 
for two suitable positive constants Ci, C2: 

CiM^'^ <Zh< 02^^^ , r? = hm log^ (l + PI) ■ dSaHS) 
The lowest order contribution to is computed in Appendix W\ and it is given by: 



,,(2) _ logM 2 
Ph - 4^2 V^^l. 



\hVh 



ih - arctan ^h 



3119) 



therefore, from (|3.4118p . (|3.4119p we get that the lowest order contribution to r] is 
positive, and it is given by r]^'^^ = Now, let us briefly comment about the 

relation between Z/i, Vh and the functions Z(k') and i'(k') appearing in the main 
result, see drm]). If |k'| = M'', we define Z{k') := Zh and v{k') := Vh] for general 
|k'| < 1, we let Z(]i') and f(k') be smooth interpolations of these sequences. Of 
course, we can choose these interpolations in such a way that, if M'* < |k'| < M'^^^, 



Z{k') 
Zh 


- 1 


< 


Zh+i 
Zh 


- 1 


= 0(r/ log M) , 


v{k') 

V-00 - 


- Vh 

- Vh 


< 


Vh+l 

V-00 


- Vh 

- Vh 


= 0(r?logM) 



3120) 



therefore, we can replace in the leading part of the 2-point Schwinger function, 
which will be computed in Section 14.21 the wave function renormahzation Zj and 
the effective Fermi velocity Vj by Z(k') and v(k'), provided that the correction term 
R{k') in (ll.4[ip is defined so to take into account higher order corrections satisfying 
the bounds (I3.4[20p . 



^Ea. (|3.4ll6|l must be understood as an order by order inequality: if we truncate the theory at 
order A'' in the bare couphng e, both sides of the inequahty in Ea. (|3.4ll6[l are verified asymptotically 
as e ^ 0, for all iV > 1. 
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3.5 Ward identities 

To conclude the discussion on the flow of the running coupHng constants, we need 
to show that our assumptions of the flow of the effective coupHngs {e^^/i, f^^h} are 
indeed vahd. To prove this, we shah exploit the gauge invariance of the model to 
derive suitable lattice Ward identities; these are nonperturbative identities between 
Schwinger functions, which will imply the cancellations that we need in order to 
control the flows of {e^^/i, f^^/i}. Our strategy follows closely the one introduced 
by Benfatto and Mastropietro in the analysis on one dimensional Luttinger liquids, 

mm- 

Strategy. Before getting into the technical details, let us first give an idea of the 
strategy. At each RG step we shall introduce a model, to be called the reference 
model in the following; this model is not Hamiltonian and it is defined only in terms 
of a functional integral. The main features of the reference model are that: (i) it is 
defined in terms of a finite number of Grassmann variables, living on a space-time 
lattice: the spatial part of the lattice is given by the physical honeycomb lattice, 
while the temporal one is introduced by hand, and its mesh will be sent to zero; (ii) 
the bosonic sector of the model is regularized by an infrared cut-off on the same 
momentum scale h that we have to integrate in the full model. The key points of 
the strategy are the following ones. 

(i) The reference model can be studied using multiscale analysis and renormaliza- 
tion group, and the running coupling constants on scales greater than h are 
equal to those of the full model. 

(ii) After the scale h, that is after having integrated out all the photon fields, the 
reference model is reduced to an effective fermionic model which turns out 
to be superrenormalizable, and can be studied along the lines of [36]; this in 
particular implies that the running coupling constants of the reference model 
"cease to fiow" on scale h. 

(iii) The bosonic infrared cut-off does not break gauge invariance; that is, it is 
possible to derive Ward identities for this model which are equal to those 
which can be formally written in the full model, except for the fact that the 
Schwinger functions involved in the identities are all evaluated in presence of a 
bosonic infrared cut-off on scale h; we stress that the presence of the cut-off on 
the photons does not produce corrections to the Ward identities, in contrast 
to what would happen with a cut-off on the fermions (see |38j ) . 

(iv) Finally, the Ward identities provide relations for the bosonic two point func- 
tion, and relating the fermionic two and three point functions; from this iden- 
tities we get that the mass of the photon field on a given scale h is bounded 
by a quantity of order M^e^ , while the effective charge on scale h is close to 
the bare charge up to corrections bounded by O(e^) uniformly in h. 

These identities hide remarkable cancellations between Feynman graphs, and it is 
interesting to check them at lowest order; see Appendix IF. 5| where it is shown that 
at one loop in naive perturbation theory the renormalized mass of the photon field 
is zero and the effective charge is equal to the bare one. 
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3.5.1 The reference model 

Here we shall define the reference model; In particular, we shall prove that the 
generating functional of the Schwinger functions is left invariant by local phase 
transformations on the fermions (the Jacobian of the transformation is equal to 
1). This is a well known property of Grassmannian functional integrals, see [HT] 
for instance. But at the same time we know that, because of the invariance of 
the model under local gauge transformations, which are nothing more than a local 
phase transformation on the fermions combined with a shift of the photon field, the 
argument of the functional integral transforms in a suitable way; in this way we 
shall get an identity between generating functionals that gives rise to an infinite set 
of Ward identities. 



Definition and properties. Consider the sets 

Bn 

A := Au{xo = ^,n = 0,l,...,A^-l} , 
T>^,L := VLu{ko = y{n+^),n = 0,l,...,N-l}, 
Vp^L := Plu{po = — , n = 0,l,...,iV-l} ; ^l) 

With each k G ^^,Li P £ "^p.L we associate a fermion or boson gaussian field, 
respectively; with a slight abuse of notation, we shall denote these two fields with 
the same symbols used for the original model. The Fourier transform of the fermion 
field lives in the space-time lattice A, and it is given by 

The fermion and boson propagator are defined as follows. 



9e{^) ■= Xx(^o) 



-iko -tn*{k) 
-tn{k) -iko 



.(>h)f ^ _ f dp3 XKiPo)X[h,o]{\p\) 
{27T) p^+pI 



that is, they are equal to the "usual" ones except for the replacement of xk, x{\p\) 
with Xx, X[/i,o](|p|)5 respectively. The function Xk is a C°° function defined so that: 
(i) lime^oXKiko) = XK{ko); (ii) for |A;o| < XkO^o) = 1 and for |A:o| > M 
xUko) <ee-l'=ol^"''. 

A very important property of the reference model is that the Jacobian of local 
phase transformations on the fermion fields, namely 

x+{p-l){0,<5i),<7,p x+{p-l){0,<5i),a,p ' ^ ' ^1 ^ 
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is equal to 1. To see this, let 

:= n n dK<r,pdK<r,p ' «5) 

|A_| 

where the sequences {xj}, {x^} must be formed by distinct space-time points be- 
longing to A (otherwise Qa-A+ = 0, since the square of a Grassmann variable is 
zero, by definition), and consider the integral / QA_A+i^~ ,^^)', our claim is 
that 

J Qa_a+ (^- , ^+) = j Qa^a+ (e-^^"^- , e^^"^+) . ([3316) 
To prove (|3.5I6|) . we rewrite Qa-A+ in momentum space; it follows that: 

jv^QA_A,{^-,^-) = • (13317) 

Notice that the integral in the r.h.s. of (|3.5i7p is vanishing unless: (i) the sequences 
{kj} and {k^} are respectively formed by distinct elements; (ii) |^_| = |^+| = 
|P/3^i|. But the space-time points {xj}, {x^} are as many as the momenta {kj}, 
{k^}, and by construction I^'/j^l] = |A|; therefore, since by assumption the se- 
quences {xj}, {x^} are formed by distinct space-time points, '^^) must 
contain two Grassmann variables for each cr, p and space-time lattice site x, namely 
, p , and , w ? ^ 1 otherwise the outcome of the integration is 

x+(p-l)(0,(5i),CT,p x+(p-l)(0,(5i),(T,p' ^ 

zero. This proves ()3.516p . 

Consequences of gauge invariance. Let us define the generating functional of 
the Schwinger functions of the reference model as: 

:= J P(dvI/)P>,(A)e^(^+-')+^(*"^) , (13318) 

where the integration measures are determined by the propagators introduced in 
(|3.5l3p , and the fermionic source term is given by 

B{^, </>):= ^ / + <.,p^k,.,p > (13319) 

where St: ^ are Grassmann variables, with Fourier transform (ly^ ^ , := 

^^'''^P ^ ' ^x+{p-l){0,<5i),a,p 

(/3|A|)~^ SkGC/j L ^^''^^'^k o-p' ^''^oi'i cumbersome notations, we will not explicitly 
write the dependence on e, N , since at the end we will be interested in the limits 
e — 7- 0, — 7- -|-oo. For convenience, we rewrite the fermionic integration measure as 
follows: 

P(d^) := -^D^ exp I - (^+, g^-) I , (123110) 
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with 

where B(k) has been defined in I2.2i4l Using that the local phase transformation 
()3.5[4p has Jacobian equal to one, we can freely perform the replacement (I3.5l2p in 
the Gaussian weight of the Grassmannian measure, in the interaction and in the 
source term; neglecting formally the presence of the ultraviolet cutoff we would get 
the relation 

>vg)+°°^ (J, (/.) = (J + da, (/.e''^") , (!33112) 

which generates an infinite set of Ward identities, by deriving the l.h.s. and the r.h.s. 
once with respect to Op (which sets to zero the l.h.s.) and an arbitrary number of 
times with respect to the external fields J^,p, 0k,(7,p; notice that because of the fact 
that J^^x and ax always appear in the combination J^^x + dfiOiyi, as guaranteed 
by gauge invariance, the derivative w.r.t. cip is equal to —\p^ times the derivative 
w.r.t. j/i,p. However, to avoid formal expressions it is necessary to take into account 
the presence of the ultraviolet cutoff; the result is that for K fixed the analogous 
of (|3.5ll2p involves a i^-dependent correction term. The fact that cutoffs on the 
fermions produce corrections to the Ward identities is well known, [HI [12]; if they 
are not vanishing in the limit of cutoff removal we say that the Wis have anomalies, 
and this is indeed the case in one dimension, [12]. Here, on the contrary, as we 
prove in Appendix [G] the corrections are exponentially vanishing as — t- +oo. For 
simplicity, let us consider a function ax of the form ax := {f3A{s)~^ Z^pe-p^ ^ e~'P''ap; 
it is straightforward to see that, after ()3.5[2p . in the limit — )■ +oo: 

y^ff\j,<t>) ^ I p(dM,)p>^(dA)e^(^+^"+-^'*)+^(*"^^'^") • (133113) 

-(<I'+e'<=",ge-'<=a<I'-) + (l'+,g*-) + ('I'+c''=°,gc-''=°'I'-)-(iI'+,g<I'-)-ic(n,aoa) 

where 

a,p •' ^ 

Ma) = E/f (^^'^,p,.,p'I'k,.,p(-iPo)ape-(-^)-. 

The first two terms in the second line of ()3.5ll3p come from the free Grassmannian 
measure; the remaining three are produced by the interaction. Being interested in 
the derivative w.r.t. ap of (13.51 13p . we can rewrite the exponent appearing in the 
second line of (|3.5ll4p as follows: 

ie[(*+,ga^^) - (^+a,g^-) - (^+,^0^^) + (^+a,g^-) - (n,5oa)] +0{a'^) 

(133115) 

and it is easy to see that (repeated indeces are summed) 

^Va^-)=/f^«P^'.\p,,,- (133116) 



(^+,ga^' 
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and that 

[^+,ga^-) - f^+a, g^-) + (n, d^a) = ^ 17) 



Therefore, deriving with respect to Up the l.h.s. and r.h.s. of (|3.5[13|) and setting 
cip = we get that, in the Hmit e — )■ 0: 

d , . Jh.K] / ^ ^ , . [h,K], 



= (/3^A)T^W^f J(J + aa,<Aen + A^\"J(p; J, 0) , ([33118) 



tp 

with 

Aff{p;J,cb) := I §(^'^+p,.,,C,^Kk,p)^'i;,^,^,)^^^ , ([33119) 

-ie-^C^,(k,p) := ([23120) 
e-iPi(p'-i)(^-i(fco +po) - 1) [B{k + p)]^^, - e"'^^(^-i)(x^i(fco) - 1) [S(k)]^^, 

where the subscript J, i;^) in (|3.5119|) means that the average is taken in presence of 
the external fields. Formula (|3.5ll8p will be starting point in order to derive the 
Ward identities that we need to control the flows of {e^^/j, i^f_i,h}- 

3.5.2 Ward identity for the photon mass 

In this Section we shall derive the Ward identity for the photon field, which will 
imply that the photon field remains massless; the starting point will be formula 
(|3.5ll9p . Let us define: 



?KUp) := (W . <f'(./,0)l , (133121) 



P,K,27,M ■■= , n' ([23122) 

([23123) 

where the first line is the definition of the bosonic two point Schwinger function. 
Using that the action of da^ on yV^^'l'^\j + (9a, 0) is equal to the one of —ipfidj^^, 
from (I3.5[19|) it follows that: 

Apparently, this relation tells us that the regularized bosonic two point function 
is not vanishing; however, as it is proven in Appendix [Gl the correction term in 
(I3.5I24P vanishes exponentially in the ultraviolet limit K — )■ +oo. Therefore, from 
([53TM1) we get: 

= P^4'i^(p)- (133125) 
Moreover, as we shall prove in Section 1121 for |p| < M'"'^ 

(P) = S,,uM'^u, + 0(M^2) ; (I23J26) 

hence, choosing p = or p = {0,M^-^,0) in (I3.5125P it follows that, as 

desired: 

i^i^,h = 0{e^) . ([33127) 
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3.5.3 Ward identity for the effective charge 



Finally, in this Section we shall conclude the analysis of the flow of the running 
coupling constants by showing that the effective charge stays close to the bare one. 
Let us define: 



4'iU.(k,P) 



92 



2 53 



(r-4A|A|) 



?i=0 
lh,K] 



(133128) 



92 



-Agf (0,0) 



<t>=0 



j=4>=o 



where the first and second lines contain respectively the definition of the two and 
three point fermionic Schwinger functions. Proceeding analogously as for the photon 
mass, from (|3.5ll9p we get that, taking one derivative with respect to ap and two 
derivatives with respect to the fermionic external fields. 



^^/^4a2p',u(k,p) 



33129) 



where again the last term is a correction due to the ultraviolet cutoff, which is 
exponentially vanishing as K — >■ +00, and the matrix Q{p) is due to the presence 
of the lattice and it is given by 



Jh,K] 



[h,K] 



Q{p) 



'1 

e-'Pi 



3130) 



therefore, taking the ultraviolet limit in K — ^ +00 in (|3.5[29|) . we find that: 



PM^a^py^k, P) = e [ [Qm ^A':o'Z' (k) " ^'2,^7 + P) [QiP)] ^ 

The Ward identity (I3.5l3ip will allow us to derive the informations that we need 
on the renormalized charges on scale h. By the analysis of Section 14.21 it follows 
that, if k = k' + for |k'| = M'^, |k' + p| < and |p| < (we will be 
interested in the limit p — ^ 0), 



[h,+oo] . 



■,[h,+oo] I 



33131) 



S. 



[h,+oo] 
2,0 

[h,+oo] 



(k) = ^^(l + r2,o(k)), 



2,i,M (k,p) = 

5^ (k' + p) 



32) 



33133) 



h-l- 



z, 



h-1 



eo,fcPor+ + Vh-iei^hP- r+ + p^r2,i,^(k, p; 



with ro,2,^,i/(p), '"2,i,^(k, p), ?'2,o(k) expressed respectively as sums of contributions 
'"a^2,/i,i'(p)5 '"2^i,At(k) P)) '''2fi{^) with N >2, admitting the A^!-bounds: 

I<2,P,.(P)| + \r^,iA^,p)\ + |r2^o(k)| < (const.)^(^)!et • 034) 
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Now, using that 

5(")(k') - 5i")(k' + p) = 5i")(k' + p) [poTt + VH-ip- f^] + p.hfiA^, P) , 

(123135) 

with f2,o,^(k,p) = 0{\p\M~^'^) , we get: 

5Sy~l(k)-5l!;^+°°^(k + p) = -^5('^)(k' + p)Ur°+t;.-ip-fJ5('^)(k') + 

' ^r2,o,^(k,p)+f2,o,;.(k,p)) , (123136) 



with f2,o,/i(k, p) expressed as a sum of contributions f2^Q^^(k, p) bounded as 

kto,;.(k,p)| < M-2''(const.)^(y)!el ; (122137) 

therefore, choosing k' = {M^,0) and p = (p, 0), the Ward identity (I3.5l3ip together 
with the properties p.5135p - (|3.5137p impHes that, taking the Umit p — )■ 0: 



eo 



h = e + Tl [eM^^f 2,0,0 + r2,i,o) := e + eAo,h , (122138) 



with Ao,h a sum of contributions on order N > 2 admitting A^!-bounds. Analogously, 
choosing k' = (0, M^, 0), p = (0,p, 0) and taking the limit p — )■ 0: 

ei,h = e + (evh^iM^'^f 2,0,1 - ^^2,1,1) := e + eAi,h , (122139) 

with A\,h admitting bounds similar to those of Ao,h- This concludes the check of 
our assumptions on the effective couplings on scale h. 

3.6 Asymptotic Lorentz invariance 

In the previous Section we have found that, thanks to the Ward indentities implied 
by the gauge invariance of our model, the flow of the effective couplings can be 
controlled at all orders in renormalized perturbation theory. In this Section we 
shall see that the Ward identities have another remarkable consequence; they imply 
an emergent Lorentz symmetry in our model. 

In fact, consider the identities (j3.5138|) . (|3.5139p : in a Lorentz invariant theory 
we would get that Ao,h = ^i.h- In our case, where Lorentz symmetry is explicitly 
broken by the presence of the lattice, repeating an argument similar to the one 
leading to ()3. 41130 . by the short memory property it follows that: 

eo,h - ei,h = ^i,h(eo,h+i - ei,h+i) + E2,hi^ - ^^h+i) + £^3,h , ([2211) 

where -E'2,/1 and E-^,h can be expressed as series in the couplings on scales > 

/i+ 1, with the A^-th orders bounded proportionally to ("f") M^^ l"^) respectively. 
Therefore, equation (I3.6lip together with the explicit expression (13.41140 for the fixed 
point f_oo of the effective Fermi velocity implies that: 



£-0, — 00 Cl^_oo — 00 (cQ, — 00 61, — 00) ; 



(12212) 
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where -E_oo = 0{e'^^). Equation (|3.6[2p implies that eo,-oo = ei^_-oo at all orders 
in renormalized perturbation theory; plugging this result into the equation ()3.4ll4p 
for ■U—oo it follows that 

which is the speed of light in our units. 



Chapter 4 

The correlation functions 



4.1 Introduction 

In this Chapter we shah adapt the Renormalization Group techniques introduced to 
study the free energy in Chapters [21 [3] to evaluate the generating functionals of the 
Schwinger functions and of the response functions. The only difference with respect 
to the computation of the free energy discussed in the previous two Chapters is the 
presence of the external fields, against which we derive in order to get the desired 
Schwinger or response functions. In the first part of the Chapter, that is in Section 
\4:.2\ we repeat the multiscale analysis for the generating functional of the Schwinger 
functions in presence of an infrared cutoff on scale h* ; the crucial properties (13.51260 , 
(|3.5[32|) required to implement Ward identities at each step of the RG will be proved 
by keeping h* fixed, while our result ()1.4lip on the two point Schwinger function 
will be recovered by taking the limit h* — t- — oo. 

Then, in the second part of the Chapter, namely in Section 14. 3| we will apply 
the multiscale analysis to compute the generating functional of the excitonic, charge 
density wave and density-density correlation functions, introduced in (I1.415P : here 
the analysis will be a bit more subtle, since we will have to introduce new running 
coupling constants associated with monomials containing the external fields, and 
we will have to control their flows. This will be done studying the Beta function, 
or alternatively by using Ward identities. Finally, as an outcome of this procedure 
and performing an explicit computation, we will get the announced results (jl.4[6p . 

4.2 The Schwinger functions 

The multiscale integration used to compute the partition function, described in 
Chapter [3l can be suitably modified in order to compute the generating functional 
yvjj (p) (introduced in Section l3.5.ip . and in particular the two and three- 

point correlation functions appearing in the Ward identities discussed in Section 
13.5.11 From now on we will assume that the ultraviolet limit K — )■ -|-oo has been 
taken. 

Without any loss of generality, let us choose the external fermionic fields as 
^k,a,p = 'Euj=± ^t',a,p,ic^ where 0k',<7,p,a; supported in the infrared region |k'| < oq; 
in this way, the integration of the ultraviolet degrees of freedom is exactly the same 
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reported in Chapter [2] and Appendix |Dj We do so since at the end we wih be 
interested to the scaHng properties of the Schwinger functions for values of the 
fermionic and bosonic external momenta close to the Fermi points and to zero, 
respectively; and clearly, these properties are not affected by the ultraviolet regime. 
Therefore, after the integration of the ultraviolet degrees of freedom we are left with: 

Wf*^+^\j,<t>) ^ ^-mFo j p(^(<0)^p^^,(^^(<0)^gV('')(*(<0),A(<0)+J)+B(vI>(<0),^) ^ 

dMli) 

where the fermionic source term is given by: 

Then, we proceed in a way analogous to the one described in Section I3.2[ we itera- 
tively integrate the fields ^^^\A^^\ . . ., ^C^+i), . . and after the integration 

of the first \h\ infrared scales we are left with a functional integral similar to (|3.2llUp . 
but now involving new terms depending on J, (p. Let us first consider the case h > h*; 
the regime h < h* will be discussed later. 



4.2.1 Multiscale analysis: case h > h* 

We want to inductively prove that 



•e 



33) 



where: S^-^\<t)) is independent of j4, J, contains terms explicitly depending 
on ^4 + J, ^ and of order > 2 in (/;, while B^'^^ is given by: 



+ 



+- 



D 

d 



+ 



where we used the notations (|3.2i4p and 



(3214) 



d 



-V('^)(^,A)<Ak,,.,. 
d 



E 



d 



k',cr,p,a; 



k',a,^^^ + 
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Moreover, the functions Q^^\ G^'^ are defined by the following recursive relations: 



i=h ^^-1 

with QL^^(k') = 1, G^J\k') = 0. Note that, if k' is in the support of gL^\k') 



36) 



G(f)(k') 



(*'^QW(k') + 5c. (k) 



QW(k') = i-z^,.A.^rii£'^+^)(k'), 

that is \\Q^^\\i!) - 1|| < (const.) 4 and \\g^J^\\i!)\\ < (const.) ZJ^^ . 
Clearly, formulas (|4.2i3p -( |4.216|) are true for h = 0, with 



37) 



5[^'+°°](,^) =0, B^^\^,(j),A) = B{^, 



Let us now assume that (|4.2l3p -( |4.2i6p are valid at scales > h, and let us prove that 
the inductive assumption is reproduced at scale h — 1. We proceed as in Section 
3.2[ first, we renormalize the free measure by reabsorbing into P{d'i/^-'^^) the term 
exp{£,i,V('')}, see (l,'12lie)D - (i:i2ll9l) . and then we rescale the fields as in ()3.2l2ip . 
Similarly, in the definition of B^^\ Ea. lK2[^ . we rewrite V^'*) = £*V(^) + V^^\ 
combine the terms proportional to C\i,V^'^^ with those proportional to Q^^~^^\ and 
rewrite 



dk' 

IT 



+ 



+ 



+- 



dk' 
IT 

d 



T 



d 



V('^)(v/z,_i^I/,yl)G(f+i)(kO<^k,.,. 



with defined by ()4.2I6I) . Finally, we rescale by defining 



W}^''i^Zh^^^,cP,A + J) :=W'^''iVZh^,cl),A + J) , 



and perform the integration on scale h: 

g-/3|A|Fh+5(''-i)(9i)+V(''-i){^/z;~r*(-''"^',^'-''"^'+^) . 



R 



■e 



(h-i) 



9) 



where S^^~^\<f)) contains terms depending on (f) but independent of ^'^-'^"^^ \^{<^|'^^) 
and J. Defining S^-^~^^ := 5^'*"^^ + S^-'^\ we immediately see that the inductive 
assumption is reproduced on scale h — 1. 
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4.2.2 Multiscale analysis: case h < h* 

For scales smaller than h* , there are no more bosonic fields to be integrated out, 
and we are left with a purely fermionic theory, with scaling dimensions 3 — 2n, 
2n being the number of external fermionic legs, see Theorem 13.3.11 and following 
lines. Therefore, once that the two-legged subdiagrams have been renormalized 
and step by step reabsorbed into the free fermionic measure, we are left with a 
superrenormalizable theory, as in |36]. In particular, the four fermions interaction 
is irrelevant, while the wave function renormalization and the Fermi velocity are 
modified by a finite amount with respect to their values at h*; that is, if Eh* = 

ma.Xk>h*{\efj.^k\,W,i,k\}- 

Zh = Zh*{l + 0{el)) , Vh = vh*{l + 0{el)) . (112110) 



4.2.3 The fermionic two point function 

In this Section we will provide an explicit formula for the interacting two point 
Schwinger function in presence of a fixed infrared bosonic cutoff; keeping h* finite 
and choosing the external momenta on scale h* we will prove (|3.5[32|) . necessary to 
implement Ward identities in the mutliscale analysis, while taking the limit h* — >■ 
— oo we will prove the announced result (|1.4[ip . 

As for the partition function (see Section 13. 2|) , the kernels of the effective poten- 
tials produced by the multiscale integration of wjj (p) can be represented 
as sums over trees, which in turn can be evaluated as sums over Feynman graphs. 
Let us consider first the expansion for the 2-point Schwinger function. After hav- 
ing taken functional derivatives with respect to (p^, (f)^, and after having set 

(h*) 

J = (/> = 0, we get an expansion in terms of a new class of trees r S with 

k £ {—oo, —1] the scale of the root and h > k; these trees are similar to the ones 
described in Chapter [21 up to the following differences. 

1. There are -|- 2 end-points and two of them, called vi,V2, are special and, 
respectively, correspond to 

Qif--^)(ko^*g;;r'^" or to ^'g;:r'^^''^Q^^^'^(k'). 

2. The first vertex whose cluster contains both vi, V2, denoted by v, is on scale 
h. No TZ operation is associated to the vertices on the line joining v to the 
root. 

3. There are no lines external to the cluster corresponding to the root. 

4. There are no bosonic lines external to clusters on scale h < h* . 

In terms of the new trees we can expand the two point Schwinger function as, taking 
k = k' + p% with |k'| < oq: 



P\A\ 
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where /i^' < is the integer such that M^^' < |k'| < M'^k'+i^ and 52(r; k) is defined 
in a way similar to V^^\t) in (I3.3l4p . modulo the modifications described in items 
(l)-(4) above. Using the bounds described immediately after (|4.217|) . which are valid 
for k' belonging to the support of ^^^^(k'), and proceeding as in Section [3^ we get 
bounds on S2{t; k), which are the analogues of Theorem I3.3.1t 

_E _E E l|52(r;k)||<(const.)^et.(y)!^- (52112) 

k,h,N 

Being the bound (j4.2[12|) uniform in h* , our result (|1.4iip on the two point Schwinger 
function is obtained after taking the limit h* — >■ — oo in (j4.2[ll|) . and taking into 
account that Zj, Vj can be replaced with their interpolations Z(k') ~ |k'|~'', f (k') ~ 
1 — (1— i')|k'|'', provided that the error term -R(k') in the main result (|1.4iip is defined 
so to take into account the corrections O(e^logM) generated by the replacements; 
see discussion after (|3.4119p . 

In order to understand (|4.2U2p . it is enough to notice that, as far as dimensional 
bounds are concerned, the vertices vi and V2 play the role of two ly vertices with 
an external line (the cp line) and an extra y^M"'*'" factor each. Moreover, 
since the vertices on the path Vr,v connecting the root with v are not associated 
with any TZ operation, we need to multiply the value of the tree r G T^'^ by 
j^{i/2){h-k) j^£{i/2){k-h) ^ and to exploit the factor M^^/^)^'^"'^) in order to renormalize 
all the clusters in Vr^v- Therefore, 

h-i ^ 
_E _E E l|52(r;k)||<(const.)^(^)!. 

h=—oo k=—oo t-ctS^_ ) 
k,h,N 

_ _ 

. ^* ^ ^ M''M''"''k'M(^/2)(''"'')M"2''i" (I12J13) 

h<hy k<h 

where: the factor is due to the fact that graphs associated to the trees r £ Tr^ L 

^ k,n,N 

have two external lines; the factor M^~^^' is given by the product of the two short 
memory factors associated to the two paths connecting v with vi and V2, respectively; 
the "bad" factor M^mh-k) jg ^he price to pay to renormalize the vertices in T^r,v'i 
the Z^J^ and the last M~^^k' are due to the fact that vi, V2 behave dimensionally 

as V vertices times an extra Z, ^^"^M^^^' factor. Performing the summation over 

k and h in ()4.2[13p . we get ()4.2[12p . Note also that, if k' and k' + p are on scale 
^k' — h* , then the derivatives of ||52(T;k)|| can be dimensionally bounded as 

_E _E E l|9kcS2(r;k)||<(const.)^et>(y)! ^ ^ , (02114) 

k,h,N 

from which the bound on r2^^^{\i,p) stated in (|3.5137p immediately follows. 
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4.2.4 The bosonic two point function 

In this Section we will prove the property (|3.5[26|) . which is necessary to extract 
informations on the effective photon mass on scale h from the Ward identity (I3.5124P . 

Taking functional derivatives with respect to J/x.pi Ju-v ^-iid setting J = = 0, 
we get an expansion in terms of a new class of trees r G "^^J^ with k < h* the 
scale of the root; these trees are similar to those described in Chapter [21 up to the 
following differences. 

1. There are N endpoints and three of them, called vi, V2, ws, are special. The 
endpoints vi, V2 correspond respectively to 

Zh,,-i e,,H ji:i''^ - M^^iz.^,,,^4,p , (132115) 

instead, corresponds to —M^^su^ f^^^S^ j^ (the only non-vanishing tree this 
endpoint can be atteched vertex is trivial, that is the vertex is attached directly 
to the root, and it has root scale h*). 

1. There are no lines external to the cluster corresponding to the root. 

3. There are no bosonic lines external to the clusters on scales h < h*. 

Pick IpI < in terms of these trees, the bosonic 2-point function can be 

written as 

. . Wg;'+-'(J,0) = -M^V,,,.^,,. + E '52(r;p) 

^ ' k,N 

(132116) 

and 52(t, p) is defined in a way similar to V^^\t), modulo the modifications de- 
scribed in items (1) - (4) above. Proceeding as in Section [3^31 we get that 

E E ||52(r;p)||<(const.)^e^*(|)!M'^\ (j32117) 

^ ' k,N 

To understand formula (|4.2117p . notice first that now all the vertices are renormal- 
ized; this is a simple consequence of the fact that the J field is completely equivalent 
to a ^4 field, because in the multiscale integration the fields A, J always appear in 
the combination A + J . Then, from the bounds derived in Section [3.31 it is easy to 
see that 

r G ||52(t;p)|| < (const.)^ef (y)!M^ , ([3318) 

and this concludes the proof of (|4.2ll7p . 

4.2.5 The three point function 

Finally, let us consider the three point Schwinger function; in this Section we shall 
prove the property (j3.5l33p . which is needed to derive the Ward identity for the 
effective charge on a given scale. 
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Let us pick k = k' + p'J^, with |k'| = |k' + p| < M^" and |p| < M^* , which 

is the condition required in order to apply Ward Identities in the form described 
in Section 13.51 The expansion of 3-point function is very similar to the one just 
described above for the 2-point function. The result can be written in the form 

(/3^iai^a) f = mm 

eM4Gr-'Hk' + p)fOr)(k')Qr)(k')+ E E E '53(r;k,p) 

h<h* hvg>h* k,h,hy.j,N 

where Tt'I ? is a new class of trees, with A; < the scale of the root, similar 

to the trees in 7^^^ ^ , up to the fact that they have + 3 endpoints rather than 
+ 2 (see item (1) in the list preceding (|4.2[11|) ): three of them are special: vi and 
V2 are associated to the same contributions described in item (1) above, while 
is associated to a contribution Zh-^^ie^^h^j'^^-^""'^^ - M^'^su^^h-^A^^p, with v^, the 
vertex immediately preceding vj, on r (which the endpoint V'i is attached to) and 
/i„3 > h* . The value of the tree, 53(t; k, p), is defined in a way similar to 52(t; k), 
modulo the modifications described above. 53(r;k, p) admits bounds analogous to 
(jMUl-dlZISI); recalling that |k'| = |k' + p| < M^" and |p| < M'**, we find: 

_E _E E E l|53(r;k,p)||<(const.)^(^)!et.^- 

h=—oo k=—oo hv2,=h*-\-l ^^j-i^* ) 

. J2 M^'^/'^^^^-^^ m'^-^* M^'^/'^^^^'-^^'s'i M-'^^* , (02120) 
h<h* 

k<h 
hv^>h* 

where: M^^/^^^'^"'*) is the short memory factor associated to the path between the 
root and v; M^~^* is the product of the two short memory factors associated to the 
paths connecting v with vi and U2, respectively; M*^^/^'*^'**~^"3) is the short memory 
factor associated to a path between h* and V'^; M~^'* is the product of two 

factors M~^' Zj^J_^ associated to the vertices vi and V2 (see the discussion following 
(j4.2il2p and recall that in this case /ik' = h*). We remark that in this case, contrary 
to the case of the 2-point function, the fact that there is no TZ operator acting on 
the vertices on the path between the root and v does not create any problem, since 
those vertices are automatically irrelevant (they behave as vertices with at least 
5 external lines, i.e., J, 4>~ , (j)'^ and at least two fermionic lines) and, therefore, 
7^ = 1 on them. Note also that the vertices of type J4>i^, which have an TZ operator 
acting on, can only be on scale /i* — 1 or h* (by conservation of momentum) and, 
therefore, the action of the TZ operator on such vertices automatically gives the 
usual dimensional gain of the form (const.) M^^~^-"' . Performing the summations 
over fe, h, hy^ in (|4.2[2U|) . we find the analogue of (|4.2ll2p : 



E l|53(r;k,p)||<(const.)^(^)!4-^^^_^ 

h=—OQ k=—oo hvQ=h*-\~l 

k,h,hv.^ ,N 



(112121) 
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from which the bound on (k, p) stated in ()3.5l34p foUows. 

4.3 Response functions 

The techniques introduced to evaluate the free energy and the Schwinger functions 
cas be used to compute the response functions^ or generalized susceptibilities; phys- 
icahy, these objects provide further understanding on which is the effect of the 
electromagnetic interaction on the macroscopic properties of the model, in particu- 
lar they suggest what possible physical instabilities are enhanced by the interaction, 
|78j . As for the two pont Schwinger functions, we shall see that the response func- 
tions as well decay with interaction dependent anomalous exponents. 
Let us define the imaginary time response functions C^'^\'x. — y) as follows: 











^ \ ^''^ x+(0,5j),(T 

(T=t4- 


x,j) it 6"*" , p , a„ „e 

x+(0,5j),<7 ^''^ 


(CDW) 
^x,j 


:= ajr„a~_ — 5"*" , -, ^ 

^ x,(T x,o- x+{0,<5j) 

(T=t4- 


cr^X+(0,5j),o- ' 


(D) 


^ x,cr x,cr x+{0,<5j) 

(T=t4- 


(T^X+(0,5j),O- ' 



1) 



-ie^{xj) 



Being p'^^^\ p^^^'^^ and p^^^ gauge-invariant operators, we can represent their 
grand-canonical correlation functions using a functional integral representation in 
the Feynman gauge. Therefore, 



Cg)(x-y) = -^|^log/p(dM/)P(d^)e^(*'^H^(*'*'^) 



x,« yj 



$=0 



x,« y,j 



a2 _ 



$=0 



^A) := / dxoE E (5312) 

where = (/3^a)"-^ EqeP^ ^ e"'''^*q"i (we choose l>|j"] G M for simpHcity), and 
p'^'^ are given by (j4.3ll|) with , 6^ , A replaced by fermionic and bosonic Gaussian 
fields. 



4.3.1 Renormalization group analysis 

All these generalized susceptibilities can be evaluated using multiscale analysis and 
renormalization group. As we are going to see, the external fields ^^"^ behave 
dimensionally as A fields; the main difference with respect to the analyses of Chapter 
13.31 and Section 14.21 is that now we will have to introduce new running coupling 
constants, associated with the relevant and marginal (^-dependent monomials. The 
flow of some of these new running coupling constants, namely those associated 
with the monomials will diverge as a power laws in the infrared, in a 
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way similar to the wave function renormalization but in general with anomalous 
exponents different from r] = + . . .. A non- vanishing difference between these 
new anomalous exponents and r] will give rise to the anomalous scaling of the 
response functions. 



Multiscale analysis 

The ultraviolet degrees of freedom can be integrated out by following a procedure 
completely analogous to the one described in Section [2^ and AppendixlDl It follows 
that: 

e^/S'^^*) = MS) 
g-/3|A|Fo+C(>o)(#) J p(^^(<0)^p(^^^(<0)^gV(0)(vI'(<o),A(<o))+©W(<I',*(^o),A(<o)) ^ 

where all the quantities appearing in (I4.3I3|) have been defined Section 13. 2| except 
for C^-^\^) and D^^^ which contain <I>-dependent terms. In particular, 

p(°)($,^,^)= V y ff]K K n^MP - 

^ ' ' ^ /-^ / 11 k' ,(Ti,p2i-l,I^2i-l k' ,o-i,p2i,I^2i 11 >^"P^ 

^ -I rr /,! n 1 1 



q>l 1=1 i=l 

n+m>l tJ,,a,j 



r=l \i=l r=l i=l j 

c^^-'^ (^) = E / n ({qJ)^ ft q^l > »4) 

g>l"' r=l \j=l / 

where the kernels are computed following the rules explained in Appendix |D] and 
Section 13. 3t by taking into account the presence of new vertices, corresponding to 
the source terms. It is easy to see that the external fields behave dimensionally 
as A fields; therefore, the kernels proportional to the monomials ^^"'^A behave as 
relevant ones, while the ones proportional to <I)(")vi/+\i/^^ ^^"'^AA, are 
marginal. 

As for the effective potential V'-'''*, we split the ^-dependent part of the interac- 
tion as CD^^^ +TZ'D^'^\ with TZ = 1 — £; the localization operator C acts linearly on 
the kernels of "D^^^ in the following way: 

^wiSC:(p,q) := W^;?f;f(0,0)=:C2^, 

'C^o,i,2,/7(P> q) := ^o,i,2,/7(0>0) =:«M,o • (3315) 

and CW2nm q ~ ^ otherwise. As for the case $ = 0, the number of running coupling 
constants will be reduced by the symmetry properties listed in the following Lemma, 
which will play the same role of Lemma 12.3.11 in the analysis of the free energy. 

Lemma 4.3.1 The Gaussian integrations P{d^), P(dA), the interaction V{'^,A) 
and the source term D{^,'^,A) are invariant under the transformations (1) - (8) 



2n 
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of Lemma \2.3.1{ provided that transforms in the following way (remember that 
1' = 1, 2' = 3, 3' = 2;.- 

l-g ^ l>g„ with q = (go, gi, -92); 

Proof. The proof is a straightforward extension of the one of Lemma I2.3.H we leave 
the details to the reader. ■ 

By symmetry (4) it follows that: 

(11316) 

51)(P2-1)}- 

while by symmetries (4) - (8) we get: 

^^2,0,1,0; l^^i Pf ~ Pf i — • ^tii,0lQ X ) ~ iii'O ' 

with Z°o = -^-^,0. ^^,0 e and Z^q = 1 + 0(6^). Therefore, from (lOTUD and 
(|¥:HT71) we find that 

f dk' (jq i^(a),|,(<0)+ „ pM*(<0)- 



'^'^2,0,l,P,t£V"' Pf Pf J — 

exp \ip'p{52 - h){pi - 1) - '^P^F^^'i ~ 

•^•^2,0,1,^,0; V"' Pf Pf J ' 
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where the matrices r^"j are easily obtained from ()4.3I6|) . (I4.317P : 



1 



Repeating the analysis described in detail in Chapter [21 we can integrate the 
fields ip^°\A^^\ ^(^), A(^), in an interative way; after the inte- 

gration of the first \h\ infrared scales we get that: 

gW^,L(*) = g-/3|A|F^+C(>'')(1>) . (113110) 



where all the objects appearing in ()4.3[10p have been defined in Section lS^ except for 
Ci>h)^ Dih)^ which are given by (|4.3[4p with replaced by /i, and the kernels w!^^^ 
will be defined inductively. We proceed as in Section [3.21 first, we renormalize the 
free measure by reabsorbing into P^dt^^-^^) the term exp{£,i,V'''^'*}, see (|3.2ll6p - 
(I3.2[19p . and then we rescale the fields as in (I3.2l2ip . Similarly, we rewrite 



V^^\^, V^^, A) =: JZh-i^, A) , ([13111) 



and we split T)^ as CD^^^ + TZ'D^^''\ where the action of C on the kernels of p(^) is 
defined as follows: 

^<o!';SJk,q) := <o!5JO,p^^^-p^^), 

^wi'^£lip,^) ■■= <2!'";^(o,o)=:c;^ 

^^o,i,2,M(P'q) •= ^o,i,2,m(0'0) =: • ([13112) 

Moreover, the symmetry properties (1) - (8) of Lemma [4.3.11 which are preserved 
by the multiscale integration, imply the analogous of (I4.3l8p . that is: 

2^ J 2) (27r)3 'I'J '^'+q'<^''"i''^i ''''<^'/'2'^2 2,0,l,t;iV"'Pi^ Pf ) 

-J ^ (27r)3 k'+q,<7,c.i^ai,/iiiiij'^k',<7,a;2 ' (|4^13j 

where Z^f^ = Z'^^y^, Z^^^ E M, and the matrices T^"] are defined in ([CT9]) . At this 
point, we are ready to integrate the scale h; proceeding as after ([3.2125p . defining 

gV(^-i)(^/z;:i7*,A)+I'('>-i)($,^/z;~r*,A)-/3|A|Ffc-C('>-i)(*) ._ ([13114) 
/ p^(^d^ih)^p^(^dA^h)yVW{^/Z^,{-i+^W),A+AW)+V^^^ 
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our inductive assumption in reproduced at the scale h — 1 with C^-^'""*^^ = C^-^'^ + 
The integration is performed using the Gallavotti - Nicolo tree expansion 
described in Section the only difference with respect to what has been discussed 
in Section [3.31 is that the resulting trees have new endpoints, corresponding to the 
various contributions to CD^^\ 

The outcome of this procedure is that the kernels of C^-^~^) are expressed in 
terms of a renormalized expansion which involves the old effective parameters plus 
the new ones {i^'^'l, i^'^'lk^ C^fc' ^/I'fc' ■^Zk\k>h'i to close our renormalization group 
analysis we have to control the flow of these quantities. 

The flow of the running coupling constants 

In this Section we show that the flows of all the new effective parameters except ^ 
are bounded; in particular, all the new running coupling constants will be bounded 
by 0(e) uniformly in h. To see this, we shall use Ward identities (for i/^'^, C^'^, 

«;~'^) and symmetry considerations (for v'^'lf^)- Regarding .^"/j, these quantities in 
general will behave as M~^°'^, that is their evolutions are governed by anomalous 
exponents^ which will be explicitly computed at lowest order; these computations 
will allow us to characterize the long distance behavior of the excitonic, charge 
density wave and density-density susceptibilities. 

The flows of t'^'/i, C^'^, As for v^^h, the "mass" of the photon field, the flow 

of these running coupling constants will be controlled implementing Ward identities 
at each RG step. Notice that, as for z/^,hi the Beta function does not provide any 
help here; for instance, the evolution equation of f^'^ is given by 

<i = + , (133115) 

and a single non-vanishing order in /3p^_,„^ would generically produce an exponen- 
tially divergent flow for f^'^ in the limit /i — )■ — oo. Similar behaviors would take 
place for C^;^, kJ^. 

The Ward identities that we need are implied the gauge invariance of our theory 
plus the fact that the source terms appearing in the generating functional Y^p^ii^) 
are separately invariant under local gauge transformations. We will proceed as for 
^fj.,h-, ^fj.,h, see Section 1231 we introduce a reference model which is regularized by a 
bosonic infrared cut-off on scale h, and we derive Ward identities for this new theory, 
exploiting the facts that: (i) the running coupling constants on scales greater than 
the one of the infrared cut-off are the same of the non-regularized model; (ii) the 
infrared cut-off does not break gauge invariance, and does not produces corrections 
to the WI; (iii) the new model can be studied repeating the multiscale integration 
described in Section 13. 5| with the only difference that after the scale h the theory 
becomes superrenormalizable, which in particular implies that the running coupling 
constants "cease to flow" on scale h. See Section 14.21 for the discussion of the 
multiscale integration in presence of a bosonic infrared cut-off in the case of the 
generating functional of the Schwinger functions; a completely analogous analysis 
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can be repeated in presence of external fields. Let: 



:= J P(dvI/)P>;,(d^)e^(*'^+-^)+^(*'*'^+-^) , 



(113116) 



where the subscript > h denotes the presence of an infrared bosonic cutoff on scale 
h, to be imposed as in (I3.5l3p . We assume that the ultraviolet limit K — >■ +oo has 
been taken; the presence of a finite u.v. cutoff can be discussed as in Appendix [Gl 
and produces corrections that vanish exponentially in the limit of cutoff removal. 
The invariance of the theory under local gauge transformations implies that 







$=J=o=0 



[h,+oo] 



<I>=J=0 



17) 



4>=J=o=0 



KQu +Pu) 



Q3 



[h,+oo] 
/3,L 



q-p 



<S>=J=0 



q3 



KQu + Pu 



^3 



$=J=Q=0 



q J " p,j' '^•^1',-q-p 



wLY°'i(a>,j) 



<I>=J=0 



(03118) 



At the same time we know that, repeating an analysis completely analogous to 
the one performed for e^^/i and t'^^/i, choosing — 

^ij.,uq, Pu = Sti,uP, such that 

\p\ < M^~^, \q\ < M^~^, \q + p\< M^~^, 



'($,j) 



^3 



5^ T-ry^.+oo] 



-wLY°^i($,j) 



$=J=0 



<I'=J=0 



^'K + 0(e^)] , 



q J p,j' '^'^A'.-q-p 



$=j=o 

putting together (|4.3[19|) with the Ward identities (|4.3ll7p we get that 
as desired. 



(113119) 



20) 



The flow of f^'^^. These running coupling constants cannot be studied using 
Ward identities, at least not for all /i, z^; therefore, to prove the boundedness of 
their flows we will rely on simple symmetry considerations. The evolution equation 



for f^'^ ^ has the form: 



o,7 
^ 1, 



a,j na,a 
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where, setting et := maxfc>h{e^,fc, iyu,k, i^^^'fc, 1"^,';^,^' C'fc' (^^^ is intended 

over all the labels), the Beta function appearing in (|4.3[21|) is expressed as a renor- 
maUzed series in the running coupling constants with the A^-th order bounded pro- 
portionally to {N/2)\e^ . However, this dimensional bound is too pessimistic in 
this case; in fact, the graphs contributing to contain an even number of 

fermion propagators, and one propagator (fermionic or bosonic) is derived with re- 
spect to the /x-th component of the momentum. Therefore, rewriting the fermion 
propagator as ^^^^(k') = (f^'^^^ii^) -hffif^(k'), where g^^l^ii^) is obtained setting 
Rk-i,Lo(y) = in ()3.2I26I) (that is, it is the linear part of the single scale propaga- 
tor) and ||g(^)(k)|| < (const. )M'^'||(7('^)(k)||, we decompose the Beta function as 

where: 

(i) the terms labelled by (re/) are obtained by considering only trees with end- 
points on scales < 0, and replacing all fermion propagators with their rela- 
tivistic counterparts; 

(ii) because of the improved dimensional bound on ^ (k) with respect to the one 
on cju} (k) and of the short memory property of the Gallavotti-Nicolo trees, 
the bounds on the terms denoted by a tilde are improved by a factor M^^ 
with respect to the basic ones. 

Then, since g^^\^ii}^) = ~ff(!fre/(~^)' w'^'^\'P) = ''^^^''("P)) the parity of the graphs 
contributing to the relativistic part of the Beta function is equal to —1; this means 
that 

= , 023) 

and the equation (I4.3l2ip can be safely iterated up to the scale 0. This shows that 
^'jj.'ih ~ ^(^)' uniformly in /i, as desired. 

The flow of Finally, we conclude the discussion of the flow of the new 

running coupling constants by analyzing the flow of Z^^] these are marginal, and 
their flows are very similar to the one of the wave function renormalization Z/^, see 
Section 13.4.21 In fact, these running coupling constants evolve according to the 
following flow equation: 

Z^,h = Z^,h+i (1 + Plh+i,i) + C,h+i,2 , (133124) 

where ^j^n are expressed as renormalized series admitting A^! bounds; in partic- 
ular, 13'^^.^ = O(e^), P^h2 — 0{e'^). The graphs contributing to /3°/ji have the 
external field emerging from a vertex of type while the ones appearing 

in the value of /3° ^ 2 have the external field emerging from a vertex of type ^^"''^AA; 
see Fig. 14.11 for an example of two such graphs. 

However, as anticipated in the caption of Fig. 14.11 the bound on /3" 2 admits 
an improvement of a factor M^'^ with respect to the basic one. To see this, simply 
decompose P^}i^i2 /^"1+12 + '^"h+i2' exactly as for the flows of the other 
marginal terms (see discussion after (|4.3l22p ): again, /3"'^'^|2 = because it is 
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Figure 4.1. a) A graph contributing to ^^^f3" ^; b) a graph contributing to /3"/j2- 
Tire grapli b) is of order e^M^^; in fact, the same graph obtained replacing the fermionic 
propagator with its relativistic approximation gl'^j (k) is vanishing by parity. 



given by a sum of graphs with an odd number of relativistic fermion propagator, 
while /32 

h+1,2 ^ dimensional bound improved by M . 
Hence, Eq. (|4.3[24|) can be iterated up to the scale and it follows that 

Z^,h - M-^'^'^ , as /i ^ -cx) , (03125) 
rjZ := lim logM {l + PZ,h+i,i) (133126) 

— ft— 5- — OO — 

the lowest order contribution to r/" are computed in Appendix IF.4| the results 
are summarized in Table UTT] (recall that by symmetry = Z"^ ^). 



^\ — 
a 


(+,-) 


(+,+) 


E+ 








e2 


e2 


CDW 






D 




e2 
127r^ 



Table 4.1. The lowest order contributions 77S to the anomalous exponents ri[ 



Therefore, the conclusion is that the most divergent are those corresponding 
to a = £"+, uji = — a;2 and a = CDW, uJi = UJ2; in particular, these four running 
coupling constants diverge with anomalous exponents bigger than the one of Zh. As 
we are going to see in the next Section, this implies that the spatial decays of the 
excitonic and charge density wave susceptibilities are depressed by the interaction. 

4.3.2 Computation of the response functions 

Finally, we have now all the elements that we need to explicitly compute the response 
function; in particular, we shall prove the results (jl.4[6p . (|1.417|) . It follows that: 
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where C^^\^) has been defined in (j4.3ll4p : the lowest order contribution to C^^ 
is represented in Fig. 14.21 




Figure 4.2. The lowest order graph contributing to C''* ^■'($). A sum over all the possible 
inner scales is understood. 



Therefore, taking i = j = 1 for simplicity, we can rewrite ()4.3l27p as follows: 

cfi(x) = -2'nYr°5("-^'-)(x)r°5("-''-)(-x)) m28i 



-2 E E p^^-(r"5£'(x)r"^;)(-x))63l29) 



h,h'<0 ^i,a;2=± ^h-lZh'^l 
h<l 



where the factor 2 is due to the sum over the spins, h \/ h' := max{h,h'}, and 



r2'g^]^(x — y) contains higher order terms. We will get that, for |x| ^ 1: 



|x| ^++ |x| 

where, for 1 > > 1/2: 

C ■■= 2(^2 - V) , rf,i(x) = r^^fCx) + E <'r(x) , »32) 

n>2 

l<f(x)l < , ,4lTlle > K'rWI < (const.)" (f)!e" '^^^ , 



X -^ii X 



and, if 6 = 87r727: 



2 2 2 

^ VXj---^> ^2 W - ^|-p , 

Gf Mx) = ^:i^^|±^ , G^(x) = -l, 

Gf - (X) = ^^^gf^ , G2^- (X) = . (S3133) 
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Remark 4 1. Notice that by symmetry (4) of Lemma \4-'J-1 we know that 

C^i,,+i(x) = q:/(xo,Tf)) ; (113134) 



the correlations Cfj(x) with i ^ j can he studied in exactly the same way as 
Cf ]^(x), and we will not discuss them here. 



2. From the explicit lowest order computations of rj, see Table 4-1 <ind Ap- 
pendix F.4 , we know that all C° = + O(e^) except for , which are 
given by: 



4e^ 
3^ 



3135) 



therefore, we conclude that the decays of the correlations Cf]^(x) are depressed 
by the presence of the electromagnetic interaction; in the spirit of 178^ . this 
suggests that the electromagnetic interaction favors instabilities of excitonic 
or charge density wave type. 



We will study the three contributions (|4.3l28p . (|4.3l29p . (|4.3l3Up separately; the 
following bounds will play a crucial role: 



T\j{2+no+ni+n2)h 
10 1 2 y<^ V ;|| - (M^|x|)^ 

/\/f (2+no+ni +n2 )fe 



(113136) 

for all /i < 1 and > , 



and with gj^\x.) := (7("''"')(x), ■u;(^)(x) = w^^''"'\'x.). The proof of these bounds for 
rij = is a consequence of the following inequalities: 



Ibi'^x) 
^i^)(x) 



< 



/^5if^(k')||<C^M2\ (113137) 
^|^e-^'-<05-^5-^^('^)(k': 



3138) 

from which (|4.3136p immediatly follow^. The case Ui can be derived in exactly 
the same way; the details are left to the reader. 



Bound for (14.3L28|> . This term can be estimated simply by using the bound 
(|4.3l36p : it follows that 



rj^/pa^(«.^;.)(x)r"5("-^-)(-x; 



< 



Cn 



^l + |x|^)2 



3139) 



for ah > 0. 



^For A'^ even use that |x|^ — (X^^_g xf) 2 and expand the square; for N odd use that |x| 
x||x|"'^~"'", that |x| < \xi\ and the fact that — 1 is even. 
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Evaluation of (j4.3l29p . This term can be evaluated explicitly, asymptotically in 
|x| —7- +00 and at leading order in e. We will rewrite ()4.3[29p as a quantity which 
can be explicitly evaluated, plus error terms which will be small, either because they 
decay faster in |x| or because they decay with the same power but are O(e^). 

Error terms. First, we will replace the ratio appearing in (j4.3[29|) with a 
function of |x|, and we will show that the bound on the remainder is improved by 
a factor O(e^); then, we will compute the leading contribution to (|4.3l29p explicitly. 
Let us rewrite: 

= (113140) 

= |x|2(^2-';)[i + (|x|Af'^v'^V"'^(|x|M'')''(|x|M^')''^^|^ - 1 

where = 1 + O(e^). Plugging (j4.3[40p . the contribution due to the error term 
SLL},h,h'{^) can be bounded using 

(i) the inequality r° — 1 < |Qlogr|(r" + r~°), valid for all a G M and r > 0; 

(ii) the bound (j4.3l39p for the single scale propagator; 

(iii) the faclH that, for some C^^a > and > a > 0: 

Using (i), (ii), (iii), is is easy to see that 

Y: ..,,,,Kx)Tr(r"gW(x)r°<7i.';^(-x)) = E ^S'^^Cx) , 

h,h'<0 N>2 

|f-^(x)| < (const.)^(f)!e^ ^^^J(^_^.) • S42) 
Then, let us rewrite 

9^^H^') = ^Skk') + Ckk') + sS(ko 

Ci(k') := ^i-^lK + k'rj, (113143) 

where: (i) s^^\{'k') contains the error committed replacing respectively i)/i_i(k'), 

A(kO with 1, A(kO in g!^\k'y, (ii) sl^ik') is equal to gl%i^k')R,,_^^U^k'). There- 
fore, 

pi;!(k')|| < (const.)e^||ff(f)(kO|| , ll^kkOII < (const.) |k'|||5i^)(k')|| . 

(101 44) 



|x|2('?2-'?) 



^The proof of (|4.3I41|) goes as follows. Let /i(x) := —[log |x|], where [.] denotes the integer part; 
split the sum X]^<;o in X/h<h(x) ~'~X/h>h(x)' ^^^^ sum can be bounded replacing the argument 
of the sum with Ix]"^'^, while the second sum is bounded replacing the argument of the sum with 
(M'-lxr. 
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Correspondingly, we rewrite (x) as 5^ ^(x) + {(x) + li^)^ where 



5S]:(x) := e-^^"^-J^e-^^'-gi;:},{k'), (03145) 



s^j(^ := e-P--y — e-'^-^-Kk'); (0346) 

The contribution of sjfi(x) to (14.31:^9^ can be bounded exactly as the corresponding 
to s/i,a;.h,ft,'(x)) and we get a bound like (|4.3[42|) : the contribution of s^^^i'^) to (|4.3129|) 
is bounded using (I4.3l4ip . ()4.3l44p . and it decays faster than \x\~^~^'^^^~^^\ 

Evaluation of the leading term. Therefore, the leading contribution to 
(lOngi) is 

Ci"i,^(x) := -2 |x|'(^^^'^)TY(r-<7c.„£(x)r"g^,,£(-x)) , (03147) 

where g^^ci'^) •= J2h<o 9^c^ explicit computation shows that, see Appendix 

EH ' 

where the last term denotes corrections which vanish as |x| — t- +oo faster than any 
power. Therefore, setting Q = 2(ry2 — ??), b := 6'^/8 = 87r^/27 we get that, see 
Appendix IF.71 

1^ I ^1 2 ™2 I ^2 ™2 

'-'l.l.^VXJ - — — ^ZTd h COSI^Pp, • , , ,f;_<;D 



1 ™2 _|_ ™2 



CDVK^ 

4_cCDW "I" ^OS^_/^^ • ,f;_cCDlV i 



6|xp-'^++ 6|x| 

E+ , . _ -xj+xl _ COs(p^ • X) 

6|x|6-^++ 6|x|^-«+- 

cfr^Cx) = ^i^^i±M + cos(p-+-f)^o^4^. (03149) 
6|x|6-«++ 6|x|^-«+- 

from which ()4.3l33p follows. 

Bound for 14.31301 Finally, we have to show that all the other Feynman graphs 
give a smaller contribution. This can be done repeating an analysis similar to the 
one of [S], up to a few differences that will be discussed here. We can have three 
situations: 

(i) all the external fields emerge from vertices of type '^^"-'^^i^j^'tJj ^^-^'^ graph 
has more than one inner loop; 

(ii) at least one external field emerges from a vertex of type ^^'^-'AA; 
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(iii) at least one external field emerges from a vertex of type ^^"^A. 

The vertices of type ${")$(")^ do not give any contribution, since the corre- 
sponding running coupling constants n^'j^'^ are vanishing by symmetries (1) - (8) of 
Lemma 14.3.11 

The case (i) can be worked out using that the minimal loop connecting the 
points x and y is formed by at least two single scale propagators, which satisfy 
the bound (I4.3136P : then, we can apply the analysis of [8], that will be not re- 
produced here. The A^-th order contribution due to these graphs is bounded by 
(const.)^(Af/2)!e^|xr^+«S, N>2. 

Similarly, using that the running coupling constants associated to vertices of 
type ^^"'^AA are 0{e'^), we get that in case (ii) the A^-th order contribution is 
bounded by (const.)^(iV/2)!e^|xr^+(^/2) max^^. 

Let us consider case (iii); recall that the vertices of type ^^"^A correspond to 
M^v^'l^ or to Pui^'^'lh, and both the running coupling constants are 0(e). Consider 
the simplest of these graphs, namely the one depicted in Fig. 14. 3[ Assume that 



0« X 

^A/WWWW 

h 



Figure 4.3. In momentum space, the vertices of the graph correspond to M^v'^'j^ or to 

a,l 

both the vertices correspond to Vv^'^^h'^ the value of the graph in Fig. 14.31 is 



Using the bounds (I4.3I36|) . (I4.3[41|) we get that 

EQ,1 Q,l 
^IJ,,u,h^fj.y,h f)^ 

" h< 



50) 



'u;('^)(x) 



< 



(const. )e2 



X 
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and the same result holds if one vertex is replaced by M^u'^'^] higher order graphs 
can be bounded repeating the analysis of and the result is that the generic graph 
of order N of type (ii) is bounded by (const. )^e^|x|~^'''^"^^^^ 



Chapter 5 

Lattice distortions 



5.1 Introduction 

In this Chapter we wih discuss the effects of lattice distortions on the model that we 
introduced and studied in the previous Chapters. In general [8^ [55] . small lattice 
distortions can be taken into account by considering a position and bond dependent 
hopping parameter, of the form t^^i := t + (p^^i, where (p^^i is proportional to the 
distortion of the bond connecting the site x to its nearest neighbour x + 6i. Again, 
we shall perform a RG analysis to characterize the low energy behavior of the model; 
as we are going to see, some special lattice distortions, the Kekule ones, will be dra- 
matically enhanced by the electromagnetic interaction. Interestingly, the amplitude 
of this distortion behaves as a bare mass for the fermion progagator. The presence 
of a non- vanishing initial distortion will give rise to a new running coupling constant 
with respect to the case studied in Chapter [3l namely the effective fermionic mass 
Ah. We shall distinguish two regimes: in the first regime the effective mass will grow 
from its bare value Aq to a dressed value A ^ Aq; in this regime the fermionic 
field will behave essentially as a massless field, and the analysis will be qualitatively 
the same as in the case Aq = 0. After the scale h such that A^ = A, the fermion 
propagator will become massive, since the typical momentum scale will be smaller 
than the mass A; as a consequence, the flow of the running coupling constants will 
remain close to the values reached at the scale h. 

After this, we will investigate the possibility of a spontaneous deformation of the 
honeycomb lattice; in particular, we will write a non-BCS self-consistence equation 
of the amplitude of the distortion, that is for the mass of the fermion field. This 
equation is similar to the one found by Mastropietro in the context of Luttinger 
superconductors, see |62]. The equation will be discussed from a qualitative point 
of view; as we shall see, already for small coupling the solution will be influenced in 
a non-trivial way by the electromagnetic interaction. For strong coupling, we will 
argue that the phenomenon of mass generation is favored by the interaction. 

The most direct physical interest of our study is connected with the possibility of 
a metal-insulator transition in graphene, an (yet open) issue on which great efforts, 
both theoretical and experimental, have been and are still devoted. However, there 
is another physical motivation for the analysis of Kekule distortions; in fact, it has 
been recently proposed, [5D1ET1[T2], that these kind of lattice deformations may be 
at the basis of a possible electron fractionalization in the low energy excitations of 
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graphene, such as the one taking place in one dimensional systems like polyacetylene, 
[801 [58]. 

The Chapter is organized in the following way: in Section [5. 21 we slightly change 
the definition of the model, in order to take into account the effect of small lattice 
distortions; in Section 15.3.11 we describe the integration of the ultraviolet degrees 
of freedom; in Section 15.3.21 we discuss the integration of the infrared degrees of 
freedom, and we prove our result (jl.4llip on the two point Schwinger function; 
finally, in Section [5.41 we discuss the issue of a spontaneous lattice distortion, and 
we derive a variational equation for the gap. 



5.2 The model in presence of Kekule distortion 

In this Section we change the definition of the model so that the presence of a small 
lattice distortions may be taken into account. The Hamiltonian of the model in 
presence of a fixed lattice distortion is given by: 

where "Ha is the Hamiltonian defined in (|1.21ip while 'H'j^^^ {{(f)}) describes the 
electron-phonon interaction and it is given by: 

XGA cr=t4- 

1=1,2,3 

where <j)g^i is a real valued function corresponding to the distortion of the bond 
which connect x to x + 5i. As for the case (p^^i = 0, this model can be studied using 
a functional integral representation; for any positive /3, L, the average of a generic 
physical observable 0{A,'^) is defined as 

:= ^ /"p(d^)P(dA)e^(*'^)+''^(^'*)0(A^) , 03) 

P(d^)P(d^)e^(*'^)+''^(^'*) , 

where all the objects appearing in (|5.2[3p have been defined in Chapter [TJ except 
for V(i,{A, ^) which is given by: 

v^{A,^):= ([5214) 

x£A 
1=1,2,3 

We shall perform a specific choice for (pg^f, namely, we shall choose (f)g^i corresponding 
to the so-called Kekule distortion of the honeycomb lattice. This is equivalent to 
the choice 

= 4y ■= cos{pp{6j - 6j, - f )) , ([S215) 

where Aq G K has to be understood as a "small" parameter. The Kekule distortion 
corresponds to a particular dimerization pattern on the hexagonal lattice; see Fig. 
[5Jl 
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5.3 Renormalization Group analysis 

In this Section we shall describe the RG analysis that will allow us to characterize 
the ground state of the model defined in Section 15.21 For simplicity, we start by 
discussing the evaluation of the free energy; as already seen in the previous Chapters 
for the case Aq = 0, the analysis of the correlation functions will be a straightforward 
extension of the procedure that we are going to discuss. 

5.3.1 The ultraviolet integration 

Proceeding in a way completely analogous to Section 12.41 it follows that after the 
integration of the ultraviolet degrees of freedom the partition function of the model 
can be rewritten as (from now on the wide tilde will recall the Ao-dependence): 

where Fq and V^^^ are given by (repeated indeces are summed): 

/n m 
i)P2i— lj'^2i — 1 k2j,(Ti,p2i,a;2i J. J. 
„ . ^ J 1 1 



n+m>l 5.,P i=l i=l 



I 2n 

\j=i i=i 



rOo) 

^ 



Fo := imr'W^T , 03) 

and the kernels W^^ln can be computed in terms of Feynman graphs following the 
rules described in Appendix [D] and Section 13.21 with the only difference that now 
new vertices proportional to Aq appear; these vertices have two fermionic lines and 
an arbitrary number > of photon lines. Clearly, if we set Aq = the kernels 
appearing in (|5.3l2p reduce exactly to the ones of Section I2.4| and in particular 
verify the properties stated Lemma I2.4.2| which are implied by the symmetries 
listed in Lemma [2.3. H for Aq ^ the role of Lemma [2.3. H is played by the following 
one. 
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Lemma 5.3.1 For (f)^i given by !i5.IA5\) . both the Gaussian integrations P{d^), 
P{dA) and the interactions V{^,A), v^{A, defined in ^KEM)> (E3I2P, ^EEM, 
respectively, are invariant under the transformations (1) - (8) of Lemma 
\2. 3.1\ provided that: 

• under (4) jo — ^ jo + 1; with 3 + 1 = 1; 

. under (6.b) jo Jq, with I' = 1,2' = 3, 3' = 2. 

Proof. The proof is a trivial extension of the one of Lemma I2.3.1| the details are 
left to the reader. ■ 

For Aq 7^ the above symmetry properties no longer exclude the presence of 
a mass in fermion propagator; however, as proven in Appendix |Hl the following 
identities are true: 

^^ifoLjO) = t^^aL,.(0) = , 
^^2AP,3-P,.,-J0) := Ao5o, <5oGK, 
^ifoi3-p,.,-JO) = e-'^^(^--^^^)Ao5o . m^) 

where 5o = J2n>2^o^^ ^ M with 6q^^ bounded proportionally to (A^/2)!, uni- 
formly in Aq. It is convenient to rewrite the partition function (|5.3[ip using the 
more compact relativistic notations introduced at the end of Section I3.2t it follows 
that: 



g-/3|A|Fo J p<,(d^(<0))p^^(^{<0))gV(0)(^(<0),A(<0)) ^ 

where the interaction V^^\ip,A) is given by: 

/c^k' f fik' 

— V^k',<xi7M2/^,0^k',a - J V^k',ai7^^°^Ao5o^k',a + 

f dp ^ ^ _ fdk' dp - , 

- J ^i^,p^i^-p^i^fi + J -^^^:^iAp,o'i/'k'+p,a7pV'k',a^p,p + 

where: 

(i) the first terms takes into account the contribution of c^p^ifp.LC^) I^^^^q- 

(ii) The second term is a mass term for the fermion field; the 4x4 matrix ^(j'o) is 
defined as: 

7^^"):= c J . (15317) 

(iii) The third term takes into account the contributions of VFq"'2^^(0) (which verify 
the same symmetry properties of their Aq = counterparts, see Appendix IH]) . 

(iv) The fourth term takes into account the contributions of 1^2"'i p ^(0) 

(v) The fifth term collects all the higher powers in the fields, and the Taylor 
remainder of the expansions of the kernels in the momenta and in Aq. 
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5.3.2 The infrared integration 

We proceed as in the Aq = case; we evaluate the partition function by integrating 
the fields in an iterative way, starting from 

^(o)^^(o)_ Yor h > h, where h = 
h{Ao) < is a suitable scale label that will be defined in the following, the multiscale 
integration will proceed in a way analogous to the case Aq = 0; after the scale h* 
the fermion propagator will become massive. 

Massless regime. We want to inductively prove that after the integration of the 
fields on scales > h + 1 > h the original functional integral can be rewritten as: 

-mh f p<,(dv.(^'^))p<,(dA(^^))e^*''(v^'^'-'"'^'-'') , 



=/3,L = 

where P<:h{dA^-^^) has the same propagator of the Aq = case, while the fermionic 
propagator is 

^ C^^-^^T^ ,2 , - n /^2l^/l2 , X n /^2 1 + ^h(k ) ; ^9) 

finally, the effective interaction on scale h can be written as 

/n m 

(m 2n \ 

^p,+^(-irkn , mm 
j=i i=i J 

and the kernels VFg^''^""'^ together with the free energy F^, will be defined recur- 
sively. In order to inductively prove formula (|5.318p . we split the kernels appearing 
into V'-'^^ as sums of local plus irrelevant parts; this procedure will be slightly differ- 
ent with respect to the one discussed in the case Aq = 0. Let us define the operators 
Ci, i = 0,1, in the following way. 



If n = 1, m = 0: 



if n = 0, 771 



if n = 2, ?Ti = 1: 



if n = 0, m 



A<3!f ^(pi, P2) := <3!f V, 0) ; mm 
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then, given a generic kernel 1^2^^^'''°^^ we define the operators Vi, i = 0, 1, in the 
fohowing way: 

^i^^in!i2M({ka,{p,}) := 



Finally, given £j, Vj^ j = 0,1 as above, we can rewrite the effective potential on 
scale h as CV^'^^ + TZV^^\ where the localization operator C is defined as follows. 

If n = 1, m = 0: 
if n = 0, m = 2: 
if n = 2, m = 1: 

^w'ffi? (P) := ^oT'oW^Sf^Cp) ; BIS) 



if n = 0, m = 3: 



'C<3!rHPi,P2) := £o7'o<3!rnpi,P2) . qio) 



The combinations of the operators Ci and Vq give rise to terms equal to those 
discussed in the case Aq = 0; therefore, it remains to discuss the combination of the 
operators Ci and Vi. Proceeding in a way completely analogous to the one followed 
in Appendix iHl to prove (|5.3i6p . the symmetry properties (1) - (8) of Lemma 15.3.11 
which are preserved by the multiscale integration, imply that: 

^k'.a-CoPiW^io^'^'^^kOVk',. = -iA/,5;,V^k',a7^'"Vk',a , 4 = SI 153120) 
^oT'i VFo^Sf (P) = -^m,M2^'*^Mi,/^ ' ^i,h = y2,h , h = i>: . 

As we are going to see, the remainders of the Taylor expansions in Afc(k') necessarily 
appearing in the definition of the TZ = 1 — C operation, see (15.31 16p - (15.31 19p . admit 
gains in their dimensional bounds with respect to the basic ones, and for this reason 
will give rise to irrelevant kernels. 

Hence, one can repeat exactly the same steps performed in Section[22]to describe 
the integration of the partition function of the Aq = case; the only difference being 
the presence of a new running coupling function A/i(k'), which evolves according to 
the following flow equation: 

Z;,_i(k')A^_i(kO = Z^k')A^k') + Z^AhShXh{^') ■ (IS3121) 

Setting A/j_i := A/i_i(0), it follows that: 

-(l + /3f), (1022) 



A. Z. 
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where = f3^'^^'^ + 0{e^M^^) + O(e^), and (see Appendix El : 
:;A,(2) _ / 2 _L 2 2 ^ arctan ^fe 1 + ^, 
3e2 



«A,(2j _ Cp2 , o 2 2 \fm:f:lii^lZJiilnp- /If 



- log M + 0(e^ log M) ; ([53123) 

therefore, the flow equation (|5.3122p and the explicit computation (j5.3123p together 
with the result ()3.4[19p for the Beta function of imply that 

A,, ~ AoM-'^^'' , r/^ = 14 + O(e') ■ (153124) 

37r^ 

As already mentioned at the beginning of this section, the multiscale integration 
goes on until h = h, where the scale h is defined as the scale on which the fermionic 
propagator becomes massive, that is when: 

A^^ = M~^^h^ logA, Ao , (153125) 

1 



which implies that A^ ~ ^o^'' ' • It remains to show that the actions of 1 — Pi 
on a relevant kernel and of 1 — on a marginal one produce irrelevant kernels; in 
particular, we shall show that the operator \ — Vi acting on a kernel on scale h > h 
is dimensionally equivalent to M^^~^^^'^~^^\ In fact, it is straightforward to see that 
the action oi 1 — Vi is dimensionally equivalent to (A/j/M'')*"'"^; and 

where we used the definition ()5.3[25p of h. 



Massive regime. At the scale h the fermion propagator has become massive; 
however, the photon propagator is singular in the infrared, therefore we still need 
the RG to analyze this regime. We could integrate the residual fermionic degrees of 
freedom in a single step and study the effective bosonic theory by multiscale analysis; 
however, in this way we would get estimates not uniform in the bare fermion mass. 
For this reason we will proceed in a way analougous to the massless regime, i.e. we 
will "slice" the fermionic and bosonic degrees of freedom simultaneously. 

The main difference with respect to the previous regime is that now we do not 
expand anymore the kernels in powers of the fermion mass; in fact, the splitting 
performed in the massless regime would be not convenient now, because, as we are 
going to see, the momentum flowing on a scale < h will be smaller than the corre- 
sponding fermion mass (and therefore A/1/7'^ will be greater than 1, see (j5.3i26p ). 
As before, we will localize the kernels corresponding to monomials of order < 3 in 
the fields; the main difference in the single scale integration with respect to the 
previous regime is that here we do not rescale the fermionic fields. In fact, because 
of the improved dimensional bound on the fermion propagator the flows of the run- 
ning coupling constants corresponding to monomials k^ipip converge to finite values. 
The same is true for the mass term A/,, which remain close to its value on scale 
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h, and for the charges e^^h- Finally, the photon mass is controlled using the same 
Ward identity exploited in the case Aq = 0, and the local parts of the remaining 
kernels vanish by symmetry. 

After the integration of the first \h\ infrared scales the partition function can be 
rewritten as 

H^^i = e-/^l^l^^ J P<^(#(^^))P<^(dA(^^))e^*''(^^'^'-'''^'-'*) , (153127) 

where the effective potential V^'*^ is given by (|5.311Up with h replaced by h. We split 
the effective potential as JOV'^'^^ + TZV^^\ where TZ = 1 — jO and the action of C on 
the kernels of V^'*^ is defined as follows: 

By the symmetry properties of Lemma l5.3.1l it follows that: 

\h = ^2,1 > h,h e ^ > mm 
dMl, = ' ^0,3,?'^ (0) = , mm 

where: (i) dOT^ is proved in Appendix EU (ii) (I5.3[3U|) . (I5.3[3l|) are proved 
in Appendix IH.2i the proof of the second identity in the third line of (I5.3132P is 
completely analogous to the one of the case Aq = 0, which has been discussed in 
Appendix IE. 41 (the only symmetry transformations involved are (7) and (8), and 
these do not involve the label jo); the proof of the first one is less trivial, and it 
goes as follows. 

First, consider the case j'o = 1; in this case the proof is completely analogous 
to the one reported in Appendix IE. 41 for the Aq = case, since the symmetry 
transformations involved do not change the label jo = 1- Let us now consider the 
case jo = 2; setting Too := 1, Tb^j = Ti^ := 0, Tij := Tij where T is the 27r/3 
rotation matrix, by symmetry (4) it follows that: 

f[{fd^)w^f^^\Q)\ T-i = d^w^f^^~^ = , mm 

which gives 

9oW(gf(0) = 0, 
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and the last two equalities imply that diWQ^2,'li^\^) — 0) i = 1,2. The same 
argument can be repeated for jo = S, and this concludes the proof of (|5.3[29p . 
Setting 

'Hih ■■= ^/^^Sf ^(0) , e^f, := W^:^^{0, 0) , ([53135) 
the partition function can be rewritten as: 

= / P<^(dV'(^'^))P<;.(dA(^'^))e^^'''(V^^'^'''^'^'') , (I0136) 

where: (i) the propagator of P^j^^dil)^-^^) is obtained by adding to the exponential 
of the gaussian weight of P^f^{dilj^-^^) the mass term of (I5.3l29p . i.e. it is given bj0 

^(^^)(kO := X^i7o^o + i^,(kOfe-- .-i.Oo>A,_,(kO 

■■= ^/.(k')A^(k') + Z^A^A^ ; ([0137) 
(ii) the interaction V^^\^l^,A) is 

VW(^,A) :=/ [e^,;,ji|^U,,p - M^P^^^i,,pi,,_p] + 

+7^V(''H^, ^) • ([01 38) 

Let us write, as usual, ^^-''^ = ^^'^^ + ^(''), A^^^'*) = A^'*) + where yl^'^), V^''^ 
have propagators given respectively by (|3.2i26|) with h replaced by h and 



)W(k) 4(k)i7oA;o + i^/,(k')^-7-i7^''''^*A/^_,(k') 



A:2 + z}^(k')2|fc|2 + A^_,(kO 



l + ^/.-i(k')) 



4(k) = S39) 

integrating the fields on scale h along the lines of Section [3.31 we recover our starting 
formula (|5.3l27p with h replaced hy h — 1. The procedure can be iterated, and to 
prove the infrared stability of the renormalized expansion we have to check the 
boundedness of the flow of the running coupling constants. First of all, the flow 
of Dfj^^h can be controlled using again the Ward identity ()3. 51240 . Regarding the 
remaining running coupling constants, assume inductively that A^ ~ A^ for all k > 

h; the flows of A^, z^J°Jf^, e^^^^^ can be controlled exploiting the fact that the fermion 

propagator on a given scale k > h is bounded by ~ M^, that is its dimensional 
bound admits a gain M^~^ with respect to the corresponding dimensional bound 
in the massless case. In fact, consider the Feynman graph expansions of the Beta 
functions of A/j, z^J°Jf^, ^^"'a'^hi there are two possibilities: either there is at least one 



^Notice that now A,~_;^ ~ '^h-iW = + /^f ) with ~ 5^; compare with (|5.3l22p 
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(massive) fermion propagator on scale k such that h < k < h, or the dimensional 
bound of the graph is depressed by a factor M^'^^~^\ because of the short memory 
property. Therefore, in both cases the bounds admit an improvement of a factor 
^e{h-~h) respect to the usual dimensional bounds, and for this reason the Beta 
functions of the remaining running coupling constants are summahle. This means 
that for h < h 

This concludes the proof of the boundedness of the flow, and our discussion on the 
multiscale analysis for the free energy of the model in presence of Kekule distortion. 



5.3.3 The two point Schwinger function 

We conclude this Section by discussing the evaluation of the two point Schwinger 
function; in particular, we will prove our main result (jl.4llip . To procedure is 
completely analogous to the one followed in the case Aq = 0; therefore, we shall 
only discuss the differences with respect to the case Aq = 0, without repeating the 
whole argument. We integrate the fields scale after scale, and the main difference 
with respect to the case Aq = is that formula (|4.2l4p is replaced by 

^ v('')(v^^,A)G('^+i)(kO</.k',. , (153141) 



where the matrix G^''-*(k') is defined as 

G(k') := ^ ^^QW(k') , mm 

i=h ^^-1 

while the matrix Q^^'^ila') is given by (repeated indeces are summed) 

QW(k') = Q('^+i)(k') - \Z^[z^^^^,k'^ + A,5,7(^°^] G('^+i)(kO {h>h*), 
Qi?a.(kO = Qtl'K^) - ^Zh^ H'^t-^ + ^^^^^^^^A] G&t'^k') [h < h*) 

The integration of the single scale and the check of the inductive assumption is 
done following the same strategy discussed in the case Aq = 0, taking into account 
the modifications of the multiscale analysis introduced above in the evaluation of 
the free energy for the case Aq 7^ 0; in particular, after the fermion propagator 
has become massive we stop to rescale the fields, and we collect all the new terms 
proportional to k^ipij}, A^ipip into new running coupling constants z^J^Jf^, ^^i^°a h "^hich 
behave dimensionally as marginal ones, but whose flows are bounded because of the 
improved dimensional bound on the fermion propagator. 

The remaining part of the discussion is the same as for the case Aq = 0, and 
we refer the reader to Section 14. 2t the two point Schwinger function is obtained 
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deriving the generating functional with respect to the external fields (p^^a: 'P'ka- 
particular, notice that the leading contribution to the two point function is given by 
the dressed propagator (j5.319|) or (|5.3137p , depending on whether the external quasi- 
momentum k' is such that |k'| > A, |k'| < A; the former propagator is essentially 
equal to the one of the case Aq = 0, while the latter is different because the running 
coupling constants appearing there are equal to their values reached on the mass 
scale h* . This concludes the proof of our result (jl.4ll0|) . 



5.4 Gap generation 

So far, we have seen that the Kekule distortion of the honeycomb lattice is strongly 
renormalized by the electromagnetic interaction; in particular, the ratio of the 
dressed and bare "Kekule masses" A/Aq blows up as Aq — )■ 0. The natural question 
is whether it is possible that, under suitable conditions, a spontaneous distortion 
(corresponding to a gap generation in the fermionic energy spectrum) may emerge. 
To understand this, we shall consider the hopping parameter as a dynamical vari- 
able depending on the lattice site and on the bond, and the resulting model will be 
studied in the Born - Oppenheimer approximation. 

5.4.1 Gap equation 

The Hamiltonian of the model in presence of a dynamical lattice distortion is 



^a({0}) := ^a({0})+/Ca(W) , 

^a({0}) ■■= J2 ^h^ dsan) 

1=1,2,3 



where /Ca({</'}) is the elastic energy of the distortion iji G M, k G M"*" is the stiffness 
constant and g' S M is the (classical) phonon coupling. The free energy of the model 
in the Born - Oppenheimer approximation is given by: 



e-/5|A|^J'£ := [[ Y[ d,/),.^,]e-/5'CA(m)-/3|A|F,,,({0}) 
xeA 

i=l,2,3 

n 



/,..]e-/3|A|S/3.L({0}) 



xsA 
i=l,2,3 



Fp,Lm) := --^logTYe-^^^««) ; 

the total "energy" Ei3^£^{{(j)}) is extremized by values of the distortion (pt^ satisfying 
the condition dcf,^ -Ei3^L{{(j)}) = 0, which is equivalent to: 
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Now, if the equation (|5.4I3|) admits a non-trivial solution, and if this solution corre- 
sponds to a global minimum of the energy Ejs l{{<P}) uniformly in /3, |A|, then from 

liin ^ Ep^Lim) , dsns) 

|A|— !>+oo |A|— S-+00 

that is the integral over all the possible distortions in (|5.4i2p is dominated by a 
single configuration (f>* , which therefore corresponds to the lattice distortion in the 
ground state of the system. 

Clearly, the problem of solving (I5.413P is very hard, and we will not be able to 
solve it in full generality. However, we will show that the equation (|5.4l3p admits as 
a solution a configuration (j)* which corresponds to a particular Kekule distortion; 
and at least for small g this solution corresponds to a local minimum of 
This last fact follows because 

-/Ca({(?!)*}) = —Sg^^6ij , 



'x,iOCPy,j 5 

F^Am) = {pff-^pffy,^^^ (15316) 



and for max^^j 14''^ i\ small enough the — correlations are bounded objects. 
The problem of showing the uniqueness of the solution and its correspondence to a 
global minimum of the total energy is beyond the purposes of this thesis. 

5.4.2 A solution to the gap equation 

As a first step, it is convenient to rewrite the r.h.s. of (j5.4l3p in the functional integral 
representation. The observable p^^^'* is manifestly gauge invariant; therefore, we 
can rewrite its statistical average as follows: 

(P^f^K'r = 5$(„.).log/p(d^)P(d^)e^(^'^)+^V*(*.^)+(*.P^"+') (Eai7) 

a,i 

where repeated indeces are summed, and as usual J = dxo J2xeA- 

We perform the following special choice for (/)^^, corresponding to a Kekule 
distortion of the honeycomb lattice: 

= 4;?°^ ■■= <^o + cos{pp{5, - 4 - x)) , ms) 

where the parameters (pQ, Aq will be determined starting from (|5.4l3p : as we are 
going to see (po will correspond to a renormalization of the bare Fermi velocity, while, 
as it follows from the discussion in Section \5.2\ Aq, the amplitude of the Kekule 
distortion, behaves as a bare mass for the fermion propagator. Let us denote by 
(• • • )^'^ the statistical average in presence of a fixed distortion and call P^^^"^ 

( F ) 

the Fourier transform of p)^^ ■, see below; the following result is true. 



5.4 Gap generation 
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Lemma 5.4.1 For any j3, L, choosing of the form ^5.4 ~8) with (j)Q G 
the gap equation ^5.4\3^ is equivalent to the following two equations: 



An G 



Proof. It follows thalU 



90 

Ao 
3 



K |A| ^^0,1 'p^L 



a" 1 



K |A| '^Pj'l 



/3,L 



(Elio) 



/3,L 



|A| A ^^fc.* 



1 e-iPF-^'(p; 



(153112) 

Oo) 



|A| 



where the last equality is implied by the fact that (p^ • 



Oo) 



reflects the periodicity 



of the distorted lattice; see Fig. 15.11 Therefore, from (|5.4[8p and (j5.4[12p we get 



Ao 
3 



K |A| ^^Oj I fi^L 



9 



1 . 



K |A|^^PfJ '/3,L 



3113) 
3114) 



from (|5.4[13p we see that a necessary condition for (j5.4[8|) to be a solution of the gap 
equation is that the r.h.s. of (|5.4ll3p and (j5.4[14|) are independent of j, jo and w, j, joi 
respectively. To prove this, we shall use that the generating functional appearing in 
(|5.4[7|) is invariant under the symmetry transformations of Lemma 15. 3. H provided 
the external field appearing in (|5.4I7|) transforms as discussed in Section 14. 3[ 

Let us first consider the first equation in (|5.4[13p . The invariance under (4) and 
(6.a) imply respectively that 



,.(i?+)\{io-i) 



(A 



/3,L 
(io) 
0,2 I Rj, 



,.{E+)Aio) ^ , (£+), (io+i) 

\^0,2 'i3^L ^^0,3 

(io) 
0,3 'fi^L ' 



15) 
16) 



■^To prove (|5.4I12|) . note that aU the sites equivalent to x, that is those with the same configura- 
tion of distorted nearest neighbours of x, can be obtained starting from x and moving of 2ai — 02, 

X + 2ai — a2, X — s> x + 202 — 0,1. This 



2a2 — Si; tlierefore, (p^^^')^"'')^' must be invariant under x - 



imphes that in only the Fourier modes verifying the following constraints can appear: 

13, L 



-fci + — 



2nni 



3%/3 



k2 



2im2 , 



m, 712 



gZ 



Elii) 



from (|5.4llip we get fci = ^{ni +712), = ~ 1^2), and the fact that k £ Vl implies < 

ni +712 < 1, —2 < 711—712 < 2. Therefore, the allowed integers are n = (1, 0), (0, 1), (—1, 1), (1, — 1): 
the first two choices correspond to p^, p^; the last two correspond to — ~ Ppi Pf ~ Pf ~ Pf 
(these identities have to be understood modulo vectors of the reciprocal lattice). This concludes 
the proof of (|5.4I12|) . 
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therefore iP^(^j^^)^^^ independent of j and jo; moreover, from symmetry (5) we 



get that 



(153117) 



Consider now the second hne of ()5.4[14p . First of all, let us prove that the r.h.s. is 
independent of j,jo; to see this, notice that because of symmetry (4) the problem 
is equivalent to show that 



iP%A)^\ ^^^^'i^Fi^n-j+i-^^) is independent of j,jo 



^ 18) 



where we used that f^Si = (T~^j^)5i = p^(5j+i, with T the 27r/3 rotation matrix. 
The claim (|5.4118p can be checked explicitly using again symmetries (4) and (6. a); 
in fact the combination (4) - (6. a) - (4) gives: 



" Pf'-^ ' fl r, " Pp'^ I3,L 



^ Pf^^'I3,L Pf''^'/3,L Pp''^ fi,L Pp'^ (3,1 



PP'^'(3,L 



where in the last equality we used again the invar iance of under T. Therefore, 
the above chain of identities proves that 

finally, from (6.b) and (5) we get: 



(p'^V) 



F' 13, L 



{P 



Pp 



(153120) 



and this concludes the proof of Lemma 15.4.11 



To conclude, we have to show that the equations (15.41 13p . ()5.4[14p admit a non- 
trivial solution; in particular, we will be interested in understanding the effect on 
the electromagnetic interaction on the solution. First of all, we rewrite the r.h.s. of 
(|5.4113p . (|5.4114p as (we set (•••) = (•• • )^^^^ for notational simphcity): 



+ 



53121) 



<T=n 



x+5i 



+ 



3122) 



■E 



(at b- . fe^^A^.i) - l)) + (6+ . aZ fe-'^^A^^i) - 1 
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the terms appearing in (j5.4l2ip . (j5.4l22p involving the exponential of the photon 
field can be studied again using multiscale analysis and RG, and it follows that they 
give rise to O(e^) corrections (we omit the details); therefore, using the functional 
integral representation in the Feynman gauge, we can rewrite (I5.4l2ip and ()5.4l22p 
as: 

mCi""^^ = -^l%[{Kt,2Kt,i) + (KuK^,2)]+0{e'), (153123) 



|A| 

= -2/— [(^^ ^ ^+ \ + (^- ^ ^+ 
|AP^4.i^ J D \y k+pj,t,2 k+p^,t,i/^\ k+p+,t,l k+p-,t,: 

+ 0(6^) 



+ 



Notice that in the integrals appearing in (j5.4[23p the main contribution come from 
A < |k|; in fact, for |k| < A the two point function is bounded by (const. )A~^, and 
therefore the integral over |k| < A is O(A^). The gap equations (|5.4113p . (|5.4il4p . 
together with (|5.4l23p and our result on the two point function in presence of dis- 
tortion, can be in principle studied numerically, and quantitative predictions can 
be made; here however we only want to get a qualitative picture of the behavior of 
the solution as a function on the interaction. For this purpose, we can approximate 
the equations for (/iq, Aq as follows: 

^ ^ ^ ' D Z(k' + p-) A;2 + 7;(k' + p^)2|k'|2 

A<|k'|<l 

^ ^ Vy- / — - ^^'^^ + (E3124) 

3 K ^ J D Z(k' + p^) A:2 + 7;(k' + p-)2|k'|2 ^ ' 

A<|k'|<l 

where 

A(k'+p^) = Aolk'r^^^ , Zik'+p'^p) = |k'|-^ , v{k') = l-{l-v)\k'\'^ ; (153125) 

the gap equation (I5.4l24p is similar to the one found first in [62], see also [16], in 
the context of Luttinger superconductors. From (j5.4[24p we get that (po = 0{g^); 
the discussion for Aq is less trivial, and it goes as follows. 

"Small" electron-electron interactions. Let us consider the case of "small" 
electron-electron interactions. Notice that in order that the effective Fermi velocity 
flows from v, the bare value, to 2v the momentum p has to be equal to: 

(1 — \ ^ 
)"=^p~e~W for 77 < v; (15.41 26) 
1 — t; / 

therefore, for fj <^ v the Fermi velocity is for all practical purposes equal to v, and 
the integral in (j5.4[24p is equivalent to 



f-\ , Aop''-''^^ arctan e 2Ao (1 - A^-^^+i 
dp dcosO — 27^^ 2"^ = 2 TT" 

where ^ = V^"^ — 1- Plugging this result into ()5.4[24p we get that 



^o^{l- ^) , for \g\ >go = 0{^v) , ([53128) 
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and in particular the critical coupling go > is essentially independent of e. How- 
ever, still the electromagnetic interactions play a non-trivial role. In fact, since 
T] — 7]K = + O(e^), the exponent appearing in (|5.4[28|) is greater than 

1; this means that close to go (I5.4[28|) behaves as Ao(5') ~ (5^ — g'o)^'^'^ with 
a = yII^- -|- O(e^) > 0, that is 9gAo(g) is continuous at \g\ = go: the transition 
is smoothened by the interaction. 

"Strong" electron-electron interactions. The discussion in this case will be 
more speculative. From the analyses of Sections 14.31 and 15. 3| one may reasonably 
guess that electron-electron interactions tend to favor excitonic instabilities, and 
in particular a Kekule distortion of the honeycomb lattice: as we have seen in 
Section 14. 3| the excitonic response function is amplified by the interaction, 
while, as discussed in Section [5?3l a preexisting Kekule distortion is greatly enhanced. 
This belief is corroborated by the fact that the integrand appearing in the second 
equation in (I5.4124P scales as Ikl^"**^"^, with rj — r]K = + O(e^); therefore, 

if rj — r]K — ?• — 1, a situation which corresponds to a strong coupling regime, the 
integral in (j5.4l24|) is bounded from below by a quantity of order |logA|. This 
means that in this strong coupling limit the gap equation (j5.4l24|) admits a solution 
for any \g\ > 0; in particular, if r]K — r] — 1 exceeds then 

Ao (1 + ^^^^~^^~^^ ) , for \g\ > 0. ([53129) 

Hence, it is reasonable to expect the critical value of the phonon coupling go to be 
generically lowered by the electromagnetic interaction. 



Chapter 6 

Conclusions and perspectives 



In this Thesis we investigated the effect of electromagnetic interaction on electrons 
on the honeycomb lattice at half- filling (1 electron per site in average), with the 
physical motivation of understanding the role of electron-electron interactions in 
the low energy physics of graphene. In particular, we introduced a new model for 
graphene, where the the coupling with a quantized electromagnetic field is defined 
so to respect gauge invariance; therefore, our model can be seen as a "lattice gauge 
theory model" for graphene. Compared to the existing literature, the model intro- 
duced here does not neglect the lattice (no Dirac approximation is involved) and 
does not neglect the presence of the vector potential, which is necessary in order to 
guarantee that the full theory is invariant under local gauge transformations. 

Using rigorous Renormalization Group methods and lattice Ward identities we 
succeeded in characterizing, at all orders in renormalized perturbation theory, many 
of its ground state and low energy properties, [39]; in particular, we investigated: (i) 
the effect of the interaction on the two point Schwinger function (the dressed prop- 
agator); (ii) the effect of the interaction on some response functions; (iii) the effect 
of special lattice distortions, the Kekule ones; (iv) a mechanism for spontaneous 
lattice distortion. 

Regarding the two point Schwinger function, its scaling properties are deeply 
modified by the presense of the interaction; in fact, interaction-dependent anoma- 
lous exponents appear. This means that the interacting system behaves as a Lut- 
tinger liquid: as far as we know, this is the first time that Luttinger liquid behavior 
is found in a two dimensional system. Moreover, there is an emergent Lorentz sym- 
metry, due to the fact that the effective Fermi velocity (that is, the one appearing in 
the two point Schwinger function) tends to the speed of light at the Fermi surface. 
After this, we studied the effect of the electromagnetic interaction on the excitonic, 
charge density wave and density-density response functions (or susceptibilites); we 
have found that, as for the two point Schwinger function, these response functions 
scale with interaction dependent anomalous exponents. In particular, some response 
functions are amplified by the interaction, namely those corresponding to Kekule or 
charge density wave pairings; the former is associated to a Kekule distortion of the 
honeycomb lattice, see Fig. 15.11 the second to a periodic alternation of excess/deficit 
of electrons on the A/B sites of the lattice. 

To get a deeper understanding of the effect of lattice distortions we changed the 
definition of the model in order to take into account a preexisting Kekule distor- 
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tion; interestingly, the amplitude of the bonds deformation behaves as a bare mass 
for the fermion propagator, which is strongly renormalized by the electromagnetic 
interaction. This is interpreted by saying that Kekule distortions are enhanced by 
the electron-electron interaction. Finally, motivated by the results of the response 
functions and on the mass renormalization, we investigated the possibility of spon- 
taneous lattice distortions. We treated the lattice distortion as a classical dynamical 
variable (a classical phonon field), in the Born-Oppenheimer approximation; by us- 
ing a variational argument we have found that the total energy is extremized in 
correspondence of a suitable Kekule distortion, whose amplitude is determined by 
a well-defined self-consistence equation. We studied this gap equation from a quali- 
tative viewpoint; in general, one sees that the solution of the equation as a function 
of the phonon coupling g becomes non-trivial for g greater than some critical value 
go- For weak electromagnetic interactions we have found that the transition from a 
gapless to a gapped phase is "smoothened" (the interaction removes the discontinu- 
ity at go of the first derivative of the gap with respect to g) , while go is essentially 
interaction-independent; however, from the same equation we got evidence that 
strong electromagnetic interactions favor the emergence of a Kekule instability by 
lowering the value of 50- 

Regarding the perspectives of our work, these are many both from a physical 
and a mathematical point of view. For instance, from the physical point of view, 
it would be very interesting to understand the issue of universality of the electric 
conductivity in graphene; in fact, it has been experimentally observed with very 
high precision that the optical conductivity is equal to Ae^ /h: remarkably, this 
value is not affected by higher order corrections and material parameters, such 
as the Fermi velocity. The explaination of this phenomenon is very challenging 
from a theoretical point of view, and some attempts have already been made in the 
literature; however, all the existing arguments are strongly based on the logarithmic 
divergence of the Fermi velocity, which is, in our opinion, an unphysical feature of 
the models considered so far. To show that the growth of the Fermi velocity does 
not play any role in this remarkable cancellation it would be very interesting to see 
first if it can be recovered in a case in which the effective Fermi velocity is weakly 
renormalized; for instance, one can consider the two dimensional Hubbard model 
on the honeycomb lattice, whose ground state has been rigorously constructed in 
[36]. In that case, the effective Fermi velocity turns out to be an analytic function 
of the coupling, close to the bare one. 

Another interesting physical problem is related to the insensitivity of graphene 
to localization effects usually induced by disorder; to understand this issue, it would 
be very interesting to add some randomness to the model defined here, and see if our 
RG analysis can be adapted to take into account the disorder. The effect of disorder 
on Dirac fermions has been investigated well before the experimental discovery of 
graphene. For instance, in a pioneering work, [251 [26] , Fradkin considered a model 
of noninter acting Dirac fermions in presence of disorder; by a replica analysis it 
was shown that 2 is the lower critical dimension for a localization/delocalization 
transition. In particular, considering different species of fermions, in the limit 
A —7- -|- 00 it was shown that in two dimensions the mean free path of the electrons 
is always finite in presence of disorder, and that it grows exponentially as the in- 
tensity of the disorder goes to zero. Moreover, in the same formalism a universal 
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nonvanishing minimal value of the conductivity was predicted. It would be very 
interesting to see if a similar analysis can be reproduced using our rigorous RG 
framework in the case of finite A^, and eventually taking into account the effect of 
electron-electron interactions. 

From the mathematical physics point of view, the model that we introduced 
here opens the way to many interesting problems; the most interesting one is the 
proof of convergence of the perturbative series. In fact, our present results are 
perturbative, in the sense that all the quantities that we compute are given by 
series in the running coupling constants with finite coefficients, where the A^-th 
order contribution is bounded proportionally to {N/2)\; the bound is a consequence 
of the factorial growth of the Feynman graphs at a given order, and it is clearly not 
enough to prove absolute convergence. In purely fermionic theories the solution of 
this combinatorial problem is well-known, |60 | I30 | [7]. and it is based on the so-called 
Gram hounds; basically, one exploits the —1 arising in the anticommutation of the 
fermionic fields to show that the A^-th order contribution to the perturbative series 
reconstructs the determinant of an x matrix, which is bounded proportionally 
to (const. )^. Here the problem is that we have both fermions and bosons, and 
the Gram trick applies only to the fermions; to control the proliferation of graphs 
obtained exchanging bosonic lines we need some other argument. We expect this to 
be a hard but doable task; for instance, one could try to adapt the cluster expansion 
ideas of [SU [32], or perhaps the recent "loop vertex expansion" of [HI [72]. If one 
succeeds in this, it would be the first rigorous construction of a two-dimensional 
Luttinger liquid. 



Appendix A 

Grassmann integration 



In this subsection we recall some basic concepts of the Grassmann calculus. 

Let us consider a finite dimensional Grassmann algebra, which is a set of anti- 
commuting Grassmann variables {'^'^ ^~}a£A, for some finite set A. This means 
that: 

:=M>^„vI>^„', + vI/^;,vI/^„ = 0, ya,a'GA, Ve,e' = ±; (01) 

in particular ^'^^'^ = Va G A and Ve = ±. 

Let us introduce another set of Grassmann variables {d^'^ , d'^~}a£A, anticom- 
muting with and a liner operation, the Grassmann integration, defined 

by 

' = 1, f dK = 0, aeA, e = ±. ^2) 



If is a polynomial in a £ A, the operation 

J n d^td^-F{^) , dAlS) 



is simply defined by iteratively applying ()A[2p and taking into account the anticom- 
mutation rules ()A[ip . It is easy to check that for all a € A and C G C 

j a^>^a^> e ^ 

in fact e-*i^*" = 1 - ^'+C^'- and by (1112]) 

J d^+d^-e-*i^*- = C , das) 

while 

More generally, if if is an |A| x |^| invertible complex matrix, the above formulas 
are generalized in the following way: 
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Again, ()A[7p can be easily verified by using ()A[2p and the anticommutation rules 
(lAlip : in particular it follows that 

/ Yl d^+d^^e"^-^6^*^*^'^*^ = detM , ^8) 

while 

/ n rf^^ti^ae-S»..sA*^^^».'f7,i;-^+ = M'^,p, , (|A19) 
where M'^/j^, is the minor complementary to the entry Ma' pi of K. 



Appendix B 

Naive perturbation theory 



In this Appendix we will derive a "naive" perturbative series in the electric charge 
e for the grand-canonical averages computed with the Hamiltonian Ha defined in 
(jl.2il|) : in particular, we shall show that the Hamiltonian model introduced in Sec- 
tion 11.21 and the quantum field theory model defined in Section 12.21 share the same 
perturbative expansion. In this sense, we say that the quantum field theory model 
introduced in Section [2. 2l is the functional integral representation of the Hamiltonian 
model introduced in Section [1.2[ 

The Appendix is organized in the following way: in Section IB. II we derive a 
perturbative expansion for the grand-canonical averages computed with T-L\; in 
Sections IB. 1.11 IB. 1.21 we compute the fermion and boson propagators, respectively; 
finally, in Section IB. 21 we show how to represent the perturbative series introduced 
in Section [B. II in terms of Feynman graphs, and we show that the same perturbative 
series is generated by the quantum field theory model introduced in Section 12.21 



B.l Trotter product formula 

In this Section we will derive a perturbative expansion for the grand-canonical 
averages of physical observables computed with the Hamiltonian Ha; for pedagogical 
purposes, we start by discussing the case of the partition function, namely: 

The averages of physical observables can be studied in a completely analogous way. 
It is convenient to define the unperturbed imaginary time evolutions of the operators 
(^Q), with i = 1, 2, and of , as follows (x = (xo,x)): 

a± := e^o.'^^"ai q-'^o.a^^ , 6^ . := e^O'^^o^i ^ g-Wo,A^o 

^.o" ' x+(0,<5i),cr x+5i,a ' ^ ' 



{xo,x+s5j) ' 
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moreover, we define the unperturbed imaginary time evolution of V\ as: 

Va{xo) := e^'>'^''°VAe-^°-^''° (PS) 



(ie) 



i=l,2,3 



x+(o,5,).<^ ml 



Am 



+ 6+ 



(-ie)' 



x+(o,<5,)"'^''^ ml 



+ 



+y E (?^x-l)¥'(x-y,0)4o,j/o(^y-l) > 

where nx is the time evolution of the density operator. Our starting point will be 
the Trotter product formula, which consists in the following identity: 



lim 

n— >+oo 



V n / 



(IB14) 



in fact, using (|B[4p it is straightforward to see that the partition function ()B[ip can 
be rewritten as follows: 



Tr{e-'^^o>A} 



1 + 



1 



Tr{e-/^^o.A} 



5^ ±1 TY{e-^^"^AVA(ti)FA(t2) • • • VA(t„)} , 

n>l 



05) 

where the integral is over all the ti variables with 1 < i < n, under the constraint 
that they decrease in their index i, and the sign it is + if the number of V is even 
and — otherwise. 

Being T^cA quadratic in the fermionic and bosonic creation/annihilation opera- 
tors, the Wick rule holds for evaluating the ratios in the r.h.s. of ()B[5P : in particular, 
setting 



± 

X,(T,1 



it is possible to express the various terms in (IBISP a suitable integrals of products 
of expressions like, see Section IB. 21 



Wi,i'{^ - y) 



±LYC ^X,CT,p^y,(T',r'/ 

TV{e-/^«o.A} 



(IB17) 

(Ps) 
dElo) 



where xq — yo > 0, and as we are going to see in subsection IB. 1.21 Wij{x — y) 
Wij{y — x); (7+ and can be combined in a a single function: 



£/a,a'(x) : = 



(7+,,(x) ifxo>0 



(Pio) 



The functions (jBiOh . (jBilOh are called the bosonic and fermionic propagators, re- 
spectively; they can be explicitly evaluated, and the computations are reported in 
Sections IBXTI IBX21 



B.l Trotter product formula 
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B.1.1 Fermion propagator 

In this subsection we shall compute explicitly the fermionic propagator (|B[10p . The 
free fermionic Hamiltonian in momentum space is given by: 

the Hamiltonian can be diagonalized by introducing the fermionic operators 



with 7^ = lb, in terms of which we can rewrite: 

= ^ E E < - + \''^^^\HA^ ■ ®i^) 

Now, for X G A we define 

moreover, we set a^^^ := e^oA^'o a^^^e'^oA^'o and = e'^*A^±^e-^o,A^o_ ^ 

straightforward computation shows that, if — /3 < — yo ^ /3, and 1( "argument") 
is equal to 1 if "argument" is true otherwise is equal to 0: 

(T{a-,a+^,}) = (Pis) 



(5rr rr' 



^{xo~yo)\n{k)\ ^{xo-yo+l3)\n{k)\ 



(T{/3-,/3+ = (Pie) 



5. 



g-{xo-j/o)|Q{fc)| g_(a;o-j/o+/3)|n(fc)| 



and (T{a-,,/3+^,})^^^ = (T{/3-,a+^,})^^^ = 0. A priori Eq. mM and Eq. jBlTS]) 
are defined only for — /3 < — yo !^ /3, but we can extend them periodically over the 
whole real axis; the periodic extension of the propagator is continuous on the time 
variable for xq — 2/o / and it has jump discontinuities at the points XQ — yo E /3Z. 
Note that at the points xq — i/q = fSn the difference between the right and left limits 
is equal to (— which means that the propagator is discontinuous only at 
X — y = f37j X 0. For x — y ^ /3Z x we can write 

(T{/3-X<x'})«r = ItT E e-^'^^^-y) i ^. (Pl8) 



108 



B. Naive perturbation theory 



In fact, consider Eq. (|B[17P : if xq — ^ 0, by Cauchy Theorem it follows that 
^ -iko - \n{k)\ ~ IcJ^ -\z - \^{k)\ e-i-/3«gM^o-w) + 1 ' ^^'^^ 

where the integral on the r.h.s. is taken counterclockwise along the path C = C+UC- , 
with C-t := 1^ G C : = ib|r2(A;)|/2|. Assume that — yo > 0; the integral over 
C_ is vanishing (the integration path can be shifted to "^z = —oo, and the integrand 
is exponentially vanishing) while the value of the integral over C+ is determined by 
the residue of the pole in z = i|0(A;)|. It follows that: 

„-iz{xo-yo) I \n{k)\{xo-yo) 

(11120) 



c+ (27r) -iz _ |0(fc)| e-i-/3 + 1 1 + eWk)\ ' 

which is precisely the first term appearing in the r.h.s. of Eq. (IBIlSp . The same 
trick can be used to evaluate the other sums in ko appearing in Eqq. (jB[15[) , (IBI16P , 
and flBTTT]) . (piSl) follow. If xq - yo = /3n and f - y / 6, the r.h.s. of pTTTll is 
equal to 



=: (T{a,,,a+,,}) ; (021) 

J' p,i^ xo-yo=l3n 

the same holds for (Tj/J^ ^,})^ ^. If we now reexpress a^.o- ^-iid /3x,ct terms of 
«x,a and b^^^^g^-^y using (|B[12p we get: 

5o(x-y):=(T{^',,,^+,,}) = E e-*(--'')5(k) (022) 

^ ^ ^-ik(x-y^/ -n*{k)\ 

It is important to stress that (jB[22p agrees with dHHSI), flBTlGll as long as x - y 7^ 
/3Z X 0; in fact, if x = y we have 

while from (IbTTs]) . (IbTIg]) we see that 

5o(0-, 0)1,1 = 5o(0-, 0)2,2 = + (IH124) 

1 . e/3|f^{fc)l e-'^l^^^)! \ 1 



2\A\ Vi_^g/3|c(fc)| i + e-/5|f^WI^ 2 



B.l Trotter product formula 
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B.1.2 Photon propagator 



In this subsection we shall compute explicitly the photon propagator (|BI8p . The 
free bosonic Hamiltonian is: 



LAa 

and the quantized vector potential is: 
1 



P ' 



(11125) 



r=l,2 



LAa 



r=l,2 



X[h'fl]{\p\) 



PC + E C 



(11126) 



Let 



^ijiP) ■= E (%r)i(4,r)j 



r=l,2 



r _ PiPj 



hj = 1,2,3 



(IB127) 



where the last equality can be read as the completeness relation of the orthonormal 
basis |e;p.i, Spi, £p^2i £p,2' °1 ® C; a straightforward computation shows 

that: 



(11128) 



-\xo-yo\ 



:,{\xo-yo\-P)\p\ 



1 - e-^ll^l 



+ 



1 - e 



To derive ()B[28P it is useful to recall that, because of the commutation relation 
(|1.2[4|) defining the creation/annihilation operators: 



^p,r I ''^P: 



= VLAA^rip^r + 1 I rip^r + 1) , Cp.^ I Hp^r) = \/ LAA^fu^ \ TLp^r " 1) , 

(P29) 

where the volume factors ensure the correct normalization of the states. Notice 
that this time the propagator is continuous in x — y = /3Z x 0. Moreover, if 
— /3 < 2;o — yo ^ /5 we can write 

(T{ix,.A,,}),^, = ^ e--('^-)^^MMA,.(p) ; (PSD) 



formula ()B[30P can be extended periodically over the whole real axis. To get (IB[30P 
we use that, by Cauchy Theorem: 



1 



-ipo(a;o-J/o). 



1 



po= 



Po + £ "'C' 27r 



^ g-i2(a;o-s/o) ^ 



1 



^2 _|_ p2 2 _ g-i2/3sgn(xo-i;o) ' 



(Pal) 

where the integral in the r.h.s. of (IB13ip is taken counterclockwise along the path 
C' = C+ U C'_ with C^t := G C : 3=2 = ±|p|/2|. The integrals along C± are 
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determined by the residues of the poles in z = ±i|p|; it is easy to see that 



^ -iz{xo-yo). 



1 



1 



which proves (|B[30p . 



-iz/3sgn{a;o-?/o) 



1 

2H 



-\xo-yo\ 



+ 



,(|a:o-yo|-/3)bl 

1 _ ^-m 



(11132) 



B.2 Feynman rules 



The various terms contributing to (|B[5P can be graphicahy represented in terms 
of Feynman graphs; here we shah briefly discuss the rules that one has to follow 
in order to draw and evaluate such graphs. The vertices of the graphs can be of 
two types, see Fig. IB.ll a) they have two external fermionic lines and an arbitrary 
number > 1 of external bosonic (wavy) lines; b) they have four external fermionic 
lines. 




b) 



*— 











3-p, o 



Figure B.l. We can have: a) vertices with two external fermionic lines and an arbitrary 
number to > 1 of external bosonic lines (in this figure to = 3); or 6) vertices with four 
external fermionic lines, joined in pairs connected by a wavy line labelled by at its extrema 
(representing the Coulomb potential (p{x — y, 0)). 



The value of the graph element in item in Fig. IB. H is, depending on the choice 
of the label p: 



r ^ - 
(p = l) 6+ , . TT-ie/ ds[5i]iA. , . 



r=l 
1 



{p = 2) 4,4 Ri^/ ^«NvA.,(.o 

r=l 



.x+sSj ) 



the value of the graph element in item b) in Fig. IB. H is: 



";7^y,o-i,pi^y,o-i,pi¥'(^ ~ 0)'^?/0,2:o^x,o-2,P2^x,cr2,P2 



dEli) 

(11112) 
(lE2l3) 



Notice that in (IB.2I3P we have neglected the factor —1/2 which originally appears 
in the definition of the model, see ()1.2[ip : this can be justified by slightly changing 
the definition of the fermion propagator, see Remark El below. 

Then one considers all possible Feynman graphs, that is all possible ways of 
joining together lines in pairs so that: 



no unpaired line is left over, and 



B.2 Feynman rules 
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22, « 



• only lines (wavy or solid) of the same type can be paired, and 

• in the case of fermionic lines, the orientations of the paired lines have to be 
consistent. 

To each graph we assign the sign (—1)'^ of the permutation vr necessary to bring 
next to each other the pairs of fermionic operators which, in the given graph, are 
paired (one says also contracted), with the ^~ to the left of the associated . 

A solid line connecting two points labelled by xi , X2 with an arrow pointing from 
the point labelled by xi to the point labelled by X2 obtained contracting half-lines 
labelled by cxi , pi and (T2 , P2 corresponds to the fermionic propagator [go-2 ,ai (x2 — 
xi)]^^ Pi Siven by (16122^ : a wavy line connecting two points labelled by xi, ji and 
X2, J2 obtained contracting half- lines labelled by ii and ^2 corresponds to 

Nx2-xi) := [ dsi [ ds2[6j^(g)5j2]i^^i2Wi2,i^{x2,o-xifi,X2-xi+S2dj2-siSj^) . 

Jo Jo 

(IE214) 

To each graph we assign a value, which is the integral over all the space-time labels 
{xj} of the product of: 

• the sign factor, times 

• the product of functions 5'(x2 — xi), ■w{x2 — yi) for every solid or wavy line 
connecting a point labelled by xi to a point labelled by X2, times 

• a factor ibie for every vertex with two incident fermionic half-lines and a 
positive number of emerging wavy lines, where the sign is -|- if the fermionic 
exiting line is of type 1 and it is — otherwise, times 

2 

• a factor —^(p{x2 — xi,0)6ti^t2 for every wavy line labelled by at its extrema. 

Remark 5 Notice that the fermionic propagator that we used to define the Feynman 
rules, given by Eq. /(E[2^) . is not equal to (T{^'~^+})^^ if x — y = (3Z x 0, see 
discussion at the end of subsection \B.1.1[ however, this discrepancy is absorbed by 
the fact that we have not considered the factor — ^ (which is precisely the value of 
the fermionic "tadpole " graph, see Eq. iE[ 24\) ) in the definition of the value of the 
graph element b) in Fig. \B.1\ see l!B.^3\) and 

Given the above rule, it is straightforward to check that the coefficients of the 
perturbative series of (IBISP are equal to those generated by the functional integral 

J P(d*)P(di')e^(*'^) , (IE15) 

where: (i) the interaction is: 

n^,A) := J dx E4^i..i^.^(o/,),.,2(«"^'^'^' - +^^(^) ' 

^^'^^■'l ^..^-W^' "''"'^=^,1$ 



112 



B. Naive perturbation theory 



(ii) the fermionic integration measure has been defined in (|2.2i5|) : (iii) the bosonic 
interaction measure is given by = 1,2): 



r=l,2 



P{dA) :-- 



1 



•exp{-(2MA)"' E Ap[^'(p)"']m4-p} • (IE27) 

Finally, the interaction V{'^,A) is recovered using a Hubbard-Stratonovich transfor- 
mation, namely the fact that: 



[ P{dAo)exp{-{PAA)-' J2 E ie^>p(P)e"'^^"'^''''} ' (E518) 



where the integration measure is defined as: 



n . [(ttba^ 



PGP 



[,l{TTl3AAmp) 



dAo pdA 



p"^0,-p 



0,L 



•exp{-(2/3^A)^i E ^,p'^(p)~'^,-p} • (1E219) 



Pe^;3,L 

Taking the product of this measure with (|B.2I7P we get the full interacting bosonic 
measure (|2.2llip with ^ = 1 ; this concludes the derivation of the functional integral 
representation for the Hamiltonian model. An analogous argument can be repeated 
for the generaing functional of the correlations, but we will not belabor the details 
here. 



Appendix C 

Equivalence of the gauges 



In this Appendix we show that the perturbative series of gauge invariant quantities 
are equal in Coulomb and Feynman gauge; this fact allows us to rewrite the pertur- 
bative series obtained with the Trotter product formula, see Appendix [HI using the 
functional integral representation in the Feynman gauge. To be concrete, we will 
perform the check on the generating functional of the generalized susceptibilities 
introduced in Section I4.3| the same argument can be repeated for the expecta- 
tion value of any gauge invariant operator which can be expressed as series in the 
fermionic and bosonic fields. 

Pick < —h,/3,L < +00, and define 



(i) The fermion field lives on a space-time lattice: the spatial part of the lattice is 
given by the physical honeycomb lattice, while the temporal part is introduced 
by discretizing the time variables with a mesh as discussed in Section 
13. 5t each finite order of the perturbative series of (jC[l[) converges uniformly 
to a limit as — >■ -|-oo. 

(ii) Pyj^ is the bosonic gaussian measure in presence of an infrared cutoff on scale 

and in the ^ gauge; the infrared cutoff is imposed as in (j2.2[9p . 

(iii) The source term D{^, ^,A + J) has been defined in (|4.3[2p . 

The grand-canonical susceptibilities can be obtained by taking derivatives with 
respect to the external fields; our purpose is to show that, order by order in 
perturbation theory, for any finite h, /3 and L, in the limit A^ — >■ -|-oo, the result 
is independent of ^. More precisely, denoting by yvj^^°°^'^'™'($, J) the m-th order 

contribution to w|3^/^°°^'^, we shall prove that 




where: 




J) = . 



(102) 



The proof is easy, and it is based on the following Ward identity: 



(03) 
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C. Equivalence of the gauges 



formula (jCISP can be proved along the lines of Section \'S.5\ see (|3.5ll2p , (|3.5[18|) . 
Let us now show how (|CI3h implies (ICI2|) . First of all, we represent (ICIip in the 
following way: 



m>0 



w;;r''^'"(ci>,j) = y: vaie(a), (04) 

where is the set of connected Feynman graphs of order m generated by (ICiip 
(which have no ^, A external lines), and Valg(^) is the value of the graph Q. Notice 
that for any fixed /3, L and h all the coefficients of the perturbative series of the 
functional integral (|Clip are finite. The photon propagator appearing in these 
graphs is: 

^J^AP'P^) = no = l,n = 0; 

therefore, it is easy to see that: 



mi+m2=m 



V + Ps P 

where / = / ^ and 

(i) F™j^(p) is the sum of all the possible Feynman graphs of order m with two 
bosonic external lines labelled by fi, v and p as external momentum; 

(ii) F^'(p) is the sum of all the possible Feynman graphs of order with one 
bosonic external line and p as external momentum. 

Now, both r^^^(p) and r™(p) can be represented as functional derivatives with 

respect to J^,p, Jv~^ and J/i,p, respectively, of wj^^°°^'^'™(<^, J); and the identity 
(ICISP tells us that if the derivative with respect the field J^^p is contracted with 
then the results is zero. But we see that the /i, v labels in (|C[5P are always 
contracted with at least one among p^, py. therefore, (jClSp is equal to zero. This 
concludes the proof of (IC[2p . 



Appendix D 



The ultraviolet multiscale 
analysis 



In this Appendix we prove Lemma I2.4.1| the proof is based on a multiscale analysis 
similar, but much simpler, to the one discussed in Chapter [3l First of all, we rewrite 
the propagators ^("■^'■^(k), w^'^''"'\p) as 

h=l h=l 

where 

g(^\k) := i/,(A;o)ff("-'^-)(k) , u;W(p) := /7,(po)tf'("-")(p) , ©2) 

with Hi{ko) := x(M-i|fco|), Hhiko) := x(M-'^|A;o|) - x{M~^+^ko\) foil<h<K, 
Hxiko) = xiM-^-'^\ko\) - x{M-^+'^\ko\); note that in the supports of ^("■''■^(k), 
^(«.i,.)(p) it follows that J2h=iHh{ko) = Y.h=iHh{pQ) = 1- Moreover note that, 
defining the norm of a ^4 x A matrix B as ||i?|| := maxj X]i<i<yi some 
positive Co, Ci: 

max||^W(k)|| < CqM-'^ , maxlluje^^fp)!! < CqM-^'* , ©3) 
k " " p 

The first and the third of (IDiSP are obvious; they simply follow from the definitions 
and from the support properties of the single scale propagators. To prove the second 
and the fourth of (ID[3P simply note that w'^^\p) can be rewritten as 

J-aoM (27r) V V^+Pi ' 

and the integrand is bounded proportionally to M~'^^ in the support of Hh{po)- Our 
goal is to compute 

g-/3|A|Fo+V(*(*-''-),A(*-'--)) ^ 

= lim / P(dvI/[i.^l)P(d^[i.^])e^(*'''-'+*'''"'.^'''-^+^'''"') , (P5) 

K^+oo J 
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D. The ultraviolet multiscale analysis 



where P((f^'[^'^^l), P{dA^^'^^) are respectively the fermionic and bosonic Gaussian 
integrations associated with the propagators (jPlip . We perform the integration in 
an iterative fashion, analogous (but much simplier) to the procedure described for 
the infrared integration. We shall inductively prove that 

g-/3|A|Fo,K+VA-(*("'--),A(-'--)) ^ 

= e-mFn J p(rf^[l,/.])p(rf^[l,/.])eV(,'"{*("-)+*[i,.],A('-'-.)+A[i.''l) ^ 

where P{dA^^'''^) are the Gaussian integrations associated with the prop- 

agators Eti^^'^Hk), ELiW^^'Hp) and 

/n m 
[ n ^^^''-t.^T, ,P2j - 1 "^Ki, ,P2 J [ n ■ 

n,m^u <L,P^l^ j=l «=1 
n+r?i>l 

n m 

i=l i=l 

and where the constant and the kernels \nm p p admit bounds analogous 
(see below) to the ones stated in Lemma 12.4.11 In order to inductively prove (|D[6p . 
()D[7p we rewrite 

X J p((i^W)p(dA('^))e^^^'^*'"''''^'^'''''~"^*"''''^*''''''^^"'''~"^^"''^ , &) 

where P{d'if^^^), P{dA^^^) are the Gaussian integrations with propagators ^^'^^(k'), 
w^^\p), respectively. After the integration of A^^^ we define 

= J p((i^W)p(dAW)e^K'(**"''-'+*"'''""+*"'''^'"''-'+^"'''""+^"'') ;(PlO) 
we have that 

- + vi^'-'\^, A) = Y, l^gT^vPi^ + A + AW); n) , Pll) 



n>l 

where 

£l{vPi^ + ¥''\A + A^^^);n) = 

= ^log[ P(d^W)P(d^W)e<''(*+*''''^+^''^) . Pl2) 

As described for the infrared integration, the iterative action of £j^, can be conve- 
niently represented in terms of Gallavotti-Nicolo trees r S TK;h,n, where TK;h,n is 
the set of labelled trees, completely analogous to the set Th,ni unless for the follow- 
ing modifications (see figure IDTT]) : (i) a tree r G TK;h,n has vertices v associated 
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^0 




V 































h h+1 h K+1 

V 

Figure D.l. A tree r G TK;h,n with its scale labels. 



with scale labels h+1 < hv < K+1, while the root r has scale h; (ii) with each 
end-point v we associate + vj/l^'^l, ^C*-''-) + A^^'^^). In terms of these trees, 

the effective potential < h < K (with vj^^ identified with V), can be written 
as 

oo 

-/3|A|F;,+i+vi?^(^[i''^Ut'''^^) = E E VW(r;M/[i''^Ut'''^^) , (Pl3) 

where if vq is the first vertex of r and ti, . . . , {s = Sy^) are the subtrees of r with 
root uo, V'''*)(r; ^[^''^l) is defined recursively as follows: 

• if s > 1, then 

vW(r;^[i''^l,^[i''^l) = (lDll4) 
= l.£T^^(vU^+^)^ri; ^Ifi^M^A^iM). .... v(^+i)(t,; vj/H-'^Uti.'^])^ 

where "pC^+^H^i; ^[I'^^+^l, ylti-'^+^l) is equal to V'^''+^\Ti; ^[L'^+i], At^-'^+i]) if the 
subtree Tj contains more than one end-point, or if it contains one end-point 
but it is not a trivial subtree; it is equal to V{'^^^'^\ A^^'^^) if Tj is a trivial 
subtree; 

• if s = 1 and ri is not a trivial subtree, then v('')(r; ^'I^'''] , A^^''']) is equal 
to £:J^^(V(''+^)(n;^'^'''+^',^f"^'^"^^^)); otherwise, if n is a trivial subtree it 
corresponds to £'^_^^[v{¥^'''+'^l A^^^^+^^) - V{¥^^^1 A^^'^'^] 

Repeating step by step the discussion made for the infrared integration and using 
analogous notations it follows that: 



VW(r;^[i'^U'' 



Val(g) 



E E 

n=iger(r) 

Val(a) = [ n [ n 



(5(t;o)Val(g) , 



/6< 



Val(g) 



ir/ n ^-in^fH n 

V not e.p. " e£v V* e.p. 



(Pl5) 
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D. The ultraviolet multiscale analysis 



the explicit form of the vertex functions K^^"^ can be derived starting from Eq. 
(|2.2115p . Note that if K^"'^ corresponds to a vertex with N external wavy lines then 
its value is bounded as (const. )^A^!~^e^. Using the bounds ()DI3P and proceeding 
as for the infrared bounds it follows that if r G TK-h,N and Q G r(T): 



|Val(g)| < (const.)^e^ TT J-A.r''"^'"~'+<^ 

S 

V not c.p. ^* 

= (const.)^e^^Af-''"o^-o TT J-M-^'^"^^')^" (Pl6) 

V not e.p. 

where Dy := n^—l+J^vyv Therefore, the scaling dimensions are always negative, 
except for the = 1, = for V > V, which corresponds to the fermionic 

tadpole graph; however, in this case the dimensional bound can be improved, noting 
that [^(''^(k)] is odd in /co and that [^(''^(k)] 3_ M^"^^ , which means that: 

1 E [5^'nk)],,, = 0, ^||^5W(k)||<(const.)M-'^ (Pl7) 

Therefore, with any loss of generality we can assume that in ()D[16P n^, > 1. Now, 
let W^^2?fmp/i be the A^-th order contribution to the kernel M^^2»fmp^' ^^^^ 

W^S™.({ka,{p,})= E E Vai(g), (P18) 

reTK;h.N Ger(T) 

l<l=2n 

where the * on the sum denotes the following constraints: UjgpAp(/) = U-!i|Pi, 
U^^pi,k{f) = Uf^iki; Uf^pAfiif) = U^i^i; U^^^^yC/) = U^^^iPi] the bound (|D[T6l) 
implies that 

sup ||^i^2^^,p,^({ki},{pj})|| < 

< (const.)'^e^ E E ^M-^^-o [] i_M-('^«-V)o^ 

|P4l=m « not e.p. 

I<l=2n 

which, after counting the number of Feynman graphs, gives the second bound in 
(12.41 lOp . The bound on the A^-th order contribution to Fh can be proved in the 
same way, and this concludes the proof of (|2.4[lUp . 

To conclude the proof of Lemma [2.4. 1 \ note that the sequences W^^2n^ F^j^ 

are Cauchy in K; in fact, the quantities VF^^^'^^^p,^ - W^^lf^^^^^^^, F^j^ - F^j^, 
with K' > K can be written as a sums over trees with at least one endpoint on scale 
K < h* < K' + \. Therefore, since the scaling dimensions are all negative, by the 
short memory property of the GN trees it follows that the bounds on p tJ.~ 
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^K'2n m p fi^ ^OK ~ ^OK' improved by a factor M with respect to the basic 
bounds (|2.4ll0p : taking the Umit K — >■ +00 (|2.4[11|) follows. This conclude the proof 
of Lemma 12.4.11 



Appendix E 

Proof of Lemma 12.4.2 



In this Appendix we shall prove Lemma I2.4.2t to simplify the notations, we shall 
drop the dependence on the infrared and ultraviolet cutoffs h* , K in the kernels 
(whose presence does not affect the validity of the symmetries). Moreover, we shall 
be concerned only with the case /3 — >■ +00, |A| — )■ +00. The proof is based on the 
fact that and are separately invariant under 

the symmetries listed in Lemma I2.3.H therefore, V^'^ '^](^,^), see (|2.4[7p is also 
invariant under the same symmetries. 

E.l Proof of ( 12.41141) 

Let us rewrite W2,o,p,p'(k' + p'p) as 

W^2,o,py(k' + Vf) = W^2,o,p,p'(Pf) + k;a„Ty2,o,py(PF) + 0(|kf ) . jEIll) 
Constant part. If p = p', by symmetry (8) it follows that 

M^2,o,p,p(k) = -T^2,o,p,p(-A:o, k) 1^2,o,p,p(Pf) = ; (11312) 

\i p ^ (1 by symmetry (4) we get, using that Tp^ = + ^(62 — 26i) (which means 
that pp is left invariant by T), 

M^2,o,P,P'(k) = e'(''-''')'^('^^-'^^)t^2,o,py(A:o, T-i^) ^ 1^2,o,p,p'(Pf) = , dEDS) 
because 1 - e^^r^^h-^^) ^ q. 

Linear part in pQ. By simmetry (6.b) it follows that 

9oVF2,o(k) = doW2fi{ko, fci, -k2) doW2,o{pt) = doW2,o{PF) . (IE314) 
By simmetry (6. a) it follows that 

5oW^2,o,i,i(k) = doW2,o,2,2iko, -ki,k2) 9oW2,o,i,i(p'f) = '9oVF2,o,2,2(p'f) , 

(11315) 

where we used that pp + Pp^ = &i + 621 while \i p ^ p' by symmetry (4) we get that 
(see dEHl) 

= (1 - e'(^-''')^>^^-^^^))aot^2,o,py(PF) (IEI16) 
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E. Proof of Lemma [2X1 



which means that 9oW^2,o,p,p'(p'f) = 0- Finahy, from (5) we find, using (|E.l[4p 

doW2,o{k) = -doW2,o{-^y doW2,o,i,i{PF) = -5oW^2,o,i,i(Pf)* • (107) 
An these properties imply that 9o1^2,o(p'f) = —^zqI with zq reaL 

Linear part in pi. By symmetry (6.b) it foUows that 

9iH^2,o(k) = diW2,o{ko, h, -k2) diW2,o{pt) = 5il^2,o(p^) • dEIiS) 
By symmetry (8) it follows that 

9iVF2,o,p,p(k) = -diW2,o,pA-ko, k) c>iW^2,o,p,p(Pf) = , (11119) 
while from (6. a) we get 

5lT^2,0,l,2(k) = -diW2,0,2,l{ko,-ki,k2) 9iH^2,0,1,2(Pf) = -5l 1^2,0,2,1 (Pf) • 

dEDio) 

Finally, from (5) we find, using (lE.lISp 

5lP^2,0(k) = -diW2,o{-'kr diW2,0,l,2{PF) = -diW2,o,i,2{P%T ■ dEUll) 
All these properties imply that 9iVF2,o(p'f) = — •2iO"2 with zi real and 

^2 = (° . (113112) 

Linear part in p2. By symmetry (6.b) it follows that 

d2W2,o{k) = -d2W2,o{ko, fcl, -k2) d2W2,o{pt) = -92W2,o(p]^) • (13113) 

By symmetry (8) it follows that, 

52W^2,0,p,p(k) = -d2W2,0,p,p{-ko, k) d2W2,0,pAPF) = , dEHW) 

while from (6. a) we get 

92W2,0,l,2(k) = «92VF2,0,2,l(^0, -^1,^:2) ^2 W2,0,1,2(Pf) = ^2 ^2,0,2,1 (Pf) • 

(113115) 

Finally, from (5) we find, using (IE.1I13P 

52VF2,o(k) = -a2T^2,0(-k)* 52VF2,o,1,2(Pf) = 02W2,0,1,2{PfT ■ dOlG) 

All these properties imply that 02^2, oCPf) = —Z2Ujai with Z2 real and 

-1 = j) ■ (1017) 
Finally, the invariance under (4) implies that, using that is left invariant by T: 



Z2 I \-Z2 



which means that zi = Z2- This concludes the proof of (I2.4114P 



E.2 Proof of (|2.4il5p 
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E.2 Proof of ( I2.4L15I ) 



Let us rewrite Wo,2,fi,u{p) as 

Constant part. By symmetry (7), (6. a), (6.b) it follows that if ^ 7^ then 
^o,2,/i,!/(0) = 0; the invariance under (4) implies that 

( Wo,2,l,lW \ / W^o,2,i,i(0) \ 1 ^ 

Wo,2,2,2m) [ 1^0,2,2,2(0) j ' ^> 

which gives M^o,2,i,i(0) = VFo,2, 2,2(0)- The reality condition Wo,2(0)* = M^o,2(0) is 
implied by (5). 

Linear part in pQ. If ^u, z/ 7^ or = 1/ = symmetry (7) implies that 

^oWo,2,^^Ap) = doWo,2,,.A-Po^P) ^ doWoxi^A^) = ; dElS) 
if (/_i, u) = (1,0), (0, 1) by symmetry (6.a) it follows that 

doWo,2,^Ap) = -doWo,2,„APo, -Pi,P2) -doWo,2,^AO) = ; ^4) 
if (/_i, u) = (2, 0), (0, 2) the same is true by symmetry (6.b). 

Linear part in pi. n, u ^ 1 or fi = v = 1 symmetry (6. a) implies that 

^lWo,2,^.AP) = diWo,2,^,APo^ -Pi,P2) 5i 1^0,2,,.,.^ (0) = ; 
if (/i, u) = (0, 1), (1,0) symmetry (7) implies that 

diWo,2,,.Ap) = -diWo,2,^,A-Po,p) diWo,2,,.Ao) = ; dEJe) 

if {n, u) = (1,2), (2, 1) the same is true by simmetry (6.b). 

Linear part in p2. If /i, z/ 7^ 2 of ^ = = 2 symmetry (6.b) implies that 

d2Wo,2,^,Ap) = d2Wo,2,^,APo^Pi^ -P2) d2Wo,2,,,AO) = ; (|El7) 
if (fiju) = (0,2), (2,0) symmetry (7) implies that 

^2WQ,2,^,AP) = -d2Wo,2,^,A-Po,P) 52W^o,2,/.,.(0) = ; (jESlS) 
if (/-i, i^) = (1,2), (2,1) symmetry (6. a) implies that 

d2Wo,2,f.Ap) = -e^'''^'d2Wo,2MPo^ -Pi,P2) ^2Wo,2,^.AO) = . (IE119) 
This concludes the proof of (j2.4[15|) . 
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E. Proof of Lemma [2X1 



E.3 Proof of ( I2.4L16I) 

Case /i = 0. If p 7^ p' by symmetry (8) it follows that 

W^2,i,o,P,P'(k,P) = -^2,i,o,P,P'((-^o,fc),(-Po,l?)) W2xo,p,p'{Pf,0) = ; dESll) 
\1 p = p' by (6. a) it follows that 

^2,1,0,1, i(k,p) = e-'P'^W2,i,o,2,2((A;o,-A;i,A;2),(po,-Pi,P2)) 

^ T^2,i,o,i,i(Pf,0) = W^2,i,o,2,2(Pf,0) , ^2) 

where we used that ff^ + pp^ = 61 + 62- The invariance under (6.b) implies that 

W2,i,o(k,p) = VF2,i,o((A;o,A:i,-A;2),(po,Pi,-P2)) 

^ Vf^2,i,o(PF>0) = M^2,i,o(Pf",0) > 03) 
and, finally, from symmetry (5) we get that, using ()E.3[3P : 

M^2,i,o(k, p) = -H^2,i,o(-k, -p)* ^ T^2,i,o(Pf, 0) = -T^2,i,o(Pf, «)* , (IE3l4) 
which means that VF2,i,o(p'f) 0) ~ "i'^o-^ with Aq real. 

Case p = 1. The invariance under (6.b) implies that 

^^2,i,i(k,p) = M/2,i,i((^o,^i,-^2),(po,Pi,-P2)) 

^ P^2,i,i(Pf,0) = P^2,i,i(Pr.O) • (!E3l5) 

Symmetry (8) implies that, 

W^2,i,i,p,p(k,p) = -W2,i,i,p,p((-^o,^), (-Po,p)) 

^ t^2,i,i,p,p(PF,o) = 0, dElle) 

while symmetry (6. a) implies that, using (lE.SlSp . 

W^2,i,i,i,2(k,p) = -e"'P'^W2,i,i,2,i((A:o,-A;i,A;2),(po,-Pi,P2)) 

^ H^2,i,i,i,2(Pf,0) = -Ty2,i,i,2,i(PF,0) . (|E3l7) 

Finally, symmetry (5) implies that, using (IE.3I5P 

P^2,i,i(k, p) = -T^2,i,i(-k, -p)* ^ W'2,i,i(Pf, 0) = -W^2,i,i(Pf, 0)* , (|E3l8) 

which means that VF2,i,i(p'p) 0) = — Ai(T2 with Ai real and (T2 given by ()E.l[12p . 

Case fi = 2. By symmetry (6.b) it follows that 

^2,i,2(k,p) = -W2,i,2{{ko,ki,-k2),{Po,Pi,-P2)) 

=> P^2,l,2(PF,0) = -T^2,l,2(pi^",0). (|E3l9) 

The invariance under (8) implies that 

W^2,i,2,p,p(k,p) = -W2,i,2,p,p{i-ko,k),{-po,p)) 

W2,i,2,pAPf,0) = 0, dEllio) 



E.4 Proof of (|2.4il7p 
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while symmetry (6. a) implies that 

W^2,l,2,l,2(k,p) = VF2,i,2,2,l((^0,-^l,^2),(P0,-Pl,P2)) 

=> W^2,1,2,1,2(Pf,0) =P^2,1,2,2,i(PF,0) . dEUll) 

Symmetry (5) implies that, using (IE.3l9|) . 

W2,i,2{K P) = -W^2,i,2(-k, -p)* P^2,1,2(Pf, 0) = TV2,1,2(Pf, 0)* , (|E3ll2) 

which means that W2,i,2{p'f,0) = —X2^o'i with A2 real and fJi given by (lE.lllTp . 
Finally, the invariance under (4) implies that, using that is left invariant by T: 

which means that Ai = A2. This concludes the proof of ()2.4ll6p . 

E.4 Proof of (I2.4L17I) 

If /i = (/ii, /X2, /is) has an odd number of labels fii = then by (8) it follows that 
Wo,3Ap,p') = -Wo,3A{-po,P), {-Po,P)) ^o,3,m(0,0) = ; (|Ell) 

the same is true if the number of labels //j = is even by (7). This concludes the 
proof of (j2.4ll7l) and of Lemma EMI 



Appendix F 

Lowest order computations 



In this Appendix we reproduce the details of ah the lowest order computations 
performed in this Thesis. 

F.l The Beta function of the wave function renormal- 
ization 

We start by performing the computation of the lowest order contribution in the 

z (2) 

running coupling constants to the Beta function of the wave function renor- 

mahzation; by definition, see (i:i2ll5p and WTI^ . = zo,h = r° ^qVF's^J^^^^IO). 

Following the rules discussed in Section I3.2| its lowest order contribution ji^ is 
graphically represented by the Feynman graph depicted in Fig. IF.li 




Figure F.l. The lowest order contribution to /?^. The wavy line labelled by a "9o" symbol 
corresponds to 9ow^''^Hp); a sum over all the possible scales /ii, /i2 > /i + 1 is understood. 

Therefore, we find (repeated [i labels are summed): 

-r°et.+2^ / ^ [5o^('^+^)(p)]r^5('^+^)(p)r^ . 

Using inductively the Beta function equations for f/i+i, .^h+ii ^fi,h+i "we find that. 
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F. Lowest order computations 



setting Fh+i := f^_^-^ + 2/^+1/^+2, 

Ph - ^f.,hj (2vr)3 2|p|3 pI + vI\p\^ ^ ' 

+0(e^) + OCe^M^^) ; 
passing to radial coordinates p = p(cos 9, sin 9 cos 93, sin sin if) , using that 

/ dp U'h+ih+i + fUih+2 + fU2h+i) = , e^^.r" r^^r" ri; = eg,, - 2vl_A,h 

we get: 



1 COS2 

a COS 9 



cos^ ^ + ^^^^ ' 



+ 



+0{eHr'') + O(e^) . (103) 



The integral over the radial coordinate p can be computed by using the definition 

dp , 
/o P 



r^(4\i + 2A+iA+2)= (104) 

JO P 



'^^-[2(X(P) - x(Mp)) - ixHp) - xHMp))] = r ^ = logM 



/op Ji p 

finally, an explicit evaluation of the integral over d cos 9 leads to (|3.4[19|) . 

F.2 The Beta function of the Fermi velocity 

By definition, see formulas ()3.2[15l) and (I3.4l2p . := z^^\ — VhZ^j^, with zi^h = 

r^(9i W2 Q^j^ j^(0). At second order, proceeding as in the derivation of ()F.l[2p . we 
find: 

^1 = -l^ I ^^^irflf ^^^^^^^^ + + ^(^'"^"^ ' 

(the term involving the derivatives /^_,_|, /^+2 support functions is zero, as 

before); using that 

passing to radial coordinates, we get that: 



(2) 2 



1 



167r2 



f°°dp,,2 ^l r /"^ J /I sin^ 

— (4+1 + 2A+iA+2)J • [y dcos 



cos^ 9 + v\ sin^ 1 



The integral over the radial coordinate can be evaluated as in (lF.ll4p . and gives 
logM; therefore, an explicit evaluation of the integral over dcos 9 leads to 



(2) _ 2 - 1 log M ( ar ctan 4 4 - ar ctan 4 
\H-%,H^h 8^2 I ^3 

which, combined with /3^'''^^ = — Vhzj^h, proves (j3.4[19|) . 



(IO2) 



F.3 The Beta function of the Kekule mass 
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F.3 The Beta function of the Kekule mass 



By definition, A;,^;, = i-f^^'^y CaViW!i^^{\^), that is, setting F/,+i(p) := /^^^(p) + 
2A+i(p)A+2(p): 



.(2) 



dp Fh+i{p) 



(27r)3 2|p| pI + vIW 
(^2 , 27;2p2 N f dp Fh+i{p) l_ 



2|p| Po + vl\p\' 



(Di) 



1 

'8^ 







p 



dcos 9- 



1 



cos^ + vf^ sin 



where in the second equality we used that 7^7^-''°^* = — 7^-"'H/i, and that 7(j'o)*^(io) = 
I. The integral over the p coordinated can be computed as in (IF.l[4p , and an explicit 
evaluation of the integral over d cos 6 leads to (|5.3l23p . 



F.4 Anomalous exponents of the generalized suscepti- 
bilities 

In this Appendix we shall compute the lowest order contribution to the Beta func- 
tions /3" /j ]^ of the running coupling constants ^, whose flow has been discussed in 
Section 14 3.11 in particular, these computations will allow us to compute the lowest 
order contributions to the anomalous exponents of the generalized susceptibilities. 
Basically, we will compute the Feynman graph depicted in Fig. IF. 21 for all possible 
values of the inner scale labels, a, uji and 002, and with the external field h. 
does not depend on the bond label a, therefore we set a = 1 for simplicity. 




Figure F.2. The lowest order eontribution to a sum over the inner scales hi > h is 

understood. 



Preliminaries. Before working out explicitly the computations, it is useful to 
recall a list of notations and identities that will be repeatedly used in the following. 
The main difference with respect to the previous lowest order computations is that 
here we will not use the compact relativistic notations. 
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F. Lowest order computations 



(i) The vertex on scale h from which emerges the wavy hne associated to the 
external field is represented as: 



where the symbol "^^^^ has been introduced in (|3.2i4p . and the matrices 
have been defined in ()4.3[7p . 



(ii) The A^^^ vertex is represented as: 

where eo,h = eo,^ and ei^h = (i2,h = Vh-i^i^h, and the matrices T^ are given by 

(iii) The fermion propagator on scale h for k = k' + p"^ and |k'| small is given by 

^^.'^k') = .^]^^'Lr.. i - r^^o + VHi^ri) (1 + o(k)) . (IE14) 



A;^^ + ?;(k')2|fc'|2 

(iv) The following identities will be useful: 

(v) Setting F^IpI) := A(p)3 + 3h{p)^h+i{p) + ^fh{p)fh+i{p?, it follows that 

= logM . Me) 

P 

The first part of the calculation is the same for all a, cv; in fact, we can write 
the value of the graph depicted in Fig. IF.2l as follows, neglecting corrections O(e^) 
and 0{e'^M^) (repeated a indeces are not summed): 

paoa,(2) _ 

2 f JP_rM n(hi)(^)Y'-a{h2)r w^^^'Hv) = I — ^^^'^'^ • 



-1 



■ log Mrii^r::^r-r;:^r^^^ + o{e'') + oCe^M'^) . dElr) 



127r^ 



F.5 Lowest order check of the Ward identities 
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Using (|F.4[5P we get that: 

-p/x -pi/ -p;/ -p/i pi?_ pA* pi/ pi?„pl' p/J 

pM pi^ p^z+p!^ pM 9r^+ pM pf p^'+pt' pM p£^+ 

p/x pi/ pCDWpi^ pM _ pCDVF p/J pv pCDWpj/ pM _ gpCDVy 

pMpi^pDpv pA^ ^ pD ^ r![r:;r^r:;r![ = ; dEls) 

plugging these identities into (IF.4I7P we get the lowest order contribution to 
f3^ hi', the results listed in Table ITT] follow from r/2 := limft_^_oo log^/ + hi)- 



F.5 Lowest order check of the Ward identities 

In this Section we shall check at lowest order in non-renormalized perturbation 
theory the validity of the Ward identities that allowed us to prove the infrared 
stability of the flows of the effective charge and of the photon mass. 



F.5.1 Ward identity for the photon mass 

We start with the check of the Ward identity for the photon mass; as we have 
discussed in Section I3.5[ the gauge invariance of the theory implies that the effective 
photon mass is vanishing, i.e. that the electromagnetic interaction is unscreened. 

At lowest order, the graphs contributing to the dressed photon mass are depicted 
in Fig. IF. 31 for fi = or v = the second graph is absent. As we are going to show 
here, the sum of the two graphs is exactly vanishing. 




Figure F.3. The lowest order contribution to the rcnormalizcd photon mass in naive 
perturbation theory. The graphs are all computed at vanishing external momenta, and for 
/.J = or = the second graph is absent. The sum of the graphs is exactly vanishing. 

These graphs are obtained by expanding the interaction V{'^,A) to second order 
in A; it follows that, see Eq. (|2.2[14|) . assuming without any loss of generality that 
the field A^^p is supported inside the first Brillouin zone (p = in the graphs 
represented in Fig. IF.Sp : 
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F. Lowest order computations 



where eo = e, ei = 62 = ve with v = |t, ei := (1,0), 6*2 = (0, 1), and 

Notice that r^(^,0) = T^, where are the matrices introduced in Eq. (|2.4[14|) . 
Before starting, it is useful to recall that the fermion propagator is given by 

+ v^\n{k)\^ \-vn{k) iko J 

with Vl{k) = I X]j=i,2,3 With these notations, the value of the sum of the 

two graphs in Fig. ()F.3p is: 

~e^e, j ^ Tv{^(k)r''(^, 0)5(k)r^(^, 0)} + dEIls) 
-2e\ j '^T,[g{\i)T^'-{k,Q,Q)} . 

We shall consider only the cases fi = 0, fi = 2; the case /i = 1 can be obtained from 
the case /u = 2 using the symmetry (4) of Lemma 12.3.11 

Case ^ = 0. Let us first consider the case ^ = 0; using that 

ff(k)rO(fe, o)ff(k) = -dom , (IEI16) 

we get that (jF.SiSP is equal to: 

ee,J ^Tr{5o5(k)r'^(fc,0)} = . ^7) 

Cases n = 2, v ^ 0. Since 

m^^\k, o)ff(k) = -d2m , dEiis) 

we can rewrite the integral in the first line of (IF.5[5P as: 



I ^Tr{5(k)a2r'^(^,0)}; dEHQ) 



using the identity 



r^'2(fc, 0, 0) = -^52r'^(A:, 0) , dEHio) 

we see that the second line of (1F.5I5P exacly cancels the first line. 

The case ^ = v = 1 can be obtained from // = 1/ = 2 by using symmetry (4) of 
Lemma [2. 3.11 it follows that even in this case the sum of the two graphs is vanishing. 
All the other cases can be obtained using the cyclicity of the trace; this concludes 
the lowest order check of the Ward identity for the photon mass. 



F.5 Lowest order check of the Ward identities 
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F.5. 2 Ward identity for the effective charge 

In this Section we shall check at lowest order in non-renormalized perturbation 
theory the Ward identity for the effective charge; this identity allowed us to prove 
at all orders in renormalized perturbation theory that the effective charge is close 
to the bare one. 

Graphically, the lowest order contribution to the renormalized charge in naive 
perturbation theory is given by Fig. IF.4i 




Figure F.4. The lowest order contribution to the renormalized charge. To avoid ambigui- 
ties, the photon propagator is regularized in the infrared by a cutoff function. The sum of 
the graphs is exactly vanishing, and the cancellation does not depend on the presence of 
the cutoff. 



notice that the first and fifth graph are logarithmically divergent; therefore, to 
avoid ambiguities we fix an infrared cutoff on the photon propagator, for instance 
by replacing w{p) with w^-^^p) defined in (13.5131) with e = 0. As we are going to 
show, the sum of these six graphs is exactly vanishing; therefore, the dressed charge 
is equal to the bare one at lowest order. Remarkably, this cancellation does not 
depend on the presence of the infrared cutoff on the photons; this fact has been 
exploited in Section 13.51 to derive a Ward identity for the effective charge on a given 
scale h. We shall only consider the cases /i = and /i = 2; the renormalization of 
the charge corresponding to /x = f is equal to the ^ = 2 one by symmetry (4) of 
Lemma 12.3.11 

The graphs in Fig. IF. 41 can be obtained by expanding the interaction V{'^,A) 
to third order in A; it follows that, calling V2{^,A) the expansion to second order 
derived in Eq. (IF.5[ip : 

V{^,A) = V2{^,A)+ dESlll) 

with 



Case /u = 0. In this case only the first and fifth graphs in ()F.4p are present; and 
using the identity (IF.5i6p we get that their sum is exactly vanishing. 
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F. Lowest order computations 



Case /X = 2. Here the situation is more complicated, because of the presence of 
the second, third, fourth and sixth graph in Fig. IF.41 and because the 82 derivative 
can act on the vertex functions. However: 

(i) the first graph is cancelled by part of the fifth, because of the identity (IF.5I9P : 

(ii) the contribution due to the action of the 82 derivative on the vertices of the 
fourth graph are exactly cancelled by the second and third graphs; this is so 
because of the identity 

r^'2(fc,p,o) = -1 ArM(^,p-) , dEllis) 

Z Ok2 

where the factor ^ is compensated by a factor 2 which counts in how many 
ways we can choose the external wavy line in the second graph; 

(iii) the sum of the fouth and sixth graphs in Fig. IF. 41 is exactly zero, because of 
the identity 

3 ok2 

where the factor | is compensated by a factor 3 which counts in how many ways 
we can choose the external wavy line in the third graph. 

The case = 1 can be obtained from the case = 2 by symmetry (4) of Lemma 
12. 3. H it follows that even in this case the sum of the graphs in Fig. IF. 41 is exactly 
vanishing. This concludes our lowest order check of the Ward identities. 



F.6 The relativistic propagator in real space 

In this Appendix we reproduce the details of the computation leading to (|4.3i48p . 
Let us rewrite ffff^Cx) as 



:]:(x) =e-'P^.-/ ^^-■^'''''^-^{-Tlko + kr.) . dEll) 



Notice that 



D |k'|2 "^"^ ~ 'dxa J D |k'|2 

. d [dp 



J ^sin0d0/;,(p)e-^''l-l-^^ , (IE612) 



dxa J B 

where to get the second line we have chosen 9 to be the angle of k' with respect to 
X. Then, performing the integral over dcosO and the derivative with respect to Xa 
we get that 



(jE6l3) 



F.7 Proof of ()4.3i49p 
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Plugging this result in ()F.6[ip and using that J2h<o fhi') = '^^ g^t: 



h<0 



(5i/2)|x| 



di^-^x(j^) + I dt sintx'(A) 



t 



X 



(lE6l5) 



finally, using integration by parts and Riemann - Lebesgue Lemma, we know that 
for any > 



lim |x|^ 
|x|— s>+oo 



^^"*^xf A) - ij = , lim Ixl^ / sintx'fAl = , 

^X/ J |x|^+oo J ^X/ 



|x|^>+oo 



|x| 



that is, we can rewrite g^^c{'^) a-s: 



^"'^^^ (i?i/2)|x|3 7 t + Vl + |x|^ 



dESle) 
dEHlT) 



e -^Pf^_Z1^1±L1_111^ + of L 

V 1 + X 



6'|xl3 



N 



b' :-- 



Bi Sir 



F.7 Proof of (I4.3L49I) 

In this Appendix we shall compute the trace 

Tr{r"5^,,£(x)r°<7^,,£(-x)} , dEHi) 

for a = E±,CDW, D. First of ah, using (I4.3l48p . we rewrite (IRGITI) as: 

y2, |6 Tr{r"(-xor°^ +x,cjr°(xor°^ -x.rLj} ; (IEI2) 



then, we use that: 

-^^1^^ ^Wi^ -^^2] 

Tr{r^±ro^r^±roJ 



^^1^^ ^wi^ ^^2/ 

^i^2Tv{r^^^r2^r^^^r2 J = -2 , (1E14) 
±Tr{r^±ri^r^±riJ 

Ta;ia;2Tv{r^±r2 ^r^±r2^ } = ^2 , (|El5) 



and that Trjr^r/^^r^r;;; J = for Va. From (I4.3[47|) . (^2]) - (If:6T5]) it is 

easy to see that the results claimed in (|4.3149p follow. 



Appendix G 

Corrections to the Ward 
Identities 



In this Appendix we wiU show that the corrections to the Ward identities for the 
photon mass (|3.5l24p and for the effective charge (|3.5129p are exponentially vanishing 
in the ultraviolet limit K — >■ +00. The corrections p^A||^2^^j^(p), p^Ag'^^^j ^(k, p) 



can be obtained as functional derivatives with respect to the external fields Jp, Ji/-p 

or Jp, 4>t+p,a,i^ 

J, cP) := log J P(dM/)P>/,(dA)e^(*'^+^)+^(*''^)+'^(*'^") , ([Gil) 



where 

and given by: 

-ie-'C^p,{k,p) := (PS) 
e-im(p'-i)(^-i(fco +p,) _ l)[i?(k + p)]^^, -e-^^(^-i)(x^Hfco) - l)[i?(k)]^^, 

As for the generating functional of the correlation functions, the functional integral 
(jGlip can be studied using multiscale analysis and renormalization group; the only 
difference is the presence of the term C{^,J). A crucial role in the analysis is playied 
by the properties of the function C'^(k, p); in fact, it is straightforward to check 
that 

5(k + p)C^^(k,p)^(k), (IG14) 

is non- vanishing only if |A;o +po| ^ , and/or |A;o| > . Moreover, when it is 
non- vanishing (IG[4p it is dimensionally bounded as (const.) |p|M~^-'^. Let us choose 
the external fermionic field as in Section I4.2t proceeding as in Appendix |Dl after 
the integration of the ultraviolet degrees of freedom we are left with: 

p(^(<O))p^^(^^(<O))gV(O)(*(<o),A(<o)+J)+i?(*(<o),</,)+C(o)(M>(<o),0,J,A(<o)+J) ^ 

where C(°) contain all the J-dependent terms, and C(°)(0, J, 0) := 0. The 

kernels of C(^°), C(°) admit a graphical representation in terms of Gallavotti-Nicolo 
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G. Corrections to the Ward Identities 



trees, following the rules discussed in Section [3.3.21 and Appendix iPj with the only 
difference that new endpoints appear, corresponding to Jp^j^_,_p o-p^/^'C^i P)^k o-p'! 
as usual, the trees can be expanded in terms of Feynman graphs, with the new 
endpoints corresponding to the vertex represented in Fig. IG.ll 




Figure G.l. The vertex corresponding to Jp^k_|_p ^C^(k, p)^^ ^. 

Because of the properties of (IG[4p , at least one of the solid lines emerging from 
the new vertex has to be contracted on scale K; therefore, using the general dimen- 
sional bound ()D[16P and the short memory property of the GN trees, the kernels of 
are bounded proportionally to (const.) I p|M ^^Y[i=o-^ '^''^ , where is the 
number of (p external lines and {kj} are their momenta, while those in C^^^ admit an 
improvement of a factor |p|M~^^ with respect to their basic dimensional bounds. 
The bound on the kernels of C^-"^ is derived from ()D[16P taking into account that 
the tadpole graph associated to the new vertex depicted in Fig. IG.ll is exactly 
vanishing. 

The single scale integration is performed as in Section 14.21 with the only dif- 
ference that the J-dependent terms independent of A^-~^^ + J are taken 
into account in the definition of C^^" J) := C^-^\(l), J) + C^^\(l), J), where C(°) 
collect the new J-dependent kernels independent of A^-~^^ + J produced 
in the integration of the scale 0, while all the other J-dependent ones are contained 
into C^^^\^^-~^\ (j), J, A^-^^^ + J), and so on. After the integration of the first 
\k\ < \h\ scales we are left with: 

where all the J-dependent terms absorbed in C^-^'^ = C^-^^ ~^J2w=k+i ^^^'^ C^^\ 
with (:('')(0,(/-, J,0) := 0, while all the other objects appearing in (jGlGl) have been 
defined in Section 14.21 Again, because of the short memory property, the dimen- 
sional bounds on the kernels of C^'^^ and of C^^^ admit respectively an improvement 
of a factor \p\M^^~-^~^^\ \p\M^^~^~^^\ with respect to the basic ones derived in 
Section [3.3.21 Finally, the integration of the scales < h proceeds in a way completely 
analogous to the one sketched in Section 14.21 The conclusion is that the kernels of 
C^"') are exponentially vanishing for K ^ +oo, uniformly in w; this implies that 

hm A[';5,,(k,p)= hm A[';i5(k,p)=0, (07) 



as desired. 



Appendix H 

Proof of ^TT^UT) 



In this Appendix we prove that in presence of Kekule distortion the effective poten- 
tial arising after the integration of the ultraviolet degrees of freedom has the form 
()5.3[6p . With respect to the case Aq = 0, the only new check to perform is to show 
that the relevant terms, that is the mass terms for the fermion and photon field, 
have the form claimed in (j5.3l4|) . (|5.316|) . The symmetry properties to which we 
shall refer here are the (1) - (8) of Lemma 15.3.11 

H.l Fermionic mass terms 

By symmetry (8) it follows that 

2,0,p,p,uj_ 



Symmetry (4) implies that 

^^S,p',.,.'(o) = (IH312) 

= ^^2Ap,pI-,-'(0) {m^2 - A){P - 1) - ■^rF{^2 - A){P' - 1)} , 

while symmetry (6.b) implies that 

therefore, from the combination (4) - (6.b) - (4) - (6.b) we find that: 

<!p,p',.,JO) = e^^^^-^'^^''^^^^^-^'^^^Sp,p',.,.(0) . (104) 
Formula ()H.1[4[) together with (IH.1I2P gives, if p / /)': 

<J,P,p',.,.(0) = = W^iJip, .,J0) . (|H315) 
Symmetry (4) implies in particular that 

w^a,2,.,- JO) = wifSt,-^"'^^ ' me) 

while symmetry (5) implies that 

^Si2,.,- JO) = JO)] * ; (jH3l7) 
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H. Proof of l5.3i6l 



finally, symmetry (6. a) implies that 



Therefore, setting i^'iiUs-P.o; -<^(0) =• ^o^^ formulas (iHlTHl) . dlLTTSl) imply that: 

moreover, formula (IH.l[7p implies that a[,^^ G R. This, together with p.llip . 
(IH.ll5jl . concludes the proof of (fOTD . 

H.2 Bosonic mass terms 

From (6. a) and (7) it follows that, if ^ 7^ z^: 

<2,k.(o) = , (Di) 

while from (4) it follows that 

^o'2,o,o(0) = W^lofiW = W^lofiW ■ (12212) 
Symmetry (6.b) implies 

and from (4) we get 

therefore, formulas p.2[2p - (IH.2I4P imply that: 

<2il(0) = <?2,2(0). (13215) 

Formulas (|H.2[ip , (1H.2I5P prove the structure of the bosonic mass terms claimed in 
()5.316p : therefore, this together with what has been discussed in the previous Section 
concludes the proof of 15.3161 
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